LPNll-94 
IFIC/11-72 
ZU-TH 27/11 



O 



^3 



Space- like (vs. time-like) collinear limits in QCD: 
is factorization violated? 



Stefano Catani Daniel de Florian and German Rodrigo ^'^'^ 



•^"^INFN, Sezione di Firenze and Dipartimento di Fisica e Astronomia, 



OO . Universita di Firenze, 1-50019 Sesto Fiorentino, Florence, Italy 

(N- 

^^•'Departamento de Fisica and IFIBA, FCEYN, Universidad de Buenos Aires, 
^ ■ (1428) Pabellon 1 Ciudad Universitaria, Capital Federal, Argentina 

^ '^'^^Institut fiir Theoretische Physik, Universitat Ziirich, CH-8057 Ziirich, Switzerland 

^ ■ ^'^^Instituto de Fisica Corpuscular, UVEG - Consejo Superior de Invest igaciones 

Cientificas, Pare Cientific, E-46980 Paterna (Valencia), Spain 

(N 

>: 

^ : Abstract 

I We consider the singular behaviour of QCD scattering amplitudes in kinematical 

^ ■ configurations where two or more momenta of the external partons become 

collinear. At the tree level, this behaviour is known to be controlled by 
factorization formulae in which the singular collinear factor is universal (process 
independent). We show that this strict (process- independent) factorization is 
^ I not valid at one-loop and higher-loop orders in the case of the collinear limit 

■ in space-like regions (e.g., collinear radiation from initial-state partons). We 

introduce a generalized version of all-order collinear factorization, in which the 
space-like singular factors retain some dependence on the momentum and colour 
charge of the non-coUinear partons. We present explicit results on one-loop 
and two-loop amplitudes for both the two-parton and multiparton collinear 
limits. At the level of squared amplitudes and, more generally, cross sections 
in hadron-hadron collisions, the violation of strict collinear factorization has 
implications on the non-abelian structure of logarithmically-enhanced terms in 
perturbative calculations (starting from the next-to-next-to-leading order) and 
on various factorization issues of mass singularities (starting from the next-to- 
next-to- next-to-leading order) . 
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1 Introduction 



A relevant topic in QCD and, more generally, gauge field theories is the structure of the 
perturbative scattering amplitudes in various infrared (soft and collinear) regions. Virtual 
partonic fluctuations (i.e. partons circulating in loops) in the infrared (IR) region lead 
to divergent contributions to scattering amplitudes in four space-time dimensions. Real 
radiation of soft and collinear partons produces kinematical singularities, which also lead 
to IR divergent contributions after integration over the phase space of the emitted partons. 

In the context of dimensional regularization, the (virtual) IR divergences of QCD am- 
phtudes have been studied [HEllSlllllSlEllTllH] at one-loop, two-loop and higher-loop 
orders. The singularity structure related to (single and multiple) soft-parton radiation has 
been explicitly worked out [HI EDI [HI [121 [131 HI] in the cases of tree-level, one-loop and 
two-loop scattering amplitudes. 

These studies and the ensuing understanding of (virtual) IR divergences and (real) 
soft-parton singularities rely on general factorization properties of QCD. The structure 
of the IR divergences and of the soft-parton singularities is described by corresponding 
factorization formulae. The divergent or singular behaviour is captured by factors that 
have a high degree of universality or, equivalently, a minimal process dependence (i.e. a 
minimal dependence on the specific scattering amplitude). To be precise, the divergent 
or singular factors depend on the momenta and quantum numbers (flavour, colour) of the 
external QCD partons in the scattering amplitude, while the detailed internal structure of 
the scattering amplitude plays no active role. Of course, in the case of soft-parton radiation, 
the singular factors also depend on the momenta and quantum numbers of the emitted soft 
partons. 

In this paper we deal with coUinear-parton singularities. The singular behaviour of QCD 
amplitudes, in kinematical configurations where two or more external-parton momenta 
become collinear, is also described by factorization formulae. 

Considering the case of two collinear partons at the tree level, the collinear-factorization 
formula for QCD squared amplitudes was derived in a celebrated paper [15]. The corre- 
sponding factorization for QCD amplitudes (rather than squared amplitudes) was intro- 
duced in Refs. [16], [17]. At the tree level, the multiple collinear limit of three, four or more 
partons has been studied [HI HHl [HI |20l [21], |22j for both amplitudes and squared ampli- 
tudes. In the case of one-loop QCD amplitudes, collinear factorization was introduced in 
Refs. [231 [211 [121 [2S], by explicitly treating the collinear limit of two partons. Explicit, 
though partial, results for the triple collinear limit of one-loop amplitudes were presented 
in Ref. [26]. The two-parton collinear limit of two-loop amplitudes was explicitly computed 
in Refs. [271 [H]. The structure of collinear factorization of higher- loop amplitudes was 
discussed in Ref. [25] . 

The collinear-factorization formulae that we have just recalled are similar to the factor- 
ization formulae that apply to virtual IR divergences and soft-parton singularities. How- 
ever, the collinear-factorization formulae are 'more universal'. Indeed, the collinear singular 
factors only depend on the momenta and quantum numbers (flavour, colour, spin) of the 
collinear partons. In other words, the collinear singular factors have no dependence on the 



1 



external non-collinear partons of the QCD amplitudes. Throughout this paper, this feature 
of collinear-parton factorization is denoted 'strict' coUinear factorization. 

Despite so many established results, in this paper we show that strict coUinear factor- 
ization of QCD amplitudes is actually not valid beyond the tree level. 

We are not going to show that some of the known results at one-loop, two-loop or higher- 
loop levels are not correct. We simply start from the observation that these results refer 
(either explicitly [26l EZl E] or implicitly) to the coUinear limit in a specific kinematical 
configuration. This is the configuration where all partons with coUinear momenta are 
produced in the final state of the physical process that is described by the QCD amplitude. 
We refer to this configuration as the time-like (TL) coUinear limit. 

In the TL coUinear limit, strict coUinear factorization is valid. In all the other kine- 
matical configurations, generically denoted as space-like (SL) coUinear limits, we find that 
strict coUinear factorization is not valid (modulo some exceptional cases) beyond the tree 
level. We also show that, in the SL coUinear limits, QCD amplitudes fulfill generalized 
factorization formulae, in which the coUinear singular factors retain some dependence on 
the momenta and quantum numbers of the external non-collinear partons of the scattering 
amplitude. 

The violation of strict coUinear factorization is due to long-range (gauge) loop interac- 
tions between the coUinear and non-collinear partons. These virtual radiative corrections 
produce absorptive contributions that, due to causality, distinguish initial-state from final- 
state interactions. In the TL coUinear region, all the coUinear partons are produced in 
the final state and strict factorization is recovered because of QCD colour coherence (i.e., 
the coherent action of the system of coUinear partons). The SL coUinear region involves 
coUinear partons in both the initial and final states and, therefore, causality limits the 
factorization power of colour coherence. 

Owing to their absorptive ('imaginary') origin, strict- factorization breaking effects partly 
cancel at the level of squared amplitudes and, hence, in order-by-order perturbative calcu- 
lations of physical observables. Indeed, we find that such a cancellation is complete up 
to the next-to-leading order (NLO). Nonetheless, strict factorization is violated at higher 
orders. For instance, the simplest subprocess in which strict coUinear factorization is defi- 
nitely violated at the squared amplitude level is 2 — > 3 parton scattering, in the kinematical 
configurations where one of the three final-state partons is coUinear or almost coUinear to 
one of the two initial-state partons. In this subprocess we find non-abelian factorization 
breaking effects that first occur at the two-loop level. Therefore, these effects contribute 
to hard-scattering processes in hadron-hadron collisions: they produce next-to-next-to- 
leading order (NNLO) logarithmic contributions to three jet production with one low-pr 
jet (the low-pT jet is originated by the final-state parton that is almost coUinear to one of 
initial-state partons), and next-to-next-to-next-to- leading order (N^LO) contributions to 
one-jet and di-jet inclusive production. 

The strict factorization breaking effects uncovered in the simple example of 2 — i- 3 
parton scattering have more general implications in the context of perturbative QCD com- 
putations of jet and hadron production in hadron-hadron collisions. Starting from the 
N^LO in perturbation theory, these effects severely complicate the mechanism of cancel- 
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lation of IR divergences that leads to the factorization theorem of mass (colhnear) singu- 
larities [29]. These complications challenge the universal (process- independent) validity of 
mass-singularity factorization, and they are related to issues that arise in the context of 
factorization of transverse- momentum dependent distributions [301 ED E2] • The perturba- 
tive resummation of large logarithmic terms produced by colhnear parton evolution is also 
affected by the violation of strict colhnear factorization: parton evolution gets tangled with 
the colour and kinematical structure of the hard-scattering subprocess, and this leads to the 
appearance of 'entangled logarithms'. An example of entangled logarithms is represented 
by the class of 'super-leading' non-global logarithms discovered [33j in the N^LO computa- 
tion of the dijet cross section with a large rapidity gap between the two jets. Indeed, the 
physical mechanism that produces those super-leading logarithms [3l] is directly related to 
the mechanism that generates the violation of strict colhnear factorization. 

The outline of the paper is as follows. Sections [2H1] are devoted to the two-parton 
colhnear limit. In Sect. |2l we consider tree- level amplitudes; we introduce our notation 
and, in particular, the colour space formulation based on the colhnear splitting matrix. In 
Sect. [31 we review the known results on the TL colhnear limit of one- loop amplitudes. The 
SL colhnear limit at one- loop level is considered in Sect. jH Here, we illustrate the violation 
of strict factorization, we introduce our generalized form of colhnear factorization, and we 
present the result of the one-loop splitting matrix to all-orders in the dimensional regular- 
ization parameter e. In Sect. [5l the study of the colhnear behaviour of QCD amplitudes is 
extended to the multiparton colhnear limit and beyond the one-loop level. In particular, 
we illustrate the violation of strict colhnear factorization by deriving the explicit expression 
of the IR divergences (i.e., the e poles) of the one- loop multiparton splitting matrix. In 
Sect. [6l we consider the all-order IR structure of the colhnear splitting matrix, we present 
the explicit IR divergent terms at the two-loop level, and we discuss the ensuing new fea- 
tures of strict colhnear factorization. In Sect. [71 we use our results on the colhnear splitting 
matrix to compute the singular colhnear behaviour of squared amplitudes. We explicitly 
show that strict colhnear factorization is violated also at the squared amplitude level, and 
we comment on the implications for QCD calculations of hard-scattering cross sections in 
hadron-hadron collisions. In Sect. [H we briefly summarize the main results. Additional 
technical details are presented in the Appendices. In Appendix \^ we illustrate the vi- 
olation of strict colhnear factorization within the (colour-stripped) formulation in terms 
of colour subamplitudes and splitting amplitudes. In Appendix [Bl we explicitly compute 
the IR divergences of the two-loop splitting matrix. In Appendix [Cl we discuss how strict 
colhnear factorization is recovered in the TL colhnear region. 



2 Collinear limit and tree-level amplitudes 

We consider a generic scattering process that involves external QCD partons (gluons and 
masslesf] quarks and antiquarks) and, possibly, additional non-QCD particles (e.g. partons 
with no colour such as leptons, photons, electroweak vector bosons, Higgs bosons and so 
forth). The corresponding ^-matrix element (i.e., the on-shell scattering amplitude) is 
denoted hj J^{pi,p2, ■ ■ ■ ,Pn), where Pi {i = 1, ... ,n) is the momentum of the QCD parton 

*The case of external massive quarks and antiquarks is not considered in this paper. 
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Ai {Ai = g,q or q), while the dependence on the momenta of additional colourless particles 
is always understood. 

The external QCD partons are on-shell (p? = 0) and with physical spin polarizations 
(thus, Ai includes the corresponding spin wave functions). Note, however, that we always 
define the external momenta pj's as outgoing momenta. In particular, the time-component 
(i.e. the 'energy') of the momentum vector (/i = 0, 1, . . . , (i — 1) in ci space-time dimen- 
sions is not positive definite. Different types of physical processes with n external partons 
are described by applying crossing symmetry to the same matrix element A4{pi,p2, ■ ■ ■ ,Pn)- 
According to our definition of the momenta, if Pi has positive energy, A4{pi,p2, . . . ,Pn) de- 
scribes a physical process that produces the parton Ai in the final state; if Pi has negative 
energy, Jli{pi,p2, . . . ,Pn) describes a physical process produced by the collision of the an- 
tiparton Ai in the initial state. 

The matrix element Ai{pi,p2, ■ ■ ■) can be evaluated in QCD perturbation theory as a 
power series expansion (i.e., loop expansion) in the QCD coupling gs (or, equivalently, in 
the strong coupling as = g^/{4:n)). We write 

M = >l(o) +7W(2) + ... , (1) 

where Ai^^'' is the tree-leve|l] scattering amplitude, Ai^^'^ is the one- loop scattering ampli- 
tude, Ai^'^^ is the two- loop scattering amplitude, and so forth. Note that in Eq. ([T]) we have 
not written down any power of gs- Thus, At^'^^ includes an integer power of gs as overall 
factor, and Ai^^'' includes an extra factor of g^ (i.e., Al^^y At^^^ oc g^). Throughout the 
first part of the paper (Sects. [2]- l5.2p . we always consider unrenormalized matrix elements, 
and gs denotes the bare (unrenormalized) coupling constant. 

Physical processes take place in four-dimensional space time. In the four- dimensional 
evaluation of the one-loop amplitude A^'-^^ one encounters ultraviolet and IR divergences 
that have to be properly regularized. The most efficient method to simultaneously reg- 
ularize both kind of divergences in gauge theories is dimensional regularization in d ^ 4 
space-time dimensions. We work in c? = 4 — 2e space-time dimensions, and the dimensional- 
regularization scale is denoted by fi. Unless otherwise stated, throughout the paper we 
formally consider expressions for arbitrary values of d = 4 — 2e (equivalently, in terms of 
e-expansions, the expressions are valid to all orders in e). 

We are interested in studying the behaviour of Ai{pi,P2, ■ ■ ■ ,Pn) in the kinematical con- 
figuration where two of the external parton momenta become (almost) collinear. Without 
loss of generality, we assume that these momenta are pi and p2- We parametrize these 
momenta as follows: 

p>; = x.p^ + k'^^-^ , ^ = l,2, (2) 

Xi 2p ■ n 

where the light-like (p^ = 0) vector p^^ denotes the collinear direction, while n'^ is an auxil- 
iary light-like (n^ = 0) vector, which is necessary to specify the transverse components k±i 
{k±i -p = k±i ■ n = 0, with k^- < 0) or, equivalently, to specify how the collinear direction is 

^Precisely speaking, Ai^^'^ is not necessarily a tree amplitude, but rather the lowest-order amplitude for 
that given process. Thus, A^*-^-* is the corresponding one-loop correction. For instance, in the case of the 
process 77 — >■ gg, A^^") involves a quark loop. 



4 



approached. No other constraints are imposed on the longitudinal and transverse variables 
Xi and k±i (in particular, we have xi + X2 7^ 1 and k±i + k±2 7^ 0). Thus, we can consider 
any (asymmetric) collinear limits at once. Note, however, that the collinear limit is invari- 
ant under longitudinal boosts along the direction of the total momentum = Pi + • 
Thus, the relevant (independent) kinematical variables are the following boost-invariant 
quantities: a single transverse-momentum variable k'^ {k^ = Z2k^^ — Zik^2-> < 0) ^ 
single longitudinal-momentum fraction, which can be either zi or Z2 (or the ratio between 
Zi and ^2)- The longitudinal- momentum fractions zi and Z2 are 

Zi + Z2 = l . (3) 



Xi+ X2 



In terms of these boost-invariant variables, the invariant mass squared S12 = (pi + P2Y 
the system of the two 'collinear' partons is written as 

S12 = 2pi ■ P2 = . (4) 

Zl Z2 

We also define the following light-like (P^ = 0) momentum P^: 

P" = (pi+P2)^- ,, ■ (5) 



In the kinematical configuration where the parton momenta pi and p2 become collinear, 
their invariant mass S12 vanishes, and the matrix element Ai{pi,P2, ■ ■ ■ ,Pn) becomes sin- 
gular. To precisely define the collinear limit, we rescale the transverse momenta kj_i in 
Eq. ([2]) by an overall factor A (namely, kj_i — )■ A kj_i with i = 1,2), and then we perform the 
limit A — 7- 0. In this limit, the behaviour of the matrix element Ai{pi,p2, ■ ■ ■ ,Pn) is propor- 
tional to 1/A. We are interested in explicitly evaluating the matrix element contribution 
that controls this singular behaviour order by order in the perturbative expansion. More 
precisely, in c? = 4 — 2e dimensions, the four- dimensional scaling behaviour in the collinear 
limit is modified by powers of (A^)~^. Since we work with fixed e, we treat the powers of 
(A^)~^ as contributions of order unity in the collinear limit. 

In summary, considering the limit S12 — )■ 0, we are interested in the singular behaviour: 

1 r 1 

M{pi,p2,...,Pn) = mod (In'^ S12) l + 0(yii;) , (6) 

^JSl2 L J 



where the logarithmic contributions In'^ S12 (/c = 0, 1, 2, . . . ) ultimately arise from the power 
series expansion in e of terms such as (512)"*^. These logarithmic contributions are taken 
into account in our calculation, while the corrections of relative order 0{^/sl2 ) are sys- 
tematically neglected. 

As is well known [161 [IZ], the singular behaviour of tree-level scattering amplitudes 
in the collinear limit is universal (process independent) and factorized. The factorization 
structure is usually presented at the level of colour subamplitudes [17], in a colour-stripped 
form. In Ref. pBj, we proposed a formulation of collinear factorization that is valid directly 
in colour space. Here, we follow this colour space formulation, which turns out to be 
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particularly suitable to the main purpose of the present paper, namely, the general study 
of the SL collinear limit at one-loop and higher-loop orders. 

To directly work in colour space, we use the notation of Ref. |35] (see also Ref. [1]). The 
scattering amplitude depends on the colour indices {ci,C2, . . . } and on the spin (e.g. 
helicity) indices {si, S2, . . . } of the external QCD partons; we write 

-M'=^'^^'-'^"^'^''^'-''"(Pl,P2,...,Pn) . (7) 

We formally treat the colour and spin structures by introducing an orthonormal basis 
{|ci, C2, . . . , Cn) ® S2, . . . , Sn)} in colour + spin space. The scattering amplitude in 
Eq. (I7j) can be written as 

-M^^''=^'-^^^'^^'-(Pl,P2,...)= ((Cl,C2,...|®(Si,S2,...|) \Mip,,P2,...)) . (8) 

Thus |A^(pi,p25 • • • 5Pn)) is a vector in colour + spin (helicity) space. 

As stated at the beginning of this section, we define the external momenta p^'s as 
outgoing momenta. The colour indices {ci,C2, . . . c„} are consistently treated as outgoing 
colour indices: q is the colour index of the parton Ai with outgoing momentum pi (if pi has 
negative energy, q is the colour index of the physical parton Ai that collides in the initial 
state). An analogous comment applies to spin indices. 

Having introduced our notation, we can write down the colour-space factorization for- 
mula [26] for the collinear limit of the tree- level amplitude A^*^°^ We have 

|A^(°)(pi,P2,...,Pn))^5p(°)(pi,j92;P) \M^^\p,...,p^)) , (9) 

which is valid in any number c? = 4 — 2e of dimensions. The only approximation (which is 
denotecii] by the symbol '~') involved on the right-hand side amounts to neglecting terms 
that are less singular in the collinear limit (i.e. the contributions denoted by the term 
(9(7^) in Eq. (ED). 

The tree-level factorization formula relates the original matrix element (on the left- 
hand side) with n partons to a corresponding matrix element (on the right-hand side) with 
n — 1 partons. The latter is obtained from the former by replacing the two collinear partons 
Ai and A2 (with momentum pi and p2, respectively) with a single parent parton A, whose 
momentum is P (see Eq. ([5])) and whose flavour is determined by flavour conservation of 
the QCD interactions. More precisely, A is a quark (an antiquark) if the two collinear 
partons are a quark (an antiquark) and a gluon, and A is a gluon otherwise. 

The process dependence of Eq. (jO]) is entirely embodied in the matrix elements on 
both sides. The tree-level factor Sp''^\pi,p2] P), which encodes the singular behaviour 
in the collinear limit, is universal (process independent) and it does not depend on the 
non-coUinear partons with momenta ps, . . . ,Pn- It depends on the momenta and quantum 
numbers (flavour, spin, colour) of the partons that are involved in the collinear splitting 
A — )• A1A2. According to our notation, Sp^^\pi,p2', P) is a matrix in colour+spin space, 
named the splitting matrix [26j. 

■'■The symbol '~' is used throughout the paper ■with the same meaning (namely, neglecting terms that 
are non-singular or vanishing in the collinear limit) as in Eq. (|9]). 
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The splitting matrix acts between the colour space of the n—1 partons of |A^(°^(P, ..,Pn)) 
and the colour space of the n partons of the original amplitude \Ai^^\pi,P2, ■■,Pn))- Because 
Sp^'^\pi,P2', P) does not depend on the non-coUinear partons, their colour is left unchanged 
in Eq. ([9]) (precisely speaking, Sp^^^ is proportional to the unit matrix in the colour subspace 
of the non-coUinear partons). The only non-trivial dependence on the colour (and spin) 
indices is due to the partons that undergo the coUinear splitting A — )■ A1A2. Making this 
dependence explicit, we have 

5p(0)(ci,c.;c)(^^^^^.p) ^ ^c^^c^l 5p(0)(p,,p2;P) |C) , (10) 

where ci, C2 and c are the colour indices of the partons Ai, A2 and the parent parton A. The 
colour indices of gluons, quarks and antiquarks are actually different; we use the notation 
c = {a} = 1, . . . , — 1 for gluons and c = {a} = 1, Nc for quarks and antiquarks, 
where N^, is the number of colours. A colour matrix of the fundamental representation 
of the gauge group is denoted by ta^Q^ , and the structure constants are fate "we use the 
following normalization: 

[t'',t']=tfatX , Tr(n^) = ^5'^^ . (11) 

The matrices t"' are hermitian {{f^y = (t")) and the structure constants fabc are real 
ifabc = fabc)- There are four different flavour-conserving configurations A — > AxA^. The 
corresponding explicit form of the tree-level splitting matrix is: 



Spf;^^^-^'^^-'-\p^,P2-,P)=fi'gstZo.—HPi)^iP2)u{P) , (12) 

-512 



q qi92 



^pSf""""^(Pi,P2;P)=/.^^?s (-aj —v{P)^{P2)v{p,) , (13) 

-512 



9 qm 



5pSy"""^^"^(Pi,P2;P) = /i^^?sC«. — ^(Pi)r(^)^(P2) , (14) 

512 

9 9i92 

SpfJf-'--''^){p,,P2- P) = f^^gsz — (15) 

S12 

e{pi) ■ e{p2) pi ■ e*{P) + e{p2) ■ e*{P) P2 ■ e{pi) - e{pi) ■ e*{P) pi ■ 5(^2) 

where u{p) and v{p) are the customary Dirac spinors and e'^{p) is the physical polarization 
vector of the gluon (e* is the complex conjugate of e^). Spin indices play no relevant 
active role in the context of the main discussion of the present paper. They are embodied 
in the parton wave functions u,v,e and are not explicitly denoted throughout the paper. 
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The explicit expressions of Eqs. flT^ - fllSp in a definite helicity basis can be found in the 
hterature (see, for instance, the Appendix A of the second paper in Ref. [12]). 

We briefly comment on the relation between Eq. (Q and the customary collinear- 
factorization formulae for colour subamplitudes (see also the Appendix|X]). The colour-space 
factorization formula (|9]) is valid for a generic matrix element \Ai^^\pi,P2, ■■,Pn))] in partic- 
ular, the factorization formula does not require any specifications about the colour structure 
of the matrix element. CoUinear factorization of QCD scattering amplitudes is usually dis- 
cussed [16] upon colour decomposition of the matrix element. The colour decomposition, 
whose actual form depends on the specific partonic content of the matrix element (e.g., on 
the number of gluons and quark- ant iquark pairs), factorizes the QCD colour from colour- 
less kinematical coefficients, which are called colour subamplitudes (see, e.g., Ref. fTT]). 
Colour subamplitudes fulfil several process-independent properties, including collinear fac- 
torization. In the region where the two parton momenta pi and p2 become collinear, the 
collinear-factorization formula of the colour subamplitudes is a colour-stripped analogue 
of Eq. (|9]): the colour vectors |A^*^°^) on both sides are replaced by corresponding colour 
subamplitudes, and the colour matrix Sp^'^\pi,p2', P) is replaced by a universal kinematical 
function, which is called splitting amplitude and is usually denoted by Split*-°''(pi,p2; P) (see 
also the Appendix |A]) . 

In the case of the tree-level collinear splitting of two partons, the relation between the 
splitting matrix Sp^^^ and the splitting amplitude Split is particularly straightforward. 
Indeed, having fixed the fiavour of the partons in the collinear splitting process A — )■ A1A2, 
the corresponding splitting matrix Sp^^^ involves a single (and unique) colour structure 
and, therefore, we have a direct proportionality relation: 



where the colour matrix on the right-hand side is obtained by simple inspection of Eqs. f|T2|) - 
f|T5|) (see the colour factors f^^, in Eqs. f fT2|) -f lT4|) and the colour factor ifaia2a in Eq. f|T5|) ). 

As discussed at the beginning of this section, the outgoing momenta pi^s of M{pi, pn) , 
depending on the sign of their energy, actually describe different physical processes, which 
take place in different kinematical regions. Correspondingly, the collinear splitting A — )■ 
A1A2 formally describes different physical subprocesses, which take place in either the TL 
(if P1P2 > 0) or SL (if PiP2 < 0) regions. In these regions, the collinear variables zi, Z2 and 
s\2 in Eqs. ([3]) and ([1]) are constrained as follows: 



S'p^°-'(pi,P2; -P) oc (colour matrix) X Split ''°'*(pi,p2; -P) , 



(16) 



TL : 



SL : 



■512 > , Z\Z2 > , 
S12 < , Z1Z2 < . 



(17) 
(18) 



The most relevant physical subprocesses are the customary subprocesses: 



• TL {p1>0 , p°> 0): 



A* ^ Ai{z) A2{1 - z) , with z = zi = 1 - Z2 , 0<;s<l, 



(19) 



• SL (j9? < , > , P' 



<0): 
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Ai ^ A\z) ^2(1 - z) , with z = — = , < z < 1 . (20) 

Zi I- Z2 

In the TL subprocess of Eq. (fT9|) . the partons Ai, A2 and the parent parton A are phys- 
ically produced into the final state; the collinear decay of the parent parton A*, which is 
slightly off-shelil (with positive virtuality) in the vicinity of the collinear limit, transfers 
the longitudinal- momentum fractions z and 1 — z to Ai and y42, respectively. In the SL 
subproces^ of Eq. ( |20l) . the physical parton Ai, which collides in the initial state, radiates 
the physical parton ^2, with longitudinal- momentum fraction 1 — in the final state; the 
remaining fraction, z, of longitudinal momentum is carried by the accompanying ('parent') 
parton A (which is slightly off-shell, with negative virtuality, in the vicinity of the collinear 
limit) that replaces Ai as physically colliding parton in the initial state. 

There are two other physical subprocesses that are kinematically allowed: the TL sub- 
process A1A2 — 7- A (parton-parton fusion into an initial-state parton) is allowed iip\ and 
P2 are both negative, and the SL subprocess AiA* — )■ A2 (parton-parton fusion into a final- 
state parton) is allowed if < , P2 > > ^'^'^ > 0- The subprocess A1A2 A* (the 
initial-state parton A is produced by the fusion of the two initial-state collinear partons 
Ai and A2) occurs if A^(pi,p2, • • • ) corresponds to a physical process with at least three 
colliding particles in the initial state (the partons Ai, A2 and, at least, one additional parti- 
cle). The subprocess AiA* — )■ A2 (the final-state parton A2 is produced by the fusion of the 
initial-state parton Ai and the final-state parton A*) occurs if A4(pi,P25 • • • ) corresponds 
to a physical process in which the initial state contains the parton Ai and, in addition, 
either one massive particle or (at least) two particles. Owing to these kinematical features, 
these subprocesses are less relevant in the context of QCD hard-scattering processes. 

The splitting matrix Sp^^\pi,p2', P) in the factorization formula applies to any 
physical subprocesses, in both the TL and SL regions. Strictly speaking, the explicit 
expressions in Eqs. f|T2|) - f|T5|) refer to the TL region where the energies of pi,p2 and P are 
positive. The corresponding expressions in other kinematical regions are straightforwardly 
obtained by applying crossing symmetry. If the energy of the momentum P (P = P11P2 
or P) is negative, the crossing relation simply amounts to the usual replacement of the 
corresponding wave function (i.e., u{P) ^ v{—P) and e{P) -H- e*{—P)). 



3 One-loop amplitudes: time-like collinear limit 

In this section we consider the collinear behaviour of the one-loop QCD amplitudes Ai^^^ 
in Eq. ([T]). We use the same general notation as in Sect. [H However, we anticipate that 
the results are valid only in the case of the TL collinear splitting (i.e., S12 > 0). 

The singular behaviour of Ai^^\pi,P2, ■ ■ ■ iPn) in the region where the two momenta pi 
and P2 become collinear is also described by a factorization formula. The extension of the 

§Wc introduce the star superscript in A* to explicitly remind the reader that the parent parton A is 
ofF-shell before approaching the collinear limit. 

^The corresponding subprocess with pi ■(-> p2 is trivially related to Eq. ((20|) by the exchange 1 -;-> 2 of 
the parton indices. 



9 



tree-level colour-space formula to one-loop amplitudes is [2S] 



\M^'\p^,P2,...,Pn)) 



+ Sp^'\p,,P2-,P) \M^'\P,...,Pn)) 



(21) 



The 'reduced' matrix elements on the right-hand side are obtained from Ai{pi,p2, . ■ ■ ,Pn) 
by replacing the two collinear partons Ai and A2 (with momentum pi and p2, respectively) 
with their parent parton A, with momentum P. The two contributions on the right- 
hand side are proportional to the reduced matrix element at the tree-level and at the 
one-loop order, respectively. The splitting matrix Sp^^"^ is exactly the tree-level splitting 
matrix that enters Eq. (Q. The one-loop splitting matrix Sp^^\pi,p2', P) encodes new (one- 
loop) information on the collinear splitting process A — )■ Analogously to Sp^'^\ the 
one-loop factor Sp^^\pi,p2] P) is a universal (process-independent) matrix in colour+spin 
space, and it only depends on the momenta and quantum numbers of the partons involved 
in the collinear splitting subprocess. 

Within the colour subamplitude formulation, the collinear limit of two partons at the 
one-loop level was first discussed in Ref. p3l| by introducing one-loop splitting amplitudes 
Split*^^^(pi,P2; -P), which are the one-loop analogues of the tree-level splitting amplitudes 
mentioned in Sect. [2J A proof of collinear factorization of one- loop colour subamplitudes 
was presented in Ref. [21]. Explicit results for the splitting amplitudes Split in ci = 4 — 2e 
dimensions (or, equivalently, the results to all orders in the e expansion) were obtained in 
Refs. [121 [25]. 

The relation between the one-loop factorization formula (l2Ti) and its colour subampli- 
tude version is exactly the same as the relation at the tree level (see also the Appendix 
The main point is that the one-loop splitting matrix Sp^^\pi,p2', P) involves a single colour 
structure (more precisely, there is a single colour structure for each flavour configuration 
of the splitting processes A — )■ A1A2), and this colour structure is the same structure that 
occurs in the tree-level splitting matrix Sp^^\pi,p2', P). In other words, the proportion- 
ality relation in Eq. f[TB]l is valid also at the one-loop level: we can simply perform the 

replacements Sp^^\pi,p2; P) Sp^^\pi,p2]P) and Split^°^(pi,p2; -P) Split^^^(pi,p2; -P)- 
Therefore, from the known Split^^^ [121 [25] we directly obtain the corresponding Sp^^\ 

We now comment on the kinematical structure of Sp^^\pi,p2] P), i.e. on the momen- 
tum dependence of Split''^^(pi,p2; -P)- Apart from the overall proportionality to the wave 
functions u,v,e of the collinear partons (which is analogous to that in Eqs. (fT2|) - (fT5|) ). 
the kinematical structure [12, 25j depends on two different classes of contributions. One 
class contains all the contributions that have a rational dependence on the momenta; the 
other class contains transcendental functions (e.g., logarithms and polylogarithms) and, in 
particular, transcendental functions of the momentum fractions Zi and Z2. 

Considering d = 4 — 2e space-time dimensions, the one-loop integrals introduce the 
dimensional factor /i^*^. Since the one- loop corrections to Sp^^\pi,p2', P) are dimensionless, 
Sp^^\pi,P2', P) necessarily includes the overall factor 




(22) 
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where the iO prescription follows from usual analyticity properties of the scattering am- 
plitudes. Apart from this overall factor, the trascendental dependence of the two-parton 
collinear limit at one-loop order turns out to be entirely captured [121 [25] by a single hyper- 
geometric function, namely, the function 2-^1(1, l — e;x). The integral representation of 
this hypergeometric function is 



2^1(1, -e;l -e;x) 



dt t 



-l-e 



We thus define the following function: 



/(e; l/x) 



2Fi(l,-e;l-e;l-x)-l 



The expansion of this function in powers of e is as follows: 



/(e; l/x) = Inx — e 
where the dilogarithm function Li2 is 



Li2(l -x) + ^e'= Life+2(1 



X] 



k=l 



'Lioix) = - 



dt 



ln(l - 1) 



and the polylogarithms Lifc+i (with /c = 2, 3, . . . ) are defined [36] by 

^ dt 



LWi(a:) = ^ J {\ntf-' ln(l - xt) 



(23) 



(24) 



(25) 



(26) 



(27) 



As recalled in Sect. [21 there are four different flavour configurations in the collinear 
splitting process A — )• A1A2. Considering the corresponding explicit results of Refs. [121 [25], 
the one-loop splitting matrix Sp^^\pi,p2; P) of Eq. (^Tt) can be written in the following 
general (and compact) form: 

Sp^'\p^,P2■,P) = Sp^i\p^,p2-P) + Ic{pi,P2;P) Sp^''\p^,p2,P) . (28) 

The factor Ic{pi,P2', P) is specified below. Having specified this factor, the term Sp^^^ 
on the right-hand side of Eq. f[28|l can be extracted, in explicit form, from the results in 
Refs. fl2\ [25] . We do not report the explicit form of Sp^^\ since it has no relevant role in 
our discussion of the relation between the TL and SL collinear limits. In this respect, the 
only relevant property of -S'p^^ (which follows from our definition of Ic) is that it contains 
only terms with rational dependenc^ on the momenta pi,p2 and P. Moreover, all the 
terms of Sp^^^ that are IR or ultraviolet divergent in c? = 4 dimensions (i.e. all the e poles) 
are collected in the factor Ic and, thus, removed from Sp^^\ Therefore, Sp^^^ is finite if 
we set e = 0. 

The term Ic x Sp^^^ on the right-hand side of Eq. f[28|) contains all the IR and ultraviolet 
divergences of and, more importantly, it collects the entire dependence of the collinear 

"Precisely speaking, this statement is true modulo the dependence on the overall factor in Eq. ([22|). 
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behaviour at one-loop order on transcendental functions (modulo the function in Eq. 
which also appears in Sp^^^). The explicit expression of the factor Ic for a generic splitting 
process A — )■ A1A2 is 

Ic{Pi,P2]P) = glcv (~^^^2~") 

X I ^ (Ci2 - - C2) + ^ (712 - 71 - 72 + &o) (29) 
(Ci2 + - C2) f{t- z^) + (Ci2 + C2 - Ci) /(e; Z2) ] I , 



where cp is the typical volume factor of d-dimensional one-loop integrals: 

^ r(i + ^)r-(i-.) , 

(47r)'"'r(l -2e) 

The coefficients Ci,C2 and C12 are the Casimir coefficients of the partons Ai,A2 and A; 
explicitly, C = Cp = {N^ — l) / (2iVc) if the parton is a quark or antiquark, and C = Ca = Nc 
if the parton is a gluon. Analogously, the coefficients 71,72 and 712 refer to the flavour of 
the partons Ai, A2 and A; explicitly, we have 

l, = lg = lcF , ^g = ^illCA-2Nf) , (31) 

where Nf is the number of flavours of massless quarks. The coefficient 60 is the first 
perturbative coefficient of the QCD /3 function, 

bo = ^{nCA-2Nf) . (32) 

Note that, in to our normalization, we have = 'jg. 

The coefficients of the e poles in Eq. fl29]) agree with those of the general structure 
presented in Eq. (11) of Ref. [26]. The single-pole term proportional to 60 is of ultraviolet 
origin; it can be removed by renormalizing the splitting matrix Sp{pi,p2] P) (we recall that 
we are considering unrenormalized matrix elements and, correspondingly, unrenormalized 
splitting matrices). The other pole terms are of IR origin. The double-pole terms (which 
are proportional to the Casimir factors Cp and Ca) originate from one-loop contributions 
where the loop momentum is nearly on-shell, soft and parallel to the momentum of one of 
the three partons involved in the collinear splitting subprocess. The single-pole terms with 
7 coefficients are produced by contributions where the loop momentum is not soft, though it 
is nearly on-shell and parallel to the momentum of one of the collinear partons. According 
to Eq. fl2^ . the e expansion of the transcendental function gives f{e;z) = — In^ + C(e); 
therefore, /(e; zi) and /(e; Z2) contribute to Eq. ( 12^ with single-pole terms. The coefficients 
of these single-pole terms are controlled by the Casimir factors Cp and Ca and, hence, 
they originate from one-loop configurations with soft momentum; they are produced by 
contributions where the loop momentum is nearly on-shell, soft and at large angle with 
respect to the directions of the collinear partons. The specific combination of Casimir 
factors in these single-pole terms (namely, C12 + Ci — C2 and C12 + C2 — Ci) originates 
from a colour coherence effect (see Eq. (jH]) and the comments below it). 
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As anticipated at the beginning of this section, the one-loop factorization formula fl2ip 
and the explicit results in Eqs. (128|) and (!29|) (or, equivalently, the one-loop splitting am- 
plitudes in Refs. [23]-[25]) are valid in the case of the TL collinear limit (see Eq. ( [T7|) ). At 
the tree level, the TL and SL collinear limits are related by exploiting crossing symmetry, 
and the corresponding splitting matrix Sp^^\pi,p2; P) is simply obtained by applying the 
(wave function) crossing relations mentioned at the end of Sect. [2J At the one- loop level, 
we have to deal with the splitting matrix Sp^^\pi,p2', P) in Eq. ( 128|) . and we can try to 
proceed in an analogous way. Using crossing symmetry, the treatment of the one-loop 
contribution Sp^^\pi,p2; P) is straightforward; Sp^^\pi,p2; P) contains (i) wave function 
factors, which are treated by the corresponding crossing relations, and [ii) rational func- 
tions of the collinear momenta, which are invariant under crossing. The one-loop troubles 
originate from the factor Ic{pi,P2', P), since it contains /(e; Zi) and /(e; ^2)- 

The function /(e; x) (see Eq. (!24|) ) has a branch-cut singularity if the variable x is real 
and negative. The branch-cut singularity arises from the corresponding singularity of the 
hypergeometric function 2-F'i(l, 1 — e; 1 — In the case of the TL collinear limit, 

zi and Z2 are both positive (see Eq. (IT7I) and recall that Zi + Z2 = 1), and the functions 
/(e; zi) and /(e; Z2) are both well-defined. In the case of the SL collinear limit, one of the 
two variables zi and Z2 necessarily has a negative value (see Eq. f|T7j) ): therefore, one of 
the two functions, either f{e;zi) or f{e;z2), in Eq. fl29|) is necessarily evaluated along its 
branch-cut singularity and, hence, it is ill-defined. 

In summary, the issue of the TL vs. SL collinear limits is as follows. The results in 
Eqs. (1211) . fl2Sl) and (12^ cannot be extended from the TL to the SL collinear limit by using 
crossing symmetry, since this leads to ill-defined results (mathematical expressions). As 
shown in the next section, the solution of the issue involves not only the (mathematical) 
definition of the function /(e; x) along (or, more precisely, in the vicinity of) its branch-cut 
singularity, but also the introduction of new physical effects. 



4 One-loop amplitudes: general (including space-like) 
collinear limit 

4.1 Generalized factorization and violation of strict collinear fac- 
torization 

The extension of the colour-space collinear formula in Eq. f l2T]) to general kinematical 
configuration which include the two-parton collinear limit in the SL region, is 

\M^'\p,,P2,...,Pn)) ^ 5p«(Pl,P2;P;P3,...,Pn) \M^'\P,...,Pn)) 

+ Sp^'\puP2;P) \M^'\P,...,Pn)) ■ (33) 

The essential difference with respect to Eq. fl^T]) is that the one-loop splitting matrix Sp^^^ 
on the right-hand side of Eq. (133|1 depends not only on the collinear partons but also on the 

**In Appendix El we illustrate the SL collinear limit of colour subamplitudes for the specific case of pure 
multigluon matrix elements at the one-loop level. 
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momenta and quantum numbers of the non-colhnear partons in the original matrix element 
\A4^^\pi,P2, ■ ■ ■ ,Pn)) ■ Thus, Sp^^'' is no longer (strictly) universal, since it retains some 
dependence on the process (matrix element) from which the splitting matrix derives. The 
reduced tree-level, |A^'^°^(P, . . . ,Pn)), and one-loop, |A^*^^^(P, . . . ,p„)), matrix elements on 
the right-hand side of Eq. fl551) are the same as those in Eq. fl^T]) : they are still related to 
the original matrix element {pi,P2, ■ ■ ■ ,Pn)) through the same factorization procedure 
that is used in Eq. (12T|) (i.e. in the case of the TL collinear limit). 

The explicit form of the general one-loop splitting matrix Sp^^^ in Eq. fl33l) is 

Sp^^\pi,P2;P;P3, - ■ ■ ,Pn) = Sp%\pi,P2]P) +Ic{Pl,P2;P3,---,Pn) Sp'^^\pi, P2] P) , 

(34) 

where Sp^^ {pi,p2] P) is exactly the same (universal) term as in Eq. fl28|) . The difference 
with respect to the TL expression in Eq. ( l28|) arises from the replacement of Ic{pi,P2', P) 
with Ic{pi,P2',P3,---,Pn)- The term Ic{pi,P2'i P) is a c-number (i.e., colourless) factor, 
while Ic{pi,P2',P3,---,Pn) is a colour matrix. Moreover, Ic{pi,P2', P) depends on the 
collinear variables zi,Z2,Si2 and the flavour of the collinear partons Ai,A2 and A (see 
Eq. (129|) ). while Ic{pi,P2',P3,---,Pn) o,lso depends on the momentum and colour of the 
non-collinear partons. 

The expression of the colour operator Ic can be presented by using the same notation 
as in Eq. ( l29l) . We can also exploit the fact that in any kinematical configurations (see 
Eqs. ( ITTl) and ( |T8|) ) one of the two collinear variables, zi and Z2, necessarily has positive 
values (recall that zi + Z2 = 1). Therefore, with no loss of generality, we can set (choose) 

^1 > , 

and write the following explicit expression of the colour operator: 

IciPl,P2]P3,---,Pn) = glcT ( ] (35) 



X 



I ^ {Ci2 - Ci - C2) + ^ {112 - 7i - 72 + 

n 

[Ci2 + Ci - /(e; ^i) - 2 ^ T2 • T, /(e; Z2 - ^Os.s) 



Here, the subscript j (j = 3, . . . , n) refers to the non-collinear parton with momentum pj, 
and Sj2 = {pj +^2)^ is the invariant mass squared of the system formed by the j-th non- 
collinear parton and the collinear parton A2. The colour charge (matrix )[1| of the parton 
with momentum pk {k = 1,2,3, ... , n) is denoted by T^, and we define -Ti = J^c'^k'^t 
(c=l,...,iV2-l). 

On the right-hand side of Eq. (!35|) . the function f{e;zi) is well-defined, since zi > 0. 
The functional dependence on Z2 is also well-defined, since it is given by either /(e; Z2 — iO) 



^^We use the same notation as in Refs. [35l [T]: more details about colour charges and colour algebra 
relations can be found therein. 
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if Sj2 > 0, or /(e; Z2 + iO) if Sj2 < 0. Owing to the iO prescription, if Z2 < 0, the function 
/(e; Z2 ± iO) is always evaluated either above or below its branch-cut singularity. The 
presence and structure of the branch-cut singularity are physical consequences of causality, 
as discussed in Sect. 14.51 

Comparing Eq. fl29l) with Eq. fl35|) . we see that the difference is due to a single con- 
tribution. In the TL case of Eq. fl2I?l) this contribution is proportional to the following 
term: 

5(pi,P2; P) = -l (Ci2 + C2- Ci) fie; Z2) , (36) 
while in the general case of Eq. (1351) this term is replaced by the following colour operator: 

2 " 

<5(pi,P2;P3,---,Pn) = + - -Tj /(e;2;2 -iOsj-a) ■ (37) 

^ i=3 

The 'analytic' continuation from the TL coUinear region to a generic coUinear region is thus 
achieved by the introduction of a colour-energy correlated iO prescription. The main new 
physical effect in Eq. fl37|) is the presence of colour correlations between the coUinear and 
non-collinear partons. This effect produces violation of strict (naive) coUinear factorization 
of the scattering amplitudes. 

A detailed derivation of the results in Eqs. (I33l) -( |35l) . including the extension to the 
multiple coUinear limit of m (m > 3) parton momenta (see Sect. 15.11) . will be presented in 
a forthcoming paper. In Sect. 15. 3[ we illustrate the explicit computation of the IR divergent 
part of the one-loop splitting matrix Sp*-^^ in Eq. ( l33l) . The result of this computation can 
be regarded consistency check of Eqs. (I33l)-( l35i) . 

In the following we discuss some consequences of the results in Eqs. (1331) -( l35l) . To this 
purpose, we first present some colour algebra relations. An important relation is colour 
conservation; we have 

n 

J] Tfc = , (38) 
fc=i 

or, equivalently, 

n 

Y,T, = -{T,+T2) . (39) 

i=3 

Precisely speaking, the relations in Eqs. f l38|) and fl39l) are valid in operator form when the 
colour charges act onto an overall colour-singlet vector, with n partons, in colour space. 
Such vectors are, for instance, the matrix element vector |A^(')(pi,p25 • • • ^Pn)) i} = 0,1) 
and the vector Sp^^\pi,p2; P) \Ai^^\P, . . . ,Pn)) in Eq- (155]) . As a consequence of colour 
conservation in the tree- level coUinear splitting process A ^ A1A2 , we hav6§ 

(T1 + T2) 5p(°)(pi,p2;P) = 5p(°)(pi,p2;P) Tp , (40) 

where Tp denotes the colour charge of the parent coUinear parton A. We also recall that 
T1 = Cfc, where Ck is the Casimir factor of the k-th parton. Therefore, Eq. fHOj) implies: 

(Ti + T2)2 Sp^'\puP2; P) = C12 5p(°)(pi,p2; P) , (41) 

^The relation in Eq. (|40| is also valid when replacing S'p^°^ with the one- loop contribution 
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or, equivalently, 



2 T2 ■ (Ti + T2) 5p(°)(pi,p2; P) = (Ci2 + C2 - Ci) 5p(°)(Pi,P2; P) , (42) 

which follows from the identity (Ti + - TI - TI = 2 Ti ■ T2 

Using simple colour algebra relations, we can easily show that the general results in 
Eqs. lead to the TL results illustrated in Sect. El Considering the collinear 

limit in the TL region, we have Z2 > (see Eq. f|T7j) ) and, therefore, the iO prescription in 
/(e; Z2 — i0sj2) is harmless. Removing the iO prescription on the right-hand side of Eq. (137]) . 
we have 

2 " 

S{puP2;P3, ...,Pn) = + - /(e; Z2) T2 ■ J] , Z2>0 , (43) 

^ i=3 

and we can perform the sum over the colour charges of the non-coUinear partons. Since 
6 acts onto the colour vector Sp^^\pi,p2; P) \M^^\P, . . . ,p„)) (see Eqs. (133])-([35D), the 
sum over the colour charges can be carried out explicitly by using Eqs. (139|) and (H2l) and, 
hence, Eq. (^3]) becomes 

S{pi,P2;P3,---,Pn) = -^{Cu + C2~Ci^ f{e;z2) , Z2>0 , (44) 
which is equal to the TL expression in Eq. (1361) . 



In summary, the absence of evident colour correlations in Eq. f l44j) or, equivalently, the 
validity of strict collinear factorization in the TL collinear limit is a physical consequence 
of colour coherence (and colour conservation). In the case of the TL collinear limit, the 
non-coUinear partons act coherently as a single parton, whose colour charge is equal to the 
total charge of the non-collinear partons (see Eq. (H3l) ). Owing to colour conservation, this 
colour charge is equal (modulo the overall sign) to the colour charge of the parent parton; 
therefore, the total contribution of the interactions (which are separately not factorized) 
of a collinear parton with the non-collinear partons is effectively equivalent to a single 
interaction with the parent parton. This interaction factorizes and produces the colour 
coefficient in Eq. 



Incidentally, exploiting Eq. (l44l) . we note that we can remove the constraint 2:1 > and 
write the colour operator of Eq. (l35l) in a form that has a manifestly symmetric dependence 
on the variables of the two collinear partons Ai and A2. This symmetric form is 

Ic{Pl,P2]P3,---,Pn) = glcT (^^) 

r 1 

X 



1 ^ {Ci2 - Ci - C-^ + ^ ^712 - 7i - 72 + 

2 " 

+ -^T, ■ T^f{e-Zi-iQsji) + T2f{e]Z2-iQsj2) 



and it is obtained from Eq. f l35|) simply by using Eqs. ( 137]) and dS]), with the replacement 
of the subscripts 2 o 1 (i.e., ^2 ^ ^1, T2 ^ Ti and so forth). 
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4.2 The collinear limit of amplitudes with n = 3 QCD partons 



The simplest case in which the two-parton colhnear hmit can be studied occurs when the 
original scattering amplitude Ai{pi,p2, ■ ■ ■ ,Pn) involves only n = 3 QCD partons and, 
necessarily (because of kinematics), additional colourless external legs (with non- vanishing 
momentum). In this case the colour algebra of the operator Ic (or, simply, d) can be 
carried out in closed form. Setting n = 3 in Eq. fl57|) and using the colour-charge relations 
([39D and (02]), we obtain 

S{puP2;P3) = +- T2-T3 f{e-Z2-i0s23) = -- T2-{Ti+T2) f{e;Z2-lOS23) 
e e 

= - ^ (C12 + C2-C1) /(e;^2-^0s23) . (46) 

We see that the operator S is proportional to the unit matrix in colour space. Nevertheless, 
the c-number function in Eq. ( I46l) still retains process-dependent features that derive from 
the violation of strict collinear factorization (i.e., from the iO prescription in /(e; ^2~^0s23)). 

To remark these process-dependent features, we consider the specific exampl^ in which 
the external legs of the matrix element M{pi,P2,P3) = -^iQiPi): g{P2),Q{P3)',l*) are a 
gluon, a quark-antiquark pair and an off-shell photon 7*. The off-shell photon can be 
coupled to a lepton pair, thus leading to the partonic subprocess of different physical 
processes, such as, hadron production in e~^e~ annihilation (e+e~) or in lepton-hadron 
deep- inelastic scattering (DIS), and the production of lepton pairs through the Drell-Yan 
(DY) mechanism in hadron-hadron collisions. These different processes simply require the 
analytic continuation of the same matrix element, A^(5'(pi), ^'(^2), ^(Ps); 7*), to different 
kinematical regions, which are specified by the sign of the 'energies' p^ of the outgoing 
momenta Pi. To be definite, we fix > and we examine the physical partonic processes 
(see Fig. [T]) 

7* ^ q{Pi) + 9{P2) + qiPs) , p? > 0, p° > (e+e-) , (47) 
7* + q{-Pi) ^ 9iP2) + q{p3) , pI<0<pI (DIS) , (48) 
qi-P3) + qi-Pi) ^ 7* + 9{P2) , < 0, < (DY) , (49) 

in the limit where the momenta pi and p2 become collinear. Using the notation of Eqs. (fT9l) 
and (120|) . the colhnear splitting q — )■ q{pi) g{p2) formally describes two different physical 
subprocesses: in e~^e~ annihilation (Eq. fl47p ). we are dealing with the TL subprocess 
q* — i- q{z) g{l — z), where the final-state gluon is collinear to the final-state quark; in 
DIS (Eq. fj48|) ) and DY (Eq. fj49l) ) production, we are deahng with the SL subprocess 
q ~^ q*{z) ~ where the final-state gluon is collinear to the initial-state antiquark. At 
the tree-level, these two physical sphtting subprocesses are related by crossing symmetry, 
namely, by the exchange of the final-state quark with the initial-state antiquark. 

The generalized factorization formula that describes the collinear limit of the one-loop 
scattering amplitude of the processes in Eqs. (H71) - (H9ll involves the 'operator' S of Eq. (H6ll . 
Using Eq. fH6|) in the specific kinematical regions of Eqs. f H7|) -f H9|) . we have 

completely analogous example is obtained by replacing the off-shell photon with a Higgs boson (in 
this case, the qq pair can also be replaced by two gluons). 
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q* q{z) g{l- z) (0 < z < 1 , Z2 = I - z > 0): 



d^'^'-\p,,P2;P3) = --CAf{e-Z2) 



(50) 



• g -> q*{z)g{l - z) 



(0 < z < 1 , ^2 = 1 - 1/^ < 0): 



CA/(e;z2-^0) , 
CA/(e;^2 + ^0) . 



(51) 



^^^"^(pi,P2;P3) 



(52) 



The differences between the expressions in Eqs. ( !50l) . ( ISTI) and (!52|) highhght the effect of 
the violation of strict colhnear factorization at the one-loop level. In going from the TL 
expression (l50l) to the related SL expressions (Eqs. ( 15T|) and (l52l)). it is not sufficient to 
use the replacement z ^ Xjz: in fact, the crossing relation z -h^ Xjz has to be supple- 
mented with an appropriate iO prescription. More importantly, the two SL expressions in 
Eqs. (!5T!) and (!52|) are different: although they are dealing with the 'same' SL splitting 
subprocess (radiation of a final-state gluon collinearly to an initial-state antiquark), the 
singular colhnear factors of the scattering amplitude are different since they refer to an 
initial-state antiquark in two different physical processes (the DIS and DY processes). 
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DIS 




Figure 1: Two-parton collinear factorization of the one-loop amplitude 
-^(Q'(Pi)) S'fe)? ^(Ps); 7*) with n = 3 QCD partons in different kinematical configu- 
rations (e~^e~, DIS, DY) related by analytic continuation. The tree-level (pointlike 
vertices) collinear splitting subprocesses q ^ qg in the TL (e^e~) and SL (DIS or 
DY) regions are related by crossing symmetry. In the SL region, the one-loop (annular 
vertices) splitting subprocesses depend on the process and are different in the DIS and DY 
configurations. 



4.3 The collinear limit of multiparton amplitudes 

If the original matrix element Ai{pi,P2, ■ ■ ■ ,Pn) has n external QCD partons with n > 4, 
the colour algebra involved in the collinear limit cannot be explicitly carried out in general 
form. The action of the colour operator Ic (or 6) onto the colour vector Sp^^^ x 
indeed depends on the colour-flow structure of the colour vector: this structure has to be 
specified in order to explicitly perform the colour algebra. Some general features of the 
collinear limit, which are independent of the specific colour structure of J^{pi,P2, ■ ■ ■ ,Pn), 
are illustrated below. 
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The function /(e; z) in Eq. fl24p is an analytic function of the complex variable z. When 
the argument z approaches the real axis, /(e; z) has a real and an imaginary part. We write 



/(e;x - iOs) = //j(e;x) - i sign(s) //(e;x) , (x real) 
where the real {fn) and imaginary (//) parts are defined as 



/^(e; x) = Re /(e; x±iO) =- /(e; x + tO) + /(e; x - tO) 



fi{t;x) = Im f{e;x + iO) 



1 



f{e;x + tO)- fie;x-tO) 



(53) 



(54) 



(55) 



We thus consider Eq. ( P7|) and we apply the decomposition in Eq. (135]) to /(e; 2:2 — iOs^^)- 
We obtain 



S{pi,p2;P3, . . . ,p„) = + - T2 ■ 2^ /^(e; ^2) - i sign(sj2) //(e; 2:2) 



(56) 



Note that the term proportional to /i?(e; Z2) does not depend on Sj2, and it can be treated 
as the right-hand side of Eq. fH3|) : we can perform the sum over the colour charges of the 
non-coUinear partons and, using Eqs. fl39l) and fj42l) . we obtain 



<5(Pl,P2;P3, . . . ,P„,) = 5r(Pi,P2; ^') + ^ Si{pi,P2;P3, ■■■,Pv 



where 



5k(Pi,P2; ^) = - 7 (Ci2 + C^2 - /^(e; ^2 

2 
e 



(57) 



(5^ 



Si{PuP2;P3, ...,Pn) = -- T2- \ ^Tj sign(sj2) ) //(e; Z2) . (59) 

0=3 



The expressions in Eqs. fl371) and fl57j) are equivalent. The latter explicitly shows that 
the 'real' (more precisely, hermitian) part, Sr, of the colour operator d is proportional 
to the unit matrix in colour space. Incidentally, the form of 6r is analogous to that of 
the corresponding term, 6, in the TL case (see Eq. fl5B]) ): the only difference is that Z2 
can have negative values in the general case of Eq. (158|) . The 'imaginary' (more precisely, 
antihermitian) part, i6i, of the colour operator 5 has instead a non-trivial dependence on 
the colour charges of the non-coUinear partons; this part is responsible for violation of strict 
collinear factorization. 



For practical computational applications of the factorization formula fl5^ . it is useful 
to expand the one-loop splitting matrix Sp^^^ in powers of e. This eventually requires the 
corresponding expansion of the function /(e; Z2—i0sj2) in Eqs. (l35l) and (!37|) or, equivalently, 
the real functions /^(e; Z2) and //(e; Z2) in Eqs. (!58|) and (!59|) . If Z2 is positive (TL collinear 
limit), the iO prescription is not needed, and Eq. (125!) explicitly gives the expansion of 
f{e;z2). In the case of the SL collinear limit, Z2 is negative and the polylogarithms in 
Eq. (1251) are evaluated close to (either above or below) their branch-cut singularity. To 
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simplify the e expansion in the SL case, we can use the relation between the hypergeometric 
functions 2F1 of complex argument z and 1/z; Eq. (!24|) can thus be written as 



f{e;x-iOs) = /(-e; + - 



r(l + e)r(l-e) 



X — iOs 
1 — X 



(x real, and x < 1) . 



(60) 



If X is negative, /(— e; 1 — x) can safely be expanded as in Eq. ( I25ll . Therefore, the right- 
hand side of Eq. fl60l) gives a simple e expansion of /(e; ^2 — ^05^2) in the SL case. The 
formula f l60|) can also be used to evaluate the real and imaginary parts, fji and //, and 
their expansion in powers of e; we have 



fnie- x) = e{x) /(e; x) + Q{-x) { /(-e; 1 - x) 



fi{e;x) = - Q{-x) it 



Tie 



tan(7re) V 1 — a; 



-X 



-X 



X 



, (61) 



(62) 



We explicitly report the expansion of the functions /ij(e; x) and //(e; x) up to order e^: 

^2" 



e(-x)/H(e;x) = e(-x) <; -ln(-x) + e 

1 



Lig 



-X 



Liq 



fl{e]x) = - 0(-x) TT 



-X 



1 — X I 6 



- In^ 



l-x ) 2 

2 



+ i In^ ^ 



1 — X 



TT 



1 — X/ 3 \1 — X 



O(e^) ,(63) 



—X 



1 _ e In ( — 1 + - In' 

1 — X 



—X 
1 — X 



0(e=^ 



If X > 0, the e expansion of /ij(e; x) is given in Eq. (!25|) . 



(64) 



4.4 The SL collinear limit in lepton— hadron and hadron— hadron 
collisions 



We present some additional general comment^ on the SL collinear limit in the kinematical 
configuration!^ of Eq. f l20|) . This kinematical configuration occurs in the case of hard- 
scattering observables in lepton-hadron and hadron-hadron collisions. 

In lepton-hadron DIS, the partonic matrix elements A4{pi,p2,P3, ■ ■ ■ ,Pn) involve the 
collision of a lepton and a parton in the initial state. The initial-state physical parton 
with 'outgoing' momentum pi and energy p^KO can radiate the collinear parton with 
momentum p2 in the final state {p2 > 0); then, the accompanying ('parent') parton with 
'outgoing' momentum P and energy < acts as initial-state physical parton in the hard- 
scattering subprocess controlled by the 'reduced' matrix elements A^(P,p3, . . . ,Pn)- All the 



§A related discussion, limited to the specific case of amplitudes with n ~ 3 QCD partons, has been 
presented in Sect. 14.21 

^Related comments apply to the SL configuration in which < < Pi ^nd P*^ > (see the final part 
of Sect. [2]). 
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other partons (i.e. the non-colhnear partons) are produced in the final state and, thus, Sj2 > 
(with j > 3). In this kinematical configuration the colour operator 6{pi,p2;p3, ■ ■ ■ ,Pn) is 

(5(pi,p2;P3, . . . ,p„) = - ^ (Ci2 + C2- Ci) fie; Z2 - lO) , s,2 > (j > 3) . (65) 

This expression is obtained from Eq. (1371) by following the same steps as in Eqs. fH5]) 
and (jH]). Actually, the DIS expression in Eq. (165|) differs from the TL expression in 
Eq. (144|) only because of the replacement /(e; Z2) — )■ /(e; Z2 — iO), which, roughly speaking, 
simply produces an additional imaginary part. In Eq. (!65l) we note the absence of explicit 
dependence on the colours and momenta of the non-collinear partons. This implies that 
the two-parton SL collinear limit in DIS configurations effectively takes (mimics) a strictly- 
factorized form. This 'effective' strict collinear factorization is eventually the result of a 
colour-coherence phenomenon, analogously to what happens for the TL collinear limit. 
The interactions (which are separately not factorized) of the collinear parton p2 with the 
non-collinear partons produce imaginary parts; however, since all these interactions involve 
final-state partons, the imaginary parts combine coherently to mimic a single effective 
interaction with the parent parton. This global final-state effect is embodied in the definite 
iO prescription of Eq. ( 165!) . 




Figure 2: Two-parton SL collinear limit [pi ~ ZiP, i = 1,2) of the scattering amplitude 
A4{pi,P2,P3, . . . ,Pn) with n > 4 QCD partons in parton-parton hard-scattering configu- 
rations. Typical (irreducible) colour structure of factorization breaking correlations at the 
one-loop level. 

In hadron-hadron collisions (see Fig. [2]), the SL collinear splitting of Eq. (!20l) takes 
place in a partonic environment that differs from that of lepton-hadron collisions. The 
main difference is due to the fact that the matrix element J^{pi,p2,P3, ■ ■ ■ ,Pn) involves 
the initial-state collision of two QCD partons, the collinear parton with momentum pi and 
another parton, and both partons carry colour charge. The real part 6r of the colour 
operator 6 is the same as in Eq. fl65|) (see also Eq. fl58|) ). while the imaginary part Sj 
is sensitive to the colour charge of the initial-state non-collinear parton. We assign the 
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'outgoing' momentum and the energy pg < to this parton, and we can write the 
following explicit form of 6i: 

Siipi,p2;P3, -^Pn) = + 7 { (^12 + - Ci) + 4 T2 ■ T3} //(e; Z2) , S23 < < s,2 (j > 4) . 

^ (66) 
Combining this imaginary part with the real part (see Eq. fl58|) ). the colour operator <5 can 
be written in the following form: 

S{pi,p2;p3, . . . ,p„) = - ^ (Ci2 + C2- Ci) /(e; z2-i0)+ Ut2- T3 fj{e; Z2) , (67) 

which clearly exhibits the difference with respect to the DIS expression in Eq. f l65|l . 

The expression in Eq. (!66|) is obtained from Eq. (!59|) by using simple colour algebra 
relations. Since S23 < and Sj2 > (j > 4), the colour charge factor in Eq. (!59|) gives 

- ^2-|xj^,sign(.,2)| = T2-|t3-^T,|=T2-(Ti+T2 + 2T3) 



^(C^2 + C2-C^^+2T2■T, , (6^ 



where we have used Eq. (JSHD (namely, J2]=4 Tj = -(T1 + T2 + T3)) and Eq. Inserting 
Eq. (EHD in Eq. ([59]), we get the result in Eq. (166|). 



The initial-state (hence, SL) coUinear splitting in physical parton-parton scattering am- 
plitudes Ai{pi,p2, ■ ■ ■ ,Pn) with n > 4 QCD partons necessarily involves colour correlations 
between the coUinear and non-coUinear partons. A representation of these correlations in 
minimal form is presented in the expression f l66|) of the imaginary part of the colour charge 
operator d (see also Fig. |2]). 



4.5 Factorization, colour coherence and causality 

The factorized structures of Eqs. (!2T|) and (133!) and, inparticular, the violation of strict 
factorization in the case of SL collinear configurationqlH deserve some physical (though 
qualitative) interpretation. We offer some comments about that. 

The universal factorization of scattering amplitudes in the collinear limit is expected on 
the basis of a simple physical picture. When the two partonfl Ai and A2 with momenta 
Pi and P2 become collinear, their invariant mass S12 is much smaller than Sj2 ~ sji, where 
the subscript j generically denotes a non-collinear parton (j = 3,4, ...,n). Therefore, 
we are dealing with a two-scale process. Interactions between the collinear partons take 
place at the small energy scale S12 and, hence, at large space-time distances; whereas, 
interactions between the non-collinear partons and interactions between the collinear and 

"a related effect, in tlie context of ttie computation of logarithmically-enhanced contributions to gaps- 
between-jets cross sections, was described as a 'breakdown of naive coherence' [33l [34] , 

**The qualitative discussion of this subsection can straightforwardly be extended to the multiple collinear 
limit of three or more parton momenta. 
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non-collinear partons require a large energy scale and, hence, they take place at small 
space-time distances. This space-time separation of large and small distances produces 
factorization. The (large-distance) physical subprocess A — )■ A1A2, which generates the 
collinear partons, is factorized from the (short- distance) scattering amplitude that involves 
the parent collinear parton A and the non-collinear partons. 

This simple space-time factorized picture is expressed, at the tree level, by the universal 
factorization formula in Eq. At the one- loop (or higher- loop) level, the same space-time 
picture is too naive, at least in the case of gauge theories. Owing to the long-range nature of 
gauge interactions, partonic fluctuations with arbitrarily-large wavelength (i.e. arbitrarily- 
small momentum) can propagate over widely-separated space-time regions, thus spoiling 
the factorization between the small-distance and large- distance subprocesses involved in 
the collinear limit. 

To be more precise, we refer to one-loop interactions due to a gluon with soft momen- 
tum q (g — )■ 0). The soft gluon produces pairwise interactions between a collinear and a 
non-collinear parton. In the kinematical regions where the soft-gluon momentum has a very 
small angle with respect to the momentum of one of the external partons (either a collinear 
or a non-collinear parton), the one- loop interaction produces factorized (though IR diver- 
gent) contribution^ Therefore, we are left to consider interactions due to a soft wide-angle 
gluon. Each of these pairwise interactions between a non-collinear (j = 3, 4, . . . , n) and a 
collinear [i = 1, 2) parton is proportional to the colour-charge term Tj ■ Ti , and it is not 
factorized (it depends on the colour charge of both the non-collinear and collinear partons). 
These interactions separately spoil collinear factorization. Nonetheless, factorization can 
be recovered because of colour coherence. 

The colour coherence mechanism that leads to factorization is as follows. We consider 
the kinematical region of wide-angle interaction. It is the region where 9jq {6jq is the angle 
between pj and the soft momentum o) is large and, specifically, the region where 6jq can 
reach values that are parametricalljo of the order of 9ji {9ji is the angle between pj and 
the collinear momentum pi). After integration over the soft-gluon momentum, the wide- 
angle interaction proportional to Tj ■ Ti depends on 9ji and, analogously, the wide-angle 
interaction with the other collinear parton is proportional to Tj ■ T2 and depends on 9j2. 
In the collinear limit we have 9ji = 9j2 = 9.p {P is the momentum of the parent parton) 
and, therefore, these two interaction contributions are combined in a single contribution 
that depends on 9jp and is proportional to the colour-charge term Tj ■ (Ti + T2) = Tj ■ Tp . 
This single contribution is exactly the contribution of the wide-angle interaction between 
the non-collinear parton pj and the parent parton P: this contribution is factorized in the 
one-loop scattering amplitude \Ai^^\P,...,Pn)) on the right-hand side of Eqs. (!2T!) and 

dMD. 

These contributions are eventually factorized since the small-angle restriction, namely, the restriction 
to be collinear to the external leg, effectively bounds the soft-gluon interaction to act in the space-time 
region of the physical subprocess that involves the corresponding external leg. 

ttTo be slightly more precise, this wide-angle region is specified by the boundary 9jq < 6ji — 612- Indeed, 
from the viewpoint of the ensemble of the two collinear partons, there is a different wide-angle region. This 
is the region where 92q (or Oiq) is parametrically of the order of 612. In this region, considering the collinear 
limit 9i2 — )■ 0, the system of the non-collinear partons coherently interacts with the collinear parton p2 (or 
Pi), as discussed below Eq. p4l) . 
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In summary, in the coUinear limit the system of the two colhnear partons acts coher- 
ently with respect to non-factorizable wide-angle interactions with each of the non-coUinear 
partons. Owing to this coherent action, these interactions are removed from the collinear 
splitting matrix Sp^^^ and re-factorized in the matrix element \Ai^^\P , . . . , Pn)) of the 
factorization formulae fl2ip and (1331) . 

The colour coherence mechanism that we have just illustrated is completely analogous 
to the mechanism that acts by considering the radiation of two collinear partons and a soft 
gluon in tree- level scattering amplitudes. In this mixed soft-collinear limit, the collinear 
splitting process is factorized from the soft-gluon current of the parent parton (see Sects. 3.4 
and 3.5 in Ref. [11]). However, there is an essential difference between the radiation of a real 
gluon in tree- level amplitudes and the interactions of a virtual gluon in one-loop amplitudes. 
This difference is eventually responsible for the violation of strict factorization in the SL 
collinear limit of one-loop amplitudes. 

A real gluon with soft momentum q is always on-shell (g^ = 0). A virtual gluon instead 
produces both a radiative (roughly speaking, real) and an absorptive (roughly speaking, 
imaginary) contribution to the one-loop amplitude. The radiative contribution is produced 
when the virtual soft gluon is on-shell (g^ = 0), while the absorptive contribution is pro- 
duced when the soft wide-angle gluon is slightly off-shell (g^ ~ — q^^, where is the 
gluon transverse momentum with respect to the direction of the momenta of the pair of 
interacting external partons). 

The colour coherence phenomenon discussed so far in this subsection refers to the ra- 
diative part of the one-loop interactions. The absorptive (imaginary) part of the one-loop 
wide-angle interaction between the non-collinear parton pj and the collinear parton pi is 
very similar to its radiative part (it is still proportional to Tj ■ Ti and it simply depends 
on 6jp in the collinear limit), but it is non- vanishing only if sji > (as recalled below, this 
constraint follows from causality). Combining the absorptive part of the interactions of pj 
with the two collinear partons (as we did previously by combining the radiative part), we 
obtain a contribution that is proportional to Tj ■ [Ti Q{sji) + T2 ©(5^2) ], and we see that 
the two collinear partons can act coherently only if their energies pi and P2 have the same 
sign (i.e., in the case of the TL collinear limit). In the absence of this coherent action, the 
collinear splitting subprocess retains absorptive interactions with the non-collinear partons, 
and strict collinear factorization is violated in the SL collinear limit. In the generalized 
factorization formula the amplitude \Ai^^\P, . . . ,Pn)) includes the absorptive part of 
the interactions with the parent parton P (namely, the terms Tj ■ Tp Q{sjp)), and the 
remaining absorptive contribution effectively included in Sp^^^ is proportional to 



T 



T 



Ti e(s,i) + T2 0(s,2) \-T,- Tp Q{s^p) 

Ti (e(.,o - Qis^p)) + T2 (e(s,2) - e(.^.^)) 



If, for example, we consider the SL case with Z2 < (thus, Zi > 0), the energies s^i and 
Sjp have the same sign, while the energies Sj2 and Sjp have opposite sign; the absorptive 

contribution to Sp*-^^ is thus proportional to Tj ■ T2 ^0(sj2) — ©(-Sj-p) j = Tj ■ T2 sign(sj2), 
which is exactly the colour-charge factor that appears in the right-hand side of Eq. flS^ . 
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The one-loop interaction between the partons j and i has a non-vanishing absorptive 
part only if the parton energies p° and have the same sign (i.e. only if sji > 0). This 
requires that the two partons belong to either the physical initial state or the physical 
final state of the scattering amplitude. In other words, the absorptive part has a definite 
causal structure (and origin): it arises from long-range interactions that takes place at 
large asymptotic times t — ±00, either in the past (considering initial-state partons) or in 
the future (considering final-state partons), with respect to the short time intervals that 
control the hard-scattering subprocess. 

In gauge theories, factorization is potentially spoiled by the long-range nature of gauge 
interactions. Colour coherence can effectively restore the space-time factorization of small- 
distance and large-distance subprocesses. Colour coherence argument requires no distinc- 
tions between large space distances and large time distances. Causality does make distinc- 
tions between large distances at t — )■ —00 and t — t- +00. Therefore, if the large-distance 
subprocess involves interactions at both t — t- —00 and t — )■ +00 (as is the case of the two- 
parton collinear splitting in the SL region), the factorization power of colour coherence is 
limited by causality. This limitation leads to violation of 'strict' factorization. 

The absorptive part of the one-loop splitting matrix Sp^^'^ in Eq. (155]) is responsible for 
violation of strict factorization in the SL collinear region. We observe that this absorptive 
part is IR divergent. Its IR divergent contribution (see e.g. Eq. fl35l) and remember that 
/(e; z — iO) = — \n{z — iO) + C(e)) is proportional to —in Tj ■ T2/e, which exactly corre- 
sponds to the contribution of the non-abelian analogue of the QED Coulomb phase. It is 
well known (see, e.g., Refs. [371[3HlEniE3lEDlEIlE21linillIlll2| and references therein) 
that 'Coulomb gluons' lead to non-trivial QCD effects, especially in relation to various 
factorization issues and in connection with resummations of logarithmically-enhanced ra- 
diative corrections. Our study of the one-loop splitting matrix Sp^^\pi,p2; P',P3, ■ ■ ■ ,Pn) 
(including its absorptive part) is performed to all orders in the e expansion. We thus note 
that our results on the SL collinear limit of two partons are not limited to the treatment 
of Coulomb gluon effects to leading IR (or logarithmic) accuracy. 



5 Multiparton collinear limit and generalized factor- 
ization at one-loop order and beyond 

5.1 Multiparton collinear limit of tree-level and one-loop ampli- 
tudes 

The definition of the collinear limit of a set {pi, . . . ,Pm} of m (m > 3) parton momenta 
follows the corresponding definition of the two-parton collinear limit (see Eqs. ([2])-([5])). 
The multiparton collinear limit is approached when the momenta of the m partons become 
parallel. This implies that all the parton subenergies 

Sie = {Pi+PeY , with i, i e C = {1, . . . ,m} , (69) 



26 



are of the same order and vanish simultaneously. In analogy with the two-parton colhnear 
hmit, we introduce the hght-hke momentum^ P^: 

In the multiparton collinear limit, the vector approaches the collinear direction and we 
have — )■ ZiP'^, where the longitudinal-momentum fractions Zi are 

Zi = ^ = ■ ■ , z G C , (71) 

n- P n ■ (pi H \-pm) 

and they fulfill the constraint Yl^i -^^i = 1- More details on the kinematics of the multipar- 
ton collinear limit can be found in Ref. [TT| . 

As in the case of the two-parton collinear limit, the dynamics of the multiparton collinear 
limit of scattering amplitudes is different in the TL and SL regions. The TL region is 
specified by the constraints su = {pi + Pi^ > 0, where {pi,Pi} refers to any pair of 
collinear-parton momenta; note that these contraints imply 1 > Zi > 0. If these contraints 
are not fulfilled, we are dealing with the SL region. 

According to this definition, in the TL case, all the partons in the collinear set are either 
final-state partons (i.e., physically outgoing partons with 'energies' p° > 0) or initial-state 
partons (i.e., physically incoming partons with 'energies' p^ < 0). In the SL case, the 
collinear set involves at least one parton in the initial state and, necessarily, one or more 
partons in the final state. 

In the kinematical configuration where the m parton momenta pi, . . . ,Pm become simul- 
taneously parallel, the matrix element M. = A4{pi, ■ ■ ■ ,Pm,Pm+i, ■ ■ ■ ,Pn) becomes singular. 
The dominant singular behaviour is M{pi, ■ ■ ■ ,Pm,Pm+i, ■ ■ ■ ,Pn) ~ (l/x/s)"*"^ mod (In*^ s) 
(see Eq. ([6]) for comparison), where the logarithmic enhancement is due to scaling violation 
that occurs through loop radiative corrections. Here s denotes a generic two-parton (i.e., 
Sii) or multiparton (e.g., (^1+^2+^3)^) subenergy of the system of the m collinear partons. 

The extension of the coUinear-factorization formulae ([9]), f l2T|) and ( 133|) to the case of 
the multiparton collinear limit is 

|Al(°))^5p(°)(Pi,---,Pr„;^) , (72) 

1^(1)) ^ Sp^'\p,,...,p^;P;p^^u...,Pn) \M^'^) 

+ 5p(°)(pi,...,p„;P) \M^') . (73) 

On the right-hand side of Eqs. (172!) and (!73|) . we have neglected contributions that are sub- 
dominant (though, still singular) in the collinear limit, and we have denoted the 'reduced' 
matrix element by introducing the following shorthand notation: 

M=M{P,Pm+l,...,Pn) . (74) 



^To be precise, a more appropriate notation should be used to distinguish the two vectors in Eqs. ([5]) 
and ([70)) . We do not introduce such a distinction, since we always use Eqs. ([S]) and ([70)1 in connection with 
the corresponding collinear limit of 2 and m partons, respectively. 
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The reduced matrix element Ai is obtained from the original matrix element M. by re- 
placing the m collinear partons Ai, . . . , Am (whose momenta are pi, ■ ■ ■ ,Pm) with a single 
parton A, which carries the momentum P. The flavour of the parent parton A is determined 
by flavour conservation in the splitting subprocess A ^ Ai . . . Am- 

The process dependence of the tree-level factorization formula (!72|) is entirely embod- 
ied in the matrix elements Ai and A4. The splitting matrix Sp^'^\pi, . . . , pm', P) , which 
captures the dominant singular behaviour in the multiparton collinear limit, is universal 
(process independent). It depends on the momenta and quantum numbers (flavour, spin, 
colour) of the sole partons that are involved in the collinear splitting A ^ Ai . . . Am- The 
colour dependence can explicitly be denoted as (see Eq. ( ITOl) in the case of m = 2 collinear 
partons) 

^^(0) (c„...,c™;c)^^^^ ...,Pm;P)^ (Ci, . . . , Cm\ Sp^'\p,, - - - , Pm', P) \c) , (75) 

where ci, . . . , Cm are the colour indices of the partons Ai, - ■ ■ , Am, and c is the colour index 
of the parent parton A. 

At the tree level, the TL and SL collinear limits have the same structure and are 
simply related by crossing symmetry relations. In going from the TL to the SL regions, 
the multiparton splitting matrix Sp^^\pi, - - - ,Pm] P) in Eq. f l72|) only varies because of 
the replacement of the wave function factors of the collinear partons (see the final part of 
Sect. [2]). The dependence of Sp^^^ on the colours and momenta of the collinear partons is 
unchanged in the TL and SL regions. 

At the one-loop order, the TL and SL collinear limits have a different structure. As 
a consequence of the violation of strict collinear factorization, Eq. f l73|) is presented in 
the form of generalized collinear factorization (see Eq. f p3|) ). In the SL collinear region, 
the multiparton one-loop splitting matrix Sp^^^ also depends on the momenta and colour 
charges of the non-collinear partons in the matrix elements Ai and Ai- Introducing the 
colour dependence in explicit form, we have 

^^(1) icu...,c„;c,c'^^„...,c'„)^p^^ .,Pm-P;Pm+l, - - - , Pn) 

= (ci,...,c„| Sp^'^\pi,---,Pm;P;Pm+i,---,Pn) |c, C^+i, . . . , , (76) 

where Ci, . . . , c„ are the colour indices of all (collinear and non-collinear) the external par- 
tons in the original matrix element Ai, and c, c^+i, • • • , are the colour indices of all the 
external partons (i.e. the parent collinear parton and the non-collinear partons) in the re- 
duced matrix element At- We remark that Sp'^^^ does not depend on the spin polarization 
states of the non-collinear partons, since the violation of strict collinear factorization origi- 
nates from soft interactions between the non-collinear and collinear partons (see Sect. 14. 5p . 
This origin of the violation of strict collinear factorization also implies that the one-loop 
multiparton splitting matrix Sp^^^ has factorization breaking terms with a linear depen- 
dence on the colour matrix of the non-collinear partons (see Eq. fj^5l) and also Sect. 15. 3p . 
In the TL collinear region [281 126] strict collinear factorization is recovered, and Sp*-^-* is 
universal (i.e., independent of the non-collinear partons); thus, we can write: 

Sp^^\pi,-.-,Pm;P;Pm+i,---,Pn) = Sp^^\pi,---,Pm;P) , (TL coll. lim.) , (77) 
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where, precisely speaking, the notation means that Sp^ ' is proportional to the unit matrix 
in the colour subspace of the non-coUinear partons. 

As recalled in Sect. [21 coUinear factorization of QCD amplitudes is usually presented 
in a colour-stripped formulation, which is obtained upon decomposition of the matrix el- 
ement Ai in colour subamplitudes. In this formulation, the singular behaviour of the 
colour subamplitudes in the multiparton collinear limit is described by splitting ampli- 
tudes Split^^\pi, . . . ,Pm', P) (at the tree level) and Split*-^^ (at the one-loop order). The 
splitting amplitudes can be regarded as matrices in helicity space, since they depend on 
the helicity states of the collinear partons. The splitting matrix Sp in Eqs. ( !72|) and (173|) 
is a generalization of the splitting amplitude, since Sp describes collinear factorization in 
colour -fhelicity space. In the case of the two-parton collinear limit, there is a straightfor- 
ward direct proportionality (through a single colour matrix) between the splitting matrix 
Sp^*^-* and the splitting amplitude Split (see Eq. f|T6l) ). Considering the collinear limit 
of m partons, with m > 3, the corresponding splitting matrix Sp*-"-* can get contributions 
from different colour structures. In general, Sp*-"-* can be expressed as a linear combina- 
tion of colour-matrix factors whose coefficients are kinematical splitting amplitudes Split^^-*. 
Equivalently, the splitting amplitudes Split "-^^ can be regarded (defined) as colour-stripped 
components of the splitting matrix Sp^^^ . This correspondence between Sp and Split func- 
tions extends from the tree-level to one-loop order in the TL collinear region. One-loop 
splitting amplitudes can be introduced also in the SL region (see Appendix [X]), by properly 
taking into account the violation of strict collinear factorization and the ensuing colour 
entanglement between collinear and non-collinear partons. 

At the tree level, the multiparton (m > 3) collinear limit has been extensively studied in 
the literature. In the case of m = 3 collinear partons, explicit results for all QCD splitting 
processes A — )• ^1^42^3 are known for both squared amplitudes [18l[T9l[TTj and amplitudes 
[201 [211 [22] . At the amplitude level, the multiparton collinear limit is explicitly known also 
in the cases of m = 4 [201 121] and m = 5 and 6 [21] gluons. Considering some specific 
classes of helicity configurations of the collinear partons, the authors of Refs. [211 122] derived 
general results for the splitting amplitude Split*-"-* that are valid for an arbitrary number m 
of gluons and of gluons plus up to four fermions (g, qq, qqq, qqqq). 

At the one-loop level, the multiparton collinear limit in the TL region was studied 
in Ref. [26]: we explicitly computed the one- loop splitting matrix for the triple collinear 
limit q — )■ {l ^"^^ l' quarks with different flavour), and we presented the general 

structure of the IR and ultraviolet divergences of the one-loop sphtting matrices. The latter 
result is recalled in Sect. 15.31 where it is also extended to the SL collinear region. 

5.2 Generalized collinear factorization at all orders 

The TL collinear limit of all-order QCD amplitudes is studied in Ref. pS], by using the 
unitarity sewing method [231 [13] . The detailed analysis of Ref. [2B] is limited to the leading- 
colour terms, but it can be extended to subleading-colour contributions as shown by the 
explicit computations of splitting amplitudes at one-loop [22] and two-loop [22] orders. The 
collinear limit can be studied by using a different method [11] that relies on the factorization 
properties of Feynman diagrams in light-like axial gauges. By exploiting colour- coherence 
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of QCD radiation, this method, which directly applies in colour space, can be extended 
from tree- level to one- loop amplitudes in the TL collinear region. In Eq. ( I78|) . we propose 
an all-order factorization formula that is valid in both the TL and SL collinear regions. 
In the TL collinear limit, Eq. (!78|) represents a colour-space generalization of the all-order 
results of Ref. [2B] or, similarly, an all-order generalization of the colour-space factorization 
of Refs. pT| • The extension from the TL to the SL collinear regions is based on the 
generalized factorization structure (and the physical insight) that we have uncovered at the 
one-loop order (see Sects. I4l and ISTTl) . In Sect. [6], we show that Eq. fl78|) is consistent with 
the all-order factorization structure of the IR divergences of QCD amplitudes. 

The generalized factorization formula for the multiparton collinear limit of the all-order 
matrix element A4 is 

\M) - Sp{pi,...,Pm;P;Pm+i,---,Pn) \M) , (78) 

where Sp is the all-order splitting matrix. The loop expansion of the unrenormalized 
splitting matrix is 

Sp{pi,...,Pm;P;Pm+l,---,Pn) = Sp^°\pi, . . . , p^] P) 

+ Sp^^\pi,...,Pm;P;Pm+l,---,Pn) + Sp^'^\pi,...,Pm;P;Pm+l,...,Pn) + ... , (79) 

where Sp^^^' is its one-loop contribution, Sp^"^^ is the two-loop splitting matrix, and so 
forth. The loop expansion of the matrix element Ai is defined in Eq. ([1]), and an analogous 

expression (it is obtained by simply replacing '^'^^ with Ai ) applies to the reduced matrix 
element Ai. Inserting these expansions in Eq. f l78p . we obtain factorization formulae that 
are valid order-by-order in the number of loops. At the tree level and one-loop order we 
recover Eqs. fl72l) and fl73|l . respectively. The explicit factorization formula for the two-loop 
matrix element Ai^'^^ is 

|A^(2)) ^ Sp^'\p,,...,p^;P;p^^u..-,Pn) \M^') 
+ Sp^^\pi,...,p^;P;p^+i,...,Pn) l-M^^^) 
+ Sp^'\p,,...,prn;P) \M^'^) . (80) 

In the case of the TL collinear limit, strict factorization is valid and the splitting matrix 
is process independent at each order in the loop expansion; we can remove the reference to 
the non-coUinear partons, and we can simply write 

Sp{pi, . . . ,pm;P;Pm+i, ■ ■ ■ ,Pn) = Sp{pi, . . . ,pra;P) , (TL coll. lim.) . (81) 

Beyond the tree level, in the case of the SL collinear limit, the splitting matrix Sp acquires 
also a dependence on the external non- collinear partons of At and At, although it is 
(expected to be) independent of the spin polarization states of these partons. 

The all-order matrix elements Ai and Ai in Eq. f lTHj) are invariant under the renormal- 
ization procedure, since they are scattering amplitudes whose external partons are on-shell 
and with physical polarization states. Therefore, the splitting matrix Sp (the remaining 
ingredient in Eq. (EH])) is also invariant. The renormalization of A^,A^ and Sp simply 
amounts to replace the bare coupling constant with the renormalized QCD coupling. The 
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perturbative (loop) expansion with respect to the renormahzed couphng is denoted as fol- 
lows: 

>l = 7W(0'«)+7W(i.«)+7W(2.R) + ... , (82) 

Sp = + + + . . . , (83) 

where the superscripts {k,R) (A; = 0, 1, 2, ... ) refer to the renormalized expansions, whereas 
the superscripts (k) refer to the unrenormalized expansions in the corresponding Eqs. ([T]) 
and f l79|) . Since perturbative renormalization commutes with the coUinear limit, the per- 
turbative factorization formulae fl72|) . fl73|) and f lHOj) can equivalently be written by using 
the renormalized expansion. We simply have: 

|_^{o,K)^ ~ , (84) 

^ ITW^^'^'V ^P^^'^'MA^^''''^) , (85) 



The main features of the results that we present in the following sections do not depend 
on the specific renormalization procedure. To avoid possible related ambiguities, we specify 
the renormalization (and regularization) procedure that we actually use in writing explicit 
expressions. The unrenormalized quantities are evaluated in c? = 4 — 2e dimensions by 
usin^ conventional dimensional regularization (CDR); in particular, the parton momenta 
are d dimensional, and the gluon has — 2 = 2 — 2e physical polarization states. The 
renormalized QCD coupling at the scale fiR is denoted by gsif^n) (as = gs/i^'^))^ 
is obtained from the unrenormalized (bare) coupling gs by a modified minimal subtraction 
(MS) procedure. We use the following explicit relation: 



2e 2 o 2e 2/ 2 



1 - 



as(/i 



271 
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l-l)+(^(«^(/^^)) 



^7) 



where 5"^ is the customary MS factor {je = — V^(l) = 0.5772 ... is the Euler number), 

5, = exp[e(ln47r + V'(l))] , (88) 
and bo (see Eq. ( I32l) ) and bi are the first two coefficients of the QCD /3 function, with 

b^ = {17 Cl - 5 CANf - 3 CFNf)/Q . (89) 



In the following, all the renormalized expressions refer to the perturbative expansion 
with respect to as(/i^) (i.e., the renormalization scale fiR is always set to be equal to the 
dimensional-regularization scale fi). For instance, considering the splitting matrix Sp of 
m collinear partons, the tree-level and one-loop relations between the renormalized (see 
Eq. (183|1 ) and unrenormalized (see Eq. ( !79|l ) contributions are as follows: 



§The relation between different regularization schemes within dimensional regularization can be found 
in Refs. [Ml Sg |46] . 



Sp 



(0) 



(90) 
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27r 



+ 



Sp 



(1) 



(91) 



The first term on the right-hand side of Eq. flQljl originates from the fact that the corre- 
sponding tree-level contribution Sp^'^^ is proportional to gl^~^- 

In the case of the two-parton collinear limit, for later purposes, we can write the renor- 
malized version of Eq. fl34p in the following form: 



Sp^''^\p^,P2■, P;P3, ...,Pn) = Ici^) Sp^''^^\p,,p^- P) + 5p^'"^(pi,p2; P) 



(92) 



Ic (e) 



Qs(/i^) 1 
27r 2 



cr 



-Si2 - «0 



X 



I ^ (C'i2 - Ci - C2) + ^ {^12 - 7i - 72) 



— z 



C12 + Ci - C2 j /,j(e; ^1) + (C12 + - Ci ) /^(e; ^2 

n 

Y^H^r sign(si,) e(-2;,) /,(e; z,; 

j=3 i=l,2 



(93) 



where the coefficient cp = 1 -l- ^(e^) is related to the volume factor cr in Eq. (!30|) and it is 
defined as 

r(i + 6)r2(i-e) 



Cr = (47r)' Cr 



r(l-2e) 



(94) 



The colour operator I^J i^) directly related to the corresponding unrenormalized 
operator Ic{pi,P2',P3, ■ ■ ■ ,Pn) (the functional dependence on the parton momenta is not 

explicitly recalled in the argument of Ic'^) in Eq. fH5|) . In the expression on the right- 
hand side of Eq. fl93l) . the colour correlation contributions proportional to Tj ■ Ti are 
written as in Eqs. (l56l) -(!59l). by using the decomposition of /(e; z) in its real and imaginary 
parts. We also note that the contribution proportional to Bq in Ic and the renormalization 
counterterm of Sp^^'^^ (i.e. the first term on the right-hand side of Eq. ( 1911) ) have been 
included in the definition of Sp^^'^K Thus, in Eq. ( |92l) . Sp^^'^^ has been defined as 



Sp^i'''\p,,P2;P) 



+ 



2tt 2e 
,(1) 



Cr 



-S12 - iO 



- 1 



5p(°'^)(pi,P2;P) 



Sp'^>{p^,P2;P) 



(95) 



The first term on the right-hand side is finite when e — )■ 0, since the ultraviolet divergence 
of Ic is cancelled by the renormalization of the one-loop splitting matrix. The remaining 



e poles of Sp^^'^^ are of IR origin, and they are included in the colour operator Ic'{^) 



(1), 



5.3 The IR divergences of the one-loop spUtting matrix 



Analogously to the case of two-collinear partons, the one-loop splitting matrix Sp^^'^^ of 
the multiparton collinear limit in Eq. (185|) (or Eq. ( !73|) ) has IR divergences. They show up 
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as double (1/e^) and single (1/e) poles in the expansion around the point e = 0. To make 
the IR behaviour explicit, we separate the divergent and finite terms as follows: 

= Sp^^^^"- + , (96) 

where Sp^^^''^"' contains the e poles, while Sp^^'^^'^' is finite when e — )■ 0. Then, we can 
present our general result for the IR divergent part. 

The computation of the divergent part shows that it can be written in a colour-space 
factorized form: 

Sp^'^''"ipU...,Pm;P;Pm+U...,Pn)=ltU^) Sp^''^\p^, . . . , p^; P) , (97) 

where S'p*-^'^'' is the tree-level splitting matrix, including its complete (i.e. exact) depen- 
dence on e. The colour space operator J^-*^ depends on the colhnear partons and, in the 
SL region, it also depends on the momenta and colour charges of the non-collinear partons 
in the original matrix element M. Analogously to the notation in Eq. fl92|) . the functional 
dependence on all partons is not explicitly denoted in the argument of /^c- present 
I^^Q in a general form, which is valid in both the SL and TL collinear limits. Considering 
the colhnear splitting process A ^ Ai . . . Am, the explicit expression of 1^^^ is 



^ i n V A'' / ^ ■ ^r< I 



where the subscript P refers to the parent collinear parton A. The sets 

C = {l,...,m} and NC = {m + 1, . . . ,n} 
denote the collinear and non-collinear partons, respectively. 

The expression (!98|l can be written in the following equivalent form: 



r{l) 



(1) 

*-mC 



[e) , (99) 



where 



A;^^^) = ^ - T, e(-z.) sign(.,,) . (100) 

^ i&C 
jdNC 



^The double pole in /[^'p interferes with the term of 0{e) in Sp^'^''^\ thus contributing to the 
single pole 1/e in Sp^^^'^"'. 
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The result in Eq. f PS]) (or Eq. fl^ ) contains double and single poles in 1/e. For the 
specific case of m = 2 collinear partons, this result agrees with the e poles of the complete 
one- loop result in Eqs. (!92|) and ( 193|) . In the general case of m collinear partons, we can 
see that I^^c embodies colour-correlation terms that produce violation of strict collinear 
factorization. In the expression fl99|) . these correlations are fully taken into account by the 
colour operator A^^(^(e). The expression fllOOp exphcitly shows that the operator A^^^^^ 
is antihermitian. In particular, A^''^(e) is proportional to a single pole 1/e and exactly 
corresponds to the colour-correlation terms that are produced by the non-abelian Coulomb 
phase (see Sect. 14.51) . 

The computation of the finite part of Eq. fl96p . in the case of the SL collinear limit of 
m = 3 partons, is in progress. In general, we can anticipate that Sp^^^^^' also contains 
factorization breaking contributions, and they are of the same type as those in A^^^^, 
namely, they have the form of two-parton correlations Tj ■ Tj with i & C and j G NC. 

In the TL collinear region, the longitudinal-momentum fractions Zi are positive. There- 
fore, A^''^ vanishes and the result in Eq. has a strict-factorization form, which agrees 
with the result in Eq. (11) of Ref. [SB]. To be precise, Eq. (11) of Ref. contains an 
additional single pole of UV origin (that expression refers to the unrenormalized splitting 
matrix) and additional IR finite term^, which include the 0(e'^) terms of Sp^^^ that ex- 
plicitly depend on the regularization scheme. In the TL collinear region, the expression 
(1^^ can also be rewritten in a slightly-simpler form, as follows: 




We can see that the antihermitian part of is simply an imaginary c-number. Therefore, 
in the TL collinear limit, the non-abelian Coulomb phase effectively takes an abelian form. 

The derivation of the results in Eqs. flU7|) and flUS|) is presented below. As for the 
rewriting in Eqs. (!99|) and fllOip . it simply follows from colour conservation properties. 
Considering the last term on the right-hand side of Eq. ( |98|) . we can write ln(zj — iOsij) = 
In \zi\ — iirQl—Zi) sign(sjj) and then, using the conservation of the total colour charge, 

E = - E (102) 

j £NC i ec 

we directly obtain Eqs. fl^^ and fllOOl) . Considering the last term in the curly bracket on 
the right-hand side of Eq. (!99|) . we can write ln(— Sj£ — iO) = Insj^ — m {su > in the TL 
collinear region), and then Eq. fllOip is directly obtained by using Xlj £ec Ti-T^ = Cp or, 

"The separation in Eq. (|96p is, in part, a matter of definition, since some IR finite terms can always be 
moved from Sp^^^^"' to Sp'^)^'^'. 
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more precisely, by using the relation 



( 5^ T, ) = Tj5 . (103) 

iec 

This relation (which generalizes Eq. fHOj) to the case of m collinear partons) follows from 
the fact that the colour charge Tp of the parent parton A is conserved in the tree-level 
collinear splitting A ^ Ai . . . Am- 

Note that, in the case of the TL collinear splitting, Eq. f ll03p remains valid when 
replacing the tree-level splitting matrix Sp^^'^^ with the all-order splitting matrix Sp. We 
remark that this all-order generalization does not apply to the SL collinear limit. In the 
SL case, owing to the violation of strict factorization, Sp^^\ Sp^'^^ and so forth, contain 
charge interactions between collinear and non-coUinear partons: therefore, only their total 
charge is conserved (and the equality in (11031) is not valid beyond the tree level). 

The results in Eqs. ( 197|1 and ( p8|l can be derived in a simple way by exploiting the known 
IR structure of generic one-loop QCD amplitudes [HI ESI [I]- Considering the scattering 
amplitude A4 with n external QCD partons (and any number of colourless external legs), 
the renormalized one- loop contribution can be written in the following factorized form [1] : 

17^(1'^)) = /«(e) + , (104) 

where 

1=1 ^ ' «,i=i 

and the one-loop term |A^*^^)^'^ ) is finite when e — )• 0. The singular dependence on e 
is embodied in the factor /^^""(e) that acts as a colour-charge operator onto the tree- level 
matrix element (which retains its complete dependence on e). The colour operator 

/^■* is equivalent, though it is not exactly equal, to the operator J*-^-* introduced in Ref. [1]; 
the expression (11051) is obtained from /*^^^ by removing all terms that are finite at e = 
(these terms are absorbed in the definition of lA^*^^)*^"^ )). 

The universality structure of the IR factorization formulae (11041) and (11051) has direct 
consequences on the collinear limit of the scattering amplitudes. Indeed, Eq. (11041) can be 
applied to the matrix element M. before performing the collinear limit, and an analogous 
relation (which is obtained by the replacements M. ^ M. and /^^■' — )■ ) applies to the 
reduced matrix element M.. Therefore, we have: 

IXT^O'^)) ^ - lAY^'"''^) (106) 

= /iy(e) - /(i)(e) + lA^(i)fi-) - 5p(°'^) (107) 

^ {l^{l{e) - J^)(6)) + - |7W^'^'°-) . (108) 

Equation ( I106p is just a rewriting of the one- loop collinear factorization in Eq. ( 185|) . Then, 
Eq. (flOTl) is obtained by using Eq. (fTOl]) for both M^^^^^ and 'J^^'^\ Finally, Eq. ( 1T08|) is 
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obtained by applying the tree-level factorization formula f l84p to Performing these 

steps, we have only neglected subdominant collinear terms (as denoted by the approximate 

equalities in Eq. (1T06|) and (fT08|) ). The last terms, IM^^^^""-) and Sp^°'^^ |AY^^^^°'), on the 
right-hand side of Eq. fllOSp do not contain any e poles, and they are finite if e — )■ 0. Thus, 
in the collinear limit, we can write 

5p(^'^) ^ Jiy(e) - 5p(°'^) I^^ie) + 0(e°) . (109) 



The IR divergent contributions in Eqs. fPTI) and fl^Sl) directly derive from the relation 
(11091) . by simply using colour-charge conservation. To illustrate the derivation in detail, 
we note that, according to Eq. (11051) . the operator I^l^ on the right-hand side of Eq. (11091) 
contains three classes of contributions: (a) terms that only depend on the non-coUinear 
partons; (6) terms that only depend on the collinear partons in Ai; (c) terms that depend on 
both the collinear and non-coUinear partons. Correspondingly, the operator contains 
the following classes of contributions: (a) terms that only depend on the non-coUinear 
partons; (b) terms that only depend on the parent collinear parton A in Ai; (c) terms that 
depend on both the parent parton and the non-collinear partons. The terms of the classes 
(a) and (a) are exactly equal and commute with Sp^^'^^: therefore, they cancel on the 
right-hand side of Eq. (I109p . The terms of the class (6) are proportional to either Cp or 
7p; they commute with Sp^^'^^ and combine with the terms of the class (6), thus leading to 
the first three contributions in the curly bracket of Eq. (I98p . The last contribution (which 
depends on Tj ■ Tj) in the curly bracket of Eq. (1^5]) is due to the terms of the remaining 
classes (c) and (c). Indeed, these terms give the following contribution to the right-hand 
side of Eq. (I109P (we omit an overall factor of as/(4vr) and use the notation Sjp = 2pj ■ P): 



E T,-T,ln 



jGNC 



jeNC 

2 , f-Sjp-iO 



jGNC 



Sp(^'^) + 5piu,Kj z J2 Tj- Tp In 



~ - - ^ T,- ■ T, m {Zi - zOs.i) Sp 



(o,R) 



iGC 



-Sjp - iO 



(110) 



(111) 



In going from Eq. (1 11 OP to Eq. (Illip . we have performed the collinear limit Pi ~ ZiP 
{i E C), and we have thus used the following coUinear approximation: 



1 Sji iO 



In 



-ZiSjp - iO 



~ In {zi — iOsji) + In 



In {zi + iOSjp ] + In 



s -5 — iO 



-Sjp - iO 



(112) 



The term in the square bracket of Eq. f illip vanishes (we recall that the non-collinear 
charge Tj commutes with Sp*-"^'^'^) because of the conservation of the collinear charge Tp 
of the parent parton (see Eq. (I103P ). Therefore, the expression (II lip exactly corresponds 
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to the contribution to Eq. fl97p from the last term in the curly bracket of Eq. 0981] . This 
completes the derivation of the results in Eqs. (197|) and (p8ll . 

The structure of Eqs. f l96|) -f l98|) and its derivation agree with the discussion in Sect. 14.51 
In particular, the pairwise interaction terms of the classes (c) (i.e. the first contribution in 
Eq. f lllOp ) and (c) (i.e. the second contribution in Eq. f lllOp ) include the IR divergent part 
of the one-loop contributions due to the soft and wide-angle region of the loop momentum. 
The terms of the class (c) correspond to one-loop contributions to Ai that have been 
removed from the one-loop splitting matrix Sp^^'^^ and re-factorized in A4 (see Eqs. (11060 
and (llOQp ). These terms, which are rewritten as second contribution in the square bracket 
of Eq. (imp , are effectively equal to part of the contribution of the class (c): this part (i.e. 
the first contribution in the square bracket of Eq. fillip ) is due to the one- loop interactions 
at wide angle with respect to the direction of the system of collinear partons in Ai. The 
vanishing of the square bracket term in Eq. (Illip is produced by the coherent action of the 
collinear-parton system with respect to non-factorizable wide-angle interactions with each of 
the non-collinear partons. This colour-coherence mechanism guarantees strict factorization 
in the TL collinear limit. The remaining part of the terms of the class (c) (i.e. the first 
term in Eq. (Illip ) includes an absorptive (imaginary) component and a radiative (real) 
component. Owing to their causality structure, these remaining absorptive contributions 
are non- vanishing only in the SL collinear limit (i.e., Zi < 0): the absorptive part of the 
one-loop interactions between collinear and non-collinear partons produces the violation of 
strict collinear factorization. The remaining radiative contributions do not violate strict 
collinear factorization, since the non-collinear partons act coherently as a single parton, 
whose colour charge is equal (modulo the overall sign) to the total colour charge of the 
collinear partons. 

To conclude this section, we rewrite the steps in Eqs. (I109p - (lll2p by using a slightly 
different and more compact notation. This rewriting anticipates the notation that we use in 
some of the following sections (e.g.. Sect. 16. H and Appendix[B]). We consider the right-hand 
side of Eq. f ll09p , and we write it in the following form: 

/2(e) Sp^'^""^ - 4He) + 0(e°) ^ I^'\e) - l^'\e) (113) 



The right-hand side of Eq. flll3p is obtained by neglecting IR finite contributions of 0{e 

performing 
notation is 



0^ 



performing the collinear limit of the operator and defining the operator J*^^^; our 

/j;(e) ^ /«(e) , (s,, ~ 2 e C, J G NC) , (115) 



, ^ ^ = ■ (116) 

The operator J^^-* in Eq. (11151) is obtained from I^^ in Eq. (11051) by implementing the 
collinear approximation pi ^ ZiP {i G C) in all the terms of I^^ that are not singular in 
the collinear limit. The explicit expression of I^^\€) is derived by using collinear relations 
as in Eqs. (11110 and (I112p . and it is presented in Eq. (I225P of Appendix [Bl The IR operator 
/il^ of the reduced matrix element A4 depends on the colour charges T,- of the non-collinear 
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partons and on the colour charge Tp of the parent colhnear parton. The definition of I 
in Eq. f lll6p simply amounts to the implementation of the colour conservation relation 

Tp = — J^jeNc'^j- "^^^ operator J ^ is obtained from by replacing the colour matrix 
—Tp with the sum of the non-collinear charges. This replacement produces the explicit 

expression of I^^^ that is presented in Eq. fl226p of Appendix [Bl According to this definition, 
the matrix structure of /*^^^ only depends on the colour charges of the non-collinear partons: 
therefore, I^^^ commutes with Sp^"'^^ (since Sp^^'^^ does not depend on the non-collinear 
partons), and the commutation leads to Eq. fllMp . The IR divergent part Sp^^^^^^' of the 
one-loop splitting matrix Sp^^'^^ in Eq. ( |96l) is obtained by performing the colhnear limit 
of Eq. f llOQp and removing the IR finite terms of 0{e^). These steps are carried out in 
Eqs. flllSp and flll4p . Equating the relations fl97|) and f lll4p . we thus obtain the following 
representation of the IR divergent operator /^^c*- 

I^U^) = ■ (117) 

The explicit expressions fl225D and fl226|) of J^^^ and F ' can be inserted in Eq. flll7p . and 
we can directly check that we reobtain the result in Eq. fl^Sj) . 



6 All-order structure of the collinear limit and 
space-like two-loop results 

6.1 The IR structure of the spUtting matrix 



The IR structure of multiloop QCD amplitudes is not independent of their collinear be- 
haviour. This fact can be exploited to extract non-trivial information [26| 1271 13 E]. In 
Sect. 15.31 the IR divergences of the one- loop splitting matrix Sp^^'^^ were extracted from 
the known IR structure of the one-loop QCD amplitudes. In this section, our study of the 
collinear limit is extended beyond the one-loop order. 

We consider a generic scattering amplitude Ai with n external QCD partons (and any 
number of colourless external legs). The IR structure of Ai at two-loop order is known in 
explicit form, and it is given by the following colour-space factorization formula [Ij: 



1^(2,/?)) 



r(l) 



{i,R)\ 



,(2) fin/ 



:ii8) 



where the contributions Ai^'''^^ [k 
scattering amplitude (see Eq. 



0,1,2) refer to the renormalized expansion of the 
). The two-loop term is finite when e — )• 0. 

The colour-charge operator is the IR divergent operator that controls the one-loop 
factorization formula fll04p . The IR divergent two- loop operator IjJ can be written in 



terms of in the following form [H [3] : 
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Note that I^l) appears on the right-hand side with two different arguments, namely, (e) 
and I^^{2e). The expression ( I119p includes IR poles of the type 1/e", with n = 4,3,2, 1. 
The dominant pole terms, and are fully controlled by I^^^ . The double-pole 

terms, also depend on the value of the coefficient K [I]: 

The control of the single-pole terms, 1/e, also requires the knowledge of the c-number 
coefficients H^"^"^ [3J, which depend on the flavour i (quark, antiquark or gluon) of the parton 
with momentum pi. The value of can be extracted [U |3] from explicit computations 
of two-loop amplitudes; the quark (antiquark) coefficient Hq^^ is [31 [H US] 

q q 4 \^ 2 ^•^ sy ^ 8 V 108 12 y 

+ i + N,C, . (121) 

8 V 54 6 y ^ ' ^ ' 

and the gluon coefficient Hg"^^ is [21 SS] 

+ ^(^Cs + ^~^^Acl . (122) 



24 V 2 12 

The expression flllQp of the two-loop operator is similar (and equivalent) to that 

of the operator J*^^^ in Ref. P: the differences eventually amount to IR finite contributions 
to Eq. (11181) that are included in the definition of A^*^^^^"'. The essential difference is due 
to the fact that I^^i^) -^m (^) include only contributions from IR poles 1/e" (n > 1), 
with no additional IR finite contributions. The key remark of Ref. [2| is that this 'minimal 
form' of the IR divergent operators I^^^ and I^^j highly simplifies the expression of I^^^ . 
The single-pole coefficients H^'^^ in Eqs. flll9p are c-numbers [3j, whereas the form of J^^-* 
in Ref. [1] includes more complex colour-matrix correlations at 0{l/e) (see Appendix A. 3 
in Ref. [27]). We note that the coefficients i/f ^ in Eqs. ffTT9|) . ffT2TD and ffT22|) are directly 
related to the coefficients E'j^^*"-^^ in the second paper of Ref. [3J (see Eq. (3.9) therein); the 
precise relation is 

^(2)^4^W(2)^ 1^^^ . (123) 

The colour-space factorization formula (11180 can directly be exploited to extract explicit 
information (see Sects. \^72\ and 1^75]) on the collinear limit of two-loop QCD amplitudes. The 
procedure (which is illustrated in Appendix [B]) is similar to that used at one-loop order in 
Sect. 15.31 This procedure can be extended beyond the two- loop order, as discussed below. 

The one-loop and two-loop IR factorization formulae ( 11040 and (I118P can be iterated 
to higher-loop orders. The all-order generalization is 

\M) = lMie) \M)+ IM^"^-) , (124) 
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where the operator Im{^) is IR divergent, while the matrix element term Ai^"^- is IR finite. 
The perturbative (loop) expansions of Im and M.^^' in terms of the renormalized QCD 
coupling are 

iMie) = I^Sie) + /£)(e) + /£)(e) + ■ ■ ■ , (125) 

_yy^fin. ^ _yy^(0,iJ)^_yy^{l)fin. ^_^(2)fin. ^_yy^{3)fin. , (126) 

where Ai^^'^^ is the (complete) tree- level expression of the matrix element Ai. Inserting 
Eqs. ( 182|) . f ll25p and f ll26p in Eq. fll24p and performing the loop expansion, we recover 
Eqs. f ll04p and flllSp . and we obtain corresponding IR factorization formulae at three-loop 
order and higher-order levels. 

The recursive structure of Eq. f ll24p can be rewritten in the following form: 

|A^)=VM(e) , (127) 

where the all-order IR factor Vjv/(e) only depends on Im{^)'i the inverse operator V^/ 
simply is 

V,-/(e) = l-/M(e) . (128) 



We can also express the IR factor Vm in exponential form: 

VM(e) =exp|/M,cor(e)} , (129) 
where the relation between Im,cot and Im is 

expj -/M,cor(e)} = 1 - /^(e) ■ (130) 
The perturbative expansion of the operator Im,cot is 

/M,cor(6) = /g,„,(6) + /g_(e) + iflje) + ■ ■ ■ , (131) 

(k) 

and the perturbative contributions Im cor directly related to the corresponding con- 
tributions in Eq. f ll25p (in particular, Im cor — ^a^i)- Owing to the definition fll30p . the 
perturbative terms of the exponentiated operator —lM,cor give the irreducible-correlation 
component (in a statistical language) of the perturbative expansion of the operator —Im- 

In the context of the IR structure of multiparton scattering amplitudes, the exponenti- 
ated representation in Eq. f ll29p exists and is particularly suitable in view of the physical 
property of exponentiation of the IR divergent contributions to QCD scattering amplitudes 
with multiple (i.e. n > 3) external legs [50| 12]. IR exponentiation means that the dominant 
IR divergences at high perturbative order are directly captured by simply exponentiating 
the IR divergent terms that appear at lower orders. This also implies that the exponent 
function is less IR divergent than the exponential function. Using dimensional regulariza- 
tion, the dominant IR divergence of the operator Im at the n-th perturbative order (see 
Eq. ([T25D) is l\^J (e) ~ Og/e^", whereas the perturbative expansion (1131 p of the exponen- 
tiated operator lM,cor has a less singular IR behaviour of the type iM^ori^) ~ ag/e"'*'^. 
The operator lM,cor{^) has a compact all-order integral representation that is given in 
terms of a perturbatively-computable kernel of soft and coUinear anomalous dimensions 
[21 [3l m ini m [TJ [8] . We do not use this integral representation in the present paper. 
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We have briefly recalled some results on the IR structure of the QCD scattering ampli- 
tudes. These results are sufficient for the following discussion of the multiparton coUinear 
limit. Considering the collinear limit of m parton momenta {pi, . . . ,Pm} of the all-order 
matrix element Ai with n external QCD partons, we obtain 



IM^""-) = y-Jie)\M)c^V],lie) Sp\M) (132) 
= V,-/(e) Sp Yjrie) \m'''-) . (133) 



The first equality on the line fll32p is just a rewriting of the IR factorization formula 
f ll27p . Then, we have applied the collinear formula f lTHj) and, finally, in the expression 
(I133P we have used Eq. ( 1127P by replacing Ai with the reduced matrix element A4. The 
collinear limit performed in Eqs. fll32p and fll33p relates the matrix element Ai with the 

reduced matrix element A4 . Since both Ai^^- and Ai are IR finite, the colour matrix 
V^/ Sp \jj in Eq. (11331) must be IR finite in the collinear limit. The collinear limit of 
this colour matrix is denoted by Sp^"', and we can write 

5p«-^V^/(e) SpVjrie) . (134) 

Note that we have not yet implemented the collinear limit in the operator Vm- We thus 
introduce the IR divergent operators V and V as follows: 

VAf(e) - V(e) , (135) 
Vm(6)1 ^ V(e) . (136) 

The relations (I135P and (I136p represent the all-order generalization of the one-loop relations 
in Eqs. flllSp and (I116p . The relation (I135P defines V(e) through the collinear limit of 
Vjv/(e), which is obtained by using the approximation pi ~ ZiP {i = l,...,m) for all 
the terms of Vjv/ that are not singular in the collinear limit. The operator V;^ depends 
on the colour charges of the partons in Ai: the non-coUinear partons with momentum 
Pj {j G NC) and the parent collinear parton with momentum P. These colour charges 
are related by colour conservation {Tp + J^jeNc'^j ~ '^^^ relation fll36p defines V(e) 
through the implementation of colour conservation: in all the terms of Vj^ that depends on 
Tp, the colour matrix —Tp is replaced by the sum of the non-collinear charges. According 
to this definition, the colour matrix structure of V only depends on the colour charges of 
the non-collinear partons; therefore, the colour operator V has a well defined action onto 
both Ai and A4, or, equivalently, onto both the right-hand and left-hand sides of the colour 
matrices Sp and Sp^^'. 

Using the definitions f ll35p and f ll36p . Eq. f ll34p gives 



Sp = V(6) Sp^'^- V (e) , (137) 

where the all-order IR finite splitting matrix Sp^'^' has the renormalized perturbative 
expansion: 

= + + ^''■ + ... . (138) 

The relation f ll37p presents the structure of the IR divergences of the all-order splitting 
matrix Sp for the multiparton collinear limit. The IR divergences are embodied in the 
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operators V(e) and V(e). The operator V depends on the colour charges and momenta of 
the colhnear and non-colhnear partons in Ai. The operator V depends on the momentum 
P of the parent colhnear parton and on the colour charges and momenta of the non-collinear 
partons. The IR dependence on the non-collinear partons implies that, in general (actually, 
in the SL colhnear case), the splitting matrix Sp violates strict colhnear factorization. 
We also note that the IR factorization formula fll37p has a non-abelian structure (e.g., 
V5p^°-V"^ ^ Sp^""- and VV"^ ^ V"W). In the SL colhnear region, this 

structure produces strict-factorization breaking effects with distinctive non-abelian features 
(see Sects. 16.21 and 16.31) . 

We have derived Eq. fll37p by exploiting the generalized colhnear factorization in Eq. f lTHj) 
and the IR factorization property (Eq. fll24p or, equivalently, Eq. fll27p ) of the multiparton 
QCD amplitudes. According to the structure of Eqs. 01281) and fll30p . the all-order splitting 
matrix operators V and V have the following equivalent representations: 



where the all-order operators / and I (/cor and /cor) are obtained from the corresponding 
amplitude operators Im and Ij^ by using relations that are analogous to those in Eqs. f ll35p 

and f ll36p . The one-loop contribution I^^^ {I^^^) to / (/) has already been introduced 
in Eq. ffTTHD (Eq. ffTTHD ). The perturbative expansion of Eqs. f lT37D . f lT39D and f lM is 
explicitly worked out in Appendix |Bl and the ensuing two-loop results are presented and 
discussed in the following Sects. 16.21 and 16.31 

As discussed in Sect. 15.31 and at the beginning of this subsection, the structure of the 
one-loop and two-loop IR factorization formulae f ll04p and flllSp does not specify in a 
unique way the explicit form of the IR operators /^^(e) and I^^{e). Indeed, an IR finite 
redefinition of M.^^^^'^- and Ai^"^^^^' can be compensated by a corresponding redefinition of 
/^j (e) and I^mI^) (^^ two- loop order this redefinition can modify even at 0(l/e)). 

This kind of invariance (and the ensuing arbitrariness) applies to the all-order formulae 
f ll24p and f ll27p . The IR factorization invariance is particularly evident in Eq. f ll27p : this 
equation is invariant under the transformations (redefinitions) |A^^°') C/fin. jA^^"^') and 
Vm{^) — ^ Vm{^) (C^fin.)"^, where f/gn. is an invertible IR finite operator. In the explicit 
expressions (11051) and f lllQp . we have used a 'minimal form' of I^^ and /^^ namely, a form 
in which /^''(e) and I^mI^) include only terms proportional to the IR poles 1/e'^ {k > 1). 
However, the other relations and derivations presented in this subsection are independent of 
this minimal form. In particular, the IR factorization formula fll37p for the splitting matrix 
does not require (or necessarily imply) that the IR divergent operators V(e) and V(e) have 
a minimal form. Incidentally, we note that, in the one- loop expression (p2|) of the two- 

parton splitting matrix Sp^^'^\ the IR divergent operator /^^e) (see Eq. fl93l) ) does not 
have a minimal form, whereas the one- loop multiparton operator /^^^'(e) in Eqs. (1971) and 
has a minimal form. An analogous comment applies to the two-loop results discussed 
in Sects. 16.21 and 16. 3[ 

In the TL colhnear region, strict factorization is valid and, therefore, the result in 
Eq. fll37p takes a simplified form. In Appendix [Cl we show that the IR structure of the 




(139) 




(140) 
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splitting matrix Sp for the multiparton TL collinear limit can be presented as follows: 



(TL coll. lim.) , (141) 



or, in the equivalent iterative form: 



Spipi, ...,p„-P)= iTL(e) Sp{p,, . . . P) + 5p'^"-(pi, . . . P) , (TL coll. lim.) , 



are strictly factorized (completely independent of the non-collinear partons). The explicit 
perturbative expression of iTL(e) up to two-loop order is given in Appendix O 

The all-order IR structure of the TL collinear limit of m = 2 partons was discussed 
in Refs. [3 [8]. Our discussion in Appendix [C] and the expressions in Eqs. (I14ip - fll43l) 
generalize the corresponding results of Refs. [7[ [8] to the case of m > 3 collinear partons. 
Both the operators V and V in Eq. fll37p depend on the non-collinear partons, while 
the operator Vtl in Eq- (11411) is independent of the non-collinear partons. Since Vtl 
eventually originates from V and V, the strictly-factorized form of Vtl implies a non- 
trivial cancellation of the combined dependence of V and V on the non-collinear partons. 
This cancellation constrains the form of V and V and, therefore, it also constrains the 
general colour and kinematical structure of the scattering amplitude operator VAf(e) in 
the IR factorization formula fll27p (we recall that V and V derive from Vm through the 
coUinear-limit procedure in Eqs. fll35p and f ll36p ). This constraint, which is a consistency 
requirement between strict collinear factorization and IR factorization, is particularly sharp 
if the IR divergent operator \m{^) is expressed in its minimal form (see Refs. [TJ |8] and 
Appendix [C|) . 

6.2 Two-parton collinear limit at two-loop order 

The TL collinear limit of two-loop QCD amplitudes was studied in Refs. [271 113] • The 
authors of Ref. [2^ considered the two-parton collinear splitting g g g and, using the 
unitarity sewing method, they performed a direct (process-independent) computation of 
the corresponding two-loop splitting amplitude Split The authors of Ref. [H] exploited 
the universality of collinear factorization to extract Split by taking the collinear limit 
of explicit two-loop results of scattering amplitudes. Considering various scattering ampli- 
tudes with three external QCD partons (and one additional colourless external particle), the 
splitting amplitudes of all the QCD subprocesses A Ai A2 were computed in Ref. [13]. 
In Refs. [271 E], the computation of Split*^^-* was exphcitly carried out up to 0(e°), i.e. by 
neglecting terms that vanish in the limit e — 0. 

In this subsection we examine the two-parton collinear limit in both the TL and SL 
regions. Since we use factorization in colour space, we consider the two-loop splitting matrix 
S'p*-^''. The study of Refs. [271 [13] mostly refers to the splitting amplitude Split which 
controls the collinear behaviour of colour subamplitudes. In the TL collinear region, the 



(142) 



where Sp °' and the IR divergent operator Vtl (or Itl), 




(143) 
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relation between Sp^'^^ and Split*-^^ is exactly the same as at the tree level and at one-loop 
order; we can simply consider Eq. ( !T6|l and perform the replacements Sp^^^ — )■ Sp*-^^ and 
Split'-^^ — )• Split The two-loop validity of the proportionality relation in Eq. ( IT6|) follows 
from the fact that Sp*-^^ only involves a single (and unique) colour matrix for each flavour 
configuration of the splitting process A — )■ A1A2. Indeed, in the case of the subprocesses 
q ^ qg, q qg and g ^ qq, the colour matrix f^^, (see Eqs. (IT^ - (IT^ ) is the sole colour 
structure that is allowed by colour conservation (this conclusion is independent of the 
perturbative order). In the case of the subprocess g — )■ gig2, two different colour matrices, 
namely, faia2a (see Eq. f lT5|) ) and daia2a {date is the fully-symmetrized trace of f^tH^), are 
allowed by colour conservation. However, the corresponding two-loop splitting amplitude 
Sp\it^'^\pi,P2; P) turns out to be antisymmetric with respect to the exchange 1^2 and, 
therefore, Sp^"^^ is necessarily proportional to faja2aSplit^'^\pi,p2',P) [27] (the presence of 



daia2a is cxcludcd, siucc Sp^ {pi,P2', P) is symmetric with respect to the exchange 1 -H- 2). 

The general coUinear limit of two- loop QCD amplitudes A^*^^^ is controlled by the gener- 
alized factorization formula in Eq. f lHOj) (or Eq. f l86|) ). The two-loop splitting matrix Sp^"^^ 
is the new (irreducible) contribution to collinear factorization. Considering the generic 
two-parton collinear splitting A — > , we can write the two-loop renormalized splitting 

matrix in the following general (i.e., valid in both the TL and SL collinear regions) form: 
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' (2) fin. 
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where Sp^^'^^ and Sp*-^'^-* are the tree-level and one-loop renormalized splitting matrices. 



respectively. The one- loop operator 1^ (e) is given in Eq. (p3|l . and the two-loop colour- 
space operator Ic\e) is 
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where the coefficients K and Hf^ are given in Eqs. f ll20p - fll22p . The last term on the 
right-hand side of Eq. (11451) has the following explicit expression: 

,2^\ 2 / \-2e 
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In Eqs. fll44p and (11451) . the terms Sp^^'^\ Sp^^'^^ and 1^ retain their complete de- 
pendence on e. In the limit e — )• 0, the two-loop splitting matrix has IR divergences that 
lead to e-poles of the types l/e"^, 1/e^, 1/e^ and 1/e. The e-poles of Sp^'^'^^ are entirely 

embodied in the first two contributions, /^^ x Sp*-*^'^-* and 1^^ x Sp^^'^\ on the right-hand 

~(2)fin. 

side of Eq. f ll44p . The third contribution, Sp , still depends on e, but it is finite in the 
limit e — )> 0. 
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In the TL collinear region, the expressions in Eqs. ( I144p and fll45p agree with the QCD 
results of Refs. [H [27]. We only note that Refs. [Ill [27] use the 't Hooft-Veltman (HV) 
variant of dimensional regularization, whereas we use the CDR scheme throughout this 
paper. The comparison between these two schemes poses no difficulties, since the results in 
the HV scheme are obtained from those in the CDR scheme by simply replacing the formal 
wave functions u{p),v{p),e{p) of the external collinear partons in Sp with spin polarization 
states of definite (positive and negative) helicity. 

The derivation of the results in Eqs. fll44p -( !T46|) is presented in Appendix [Bl We recall 
that the derivation is based on two input results: the IR factorization formula (I137P (and 
the knowledge in explicit form [H |3] of the two-loop IR factorization formula flllSp of the 
multiparton QCD amplitudes) and the explicit knowledge (see Sect. 14. ip of the one- loop 
splitting matrix Sp^^'^^ to all orders in e in both the TL and SL collinear regions. In 
particular (see Eq. (1240 p and the related discussion), the coefficients of the single-pole 
terms, 1/e, in Eqs. (I144p and (I146P depend on the terms of 0{e^) in Sp^^'^\ 

In the SL collinear region (i.e., su < 0), strict collinear factorization is violated, and 
all the contributions (with the sole exception of Sp^^'^^) to Eq. fll44p depend on the non- 
collinear partons. We recall (see Eqs. (1921) and (195]) ) that, at one- loop order, the factoriza- 
tion breaking terms are due to two-parton colour correlations of the type Tj ■ Tj, where i 
(i = 1 or 2) labels the collinear parton with Zi < and the label j refers to a non-coUinear 
parton. Owing to the iterative structure of Eq. (11440 . Sp^'^'^^ contains factorization break- 
ing terms of the same type as at the one- loop level, and it also contains the 'square' of these 

terms (see the contributions I^J x Sp*-^'^-* in Eq. (11441) and (I^c^)'^ in Eq. (I145p ). Moreover, 
the result in Eq. (I145p shows new aspects of the violation of strict collinear factorization. 
These aspects are clearly illustrated by the main features of the two-loop colour operator 
Ag'(<). 

^ (2') 

The operator (e), which contains double and single poles in 1/e, produces violation 
of strict collinear factorization since it depends on the non-collinear partons. The expression 

_______ - — ' ('2\ 

(11461) shows that is non- vanishing only in the SL collinear region (it requires 2;^ < 0), 

^ (2) 

and that A^ is definitely non-abelian (it is proportional to fabc)- 

The factorization breaking terms of Sp*-^^ are due to the absorptive part of one-loop 
contributions that are present in both abelian and non-abelian gauge theories. The non- 
abelian character of A^?'' originates from a two- loop interference effect (see, e.g., Eqs. ( 12350 
and (12360 ). The absorptive part of the interactions in a loop interferes with the radiative 
and absorptive part^of the interactions in the other loop. Owing to the causality structure 
of the absorptive part, the interferences between the two loops occur at different (asymp- 
totic) times: non-abelian interactions with different (time) orderings do not commute, and 
this produces new factorization breaking terms at the two-loop level. 

^ (2) 

The colour structure of A^ involves three-parton correlations and, specifically, correla- 
tions between a collinear parton and two non-collinear partons (see, e.g.. Fig. [3]-left). Since 
there are two collinear partons, to detect the effect of A^^ we have to consider the collinear 

■l-The absorptive correlations in a loop spoil colour coherence of the interactions in the other loop, so that 
both the radiative and absorptive parts of these interactions contribute to factorization breaking effects. 
As a consequence, has a hermitian and an antihermitian component. 
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limit of amplitudes Ai{pi,p2, ■ ■ ■ ,Pn) with n > 4 QCD partons (if n = 4, Ai necessarily 
involves additional colourless external legs, otherwise the corresponding reduced matrix 
element Ai vanishes). We also note that the two non-collinear partons j and k that are 
colour-correlated by A^"^ must have energies with opposite signs (i.e., Sjk < 0): this pair 
of non-collinear partons consists of a physical initial-state parton and a physical final-state 
parton. In particular, this energy constraint implies that AJ^ does not contribute to the 
two-parton SL coUinear limit of the amplitudes that are involved in lepton-hadron DIS. 

^ (2) 

The operator A^ typically contributes to the SL collinear limit in hadron-hadron hard- 
scattering processes. A simple example with n = 4 QCD partons is the collinear limit of 
the scattering amphtude of the process 'parton + parton — )■ vector boson + 2 partons', 
where one of the final-state partons is collinear to one of the initial-state partons. 




Figure 3: Two-parton SL collinear limit {pi ~ ZiP, i = 1,2) of two-loop QCD amplitudes 
with n = A partons and additional colourless external legs (denoted by the dashed line) 
in parton-parton hard- scattering configurations. Representative colour structure of non- 
abelian factorization breaking correlations according to (left) Eq. Ij^TJ^ o-nd (right) Eg. ( [i^<S| j. 

~ (2) 

We present the explicit expression of A^^ for a generic matrix element with n = A QCD 
partons (see Fig. |3]). We consider only the SL kinematical configurations in which A^^ 7^ 0. 
With no loss of generality, we can limit ourselves to the following kinematical region (the 
other kinematical regions are obtained by re-labeling the parton momenta): 

2:2 < , S34 < , S23 < , 
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and, using Eq. f ll46p . we have 
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Note that the expressions (11471) and (I148P are equivalent in the colhnear hmit. To obtain 
Eq. (I148p . we have implemented the collinear approximation pi ~ ZiP {i = 1,2), and we 
have used the n = 4 equality 



T'a rph rpc 
2 3 4 



f rpa rph rp 

J ahc 1 2 ; 



3 ' 



[n 



This equality derives from charge conservation (i.e., T4 = 
colour algebra, namely, from the general algebraic identity 



4) . 
-(T1 + T2 



(149) 
and from 



iUcT^TUT^ + T't 



(150) 



The factorization breaking correlations depend on the momenta of the non-collinear 
partons. At the one-loop order (see Eq. fl93|) ). this dependence only involves the sign of the 
energy of the non-collinear parton (what matter is simply the physical distinction between 
initial-state and final-state partons). The two-loop operator A^?-* instead depends also on 
the actual size of the momenta of the non-collinear partons. This dependence appears (see 
Eq. (11461) ) through the scale 

We notice that q'^ ~ {(l±<l±fi = ~ Q^) is the component of the momentum P'^ that is 
transverse to the direction of the momenta pj and pk in the reference frame where these 
two momenta are back-to-back. Since in the collinear limit we have Pi ~ ZiP {i = 1,2), 
the scale fllSip can be expressed in terms of the transverse momentum • of one of the 
collinear partons: 

Sjk 

The 'non-collinear' scale in Eq. f ll5ip and the 'collinear' scale si2/{ziZ2) are always 
positive quantities. Thus, we have 

In - .J£_J!L iO = In -JJ-^ ] - in , (153) 

where the logarithm on the right-hand side is a real number (with no imaginary part). 

_____ ^ 

Inserting Eq. f ll53p in Eq. f ll46p . we see that Aj^ has a hermitian and an antihermitian part 
(we recall that at one-loop order the factorization breaking term is purely antihermitian). 
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The hermitian part depends on the non-colhnear scale fllSip . In the antihermitian part, 
the momenta of the non-colhnear partons appear only through the sign of their energies, 
and we explicitly have 



A<?'(.) 



2/ 2^, 
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where the equivalence between the expressions (11541) and f ll55p follows from the conservation 
of the colour charge (see Eq. f ll62p ). We note that the expression f ll55p involves colour 
correlations between a non-collinear parton and the two collinear partons. This type of 
three-parton colour correlations (see, e.g.. Fig. [3]-right) completes the class of factorization 
breaking structures that can appear in the two-loop splitting matrix Sp^'^'^\ 



(2) fin. 



We briefly comment on the finite contribution Sp in Eq. fll44p . In the TL collinear 



' (2) fin. 



at 0(e° 



can be extracted from the direct compar- 

(2) fin. 



region, the explicit expression of Sp 

ison with the results of Refs. [TH ^7\. We have not computed Sp ' in the SL collinear 
region. Its expression contains terms that violate strict collinear factorization. These terms 
produce all types of colour correlations (including three-parton correlations as in Eq. (11550 ) 
that contribute to the IR divergent part of Sp^'^'^\ 



6.3 Multiparton collinear limit at two-loop order 

To our knowledge, no results on the multiparton collinear limit (with m > 3 collinear 
partons) of the two-loop QCD amplitudes A^^^^ are available in the literature. To illustrate 
some main features of the corresponding two-loop splitting matrix Sp^'^^ in the factorization 
formula (Mil) (or (I5B1) ). we proceed analogously to Sect. 15. 3[ Considering the renormalized 
splitting matrix, we introduce the following decomposition in IR divergent and IR finite 
tjGrixis' 

= 5p(2)'i'^- + 5p(2)'^^- . (156) 
All the e-pole contributions to Sp^'^'^^ are included in Sp^'^^'^"', which can be written as 

5p(2)d.. ^ spa,n) + j^Ue) + S^^'^'-- , (157) 

where Sp*-^'^-* is the one- loop splitting matrix, and /^^(^(e) is given in Eq. (pHI) . The new 
two- loop colour operator /^-"^(e) contains e poles, and additional IR divergent terms are 

included in Sp^'^^'^"' . Note, however, that Sp^'^^'^"' is purely non-abelian, it is non- vanishing 
only in the SL collinear regions, and it contains only single poles (i.e., its IR divergences 
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appear at 0{l/e)). All the higher-order poles (starting at order 0{l/e^)) are included in 
the first two terms on the right-hand side of Eq. (11571) . 



The result in Eq. (I157P is derived in Appendix |Bl The derivation is analogous to that 
of the corresponding result (see Eq. (11441) ) for the two-parton collinear limit. In the case 
of the multiparton collinear limit, we do not know the explicit expression of the one- loop 
splitting matrix Sp^^'^^ at 0{e^) in the SL region: this prevent us from computing the 

explicit expression of Sp ^ through the method used in Appendix [Bl 
The explicit expression of the colour operator I^j^c is 



+ 




where 



^(2; 2), 



2tt 



ln(-^^-.0 



-Zi sign s 



iGC j,k€NC 



-Sjk) 



(159) 



Since the flavour-dependent coefficients H^'^'^ (the subscript P in Hp' refers to the flavour of 
the parent collinear parton) are c-numbers, the non-trivial colour-charge structure of I'^^q 

(1) (2- 2) 

is due to IIJq and A^'^ . This structure produces violation of strict factorization in the 
SL collinear regions. In the two-loop splitting matrix Sp^'^'^\ the factorization breaking 
terms start to contribute at 0{l/e^) (we recall that -T^''^ includes the factorization breaking 



(2) 



operator A 



(1) 



of Eq. inm ). 



(2- 2) 

The colour-charge operator A^'^ is responsible for some distinctive features of the 
multiparton collinear limit at two-loop order. This operator is closely analogous to the 
corresponding operator A[?^ (see Eq. fll46p ) for the two-parton collinear limit. By direct 
inspection of Eq. fll59p . we see that the operator A^'^"^ produces IR divergences at the 
level of double poles 1/e^, it is non- vanishing only in SL collinear regions, it is definitely 
non-abelian, and it leads to violation of strict collinear factorization. These factorization 
breaking terms have the form of three-parton correlations: the colour charge of a collinear 
parton is correlated to the colour charges of two non-coUinear partons (the partons j and 
k) that have energies with opposite sign {sjk < 0), and the intensity of the correlation is 
controlled by the transverse- momentum scale -^^ (see Eq. (11511) ). Since 9^p^.^ > 0, 
the logarithm on the right-hand side of Eq. (11591) has a real and an imaginary part (i.e. 
ln(— qr^~ — iO) = ln(qr^~ _^.^//x^) — m) and, correspondingly, the operator A|^'^^ has a 
hermitian and an antihermitian part. 
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Using Eq. f ll59p . the explicit expression of the antihermitian part of ^\y^c 



ieC j,k£NC 

(160) 

^ ^'''^ s s '^^''^'^•^^ • 

(161) 

The expressions f ll60p and f ll6ip are equivalent in the collinear limit (see Eq. f ll62p and the 
accompanying comments in the final part of this subsection). In particular, Eq. fll6ip shows 
that the antihermitian part of A^'^ can be expressed in terms of three-parton correlations 
that involve a non-collinear parton and two collinear partons (the partons i and i) that 
have energies with opposite sign {sa < 0). 

We add some brief comments on the two-loop contributions Sp^'^^'^"' (see Eq. 01570 ) 
and S'p'-^^^'^' (see Eq. fll56p ). The term Sp ^ which is proportional to the single pole 
1/e, is non- vanishing only in SL collinear regions and (analogously to A^'^ ) it involves 
non-abelian factorization breaking correlations between a collinear parton and two non- 
collinear partons (see Eq. (12330 ). The IR finite term Sp^"^^^^' is non- vanishing in both 
the TL and SL collinear regions. In the SL collinear limit, Sp^"^^^^' receives factorization 
breaking contributions from all types of colour correlations that appear in Sp^^^'''^-. In 
particular, Sp^'^^^^' also includes non-abelian correlations between a non-collinear parton 
and two collinear partons (see, e.g., Eq. fll6ip ). 

Factorization breaking terms that correlate the colour matrices of three partons can be 
related by using the following identity: 

fabc Yl E Q^-^^'^) = f-'^ E E ^"^^'^^ ' 

iGC j,k£NC i,eeC j&NC 

j^k i^e 

(162) 

where hi is an arbitrary c-number function that depends on the collinear parton i. The 
identity (11620 relates terms that involve one collinear parton and two non-collinear partons 
to terms that involve two collinear partons and one non-collinear parton. Note that this 
relation requires that the kinematical function hi is independent of the non-collinear partons 
(e.g., of the momenta of the non-collinear partons). If the function hi is simply hi = Q{—Zi), 
we can implement the collinear limit pi ~ ZiP and we can rewrite the right-hand side of 
Eq. (fT62D as follows: 



fabcYl E E E ^"^^'^'^^ig^M®(-^^)®(- 



i,e&cj€Nc i,e€Cj€Nc 



(163) 

The proof of Eq. (11620 (which follows from colour conservation) and the derivation of the 
collinear relation fll63p are presented in Appendix[B](see Eqs. fl2450 - fl249p ). The expression 
(11610 is obtained from (I16O0 by simply using Eqs. (11620 and fll63p . The expression (11550 
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is obtained from fll54p by directly using Eq. fll62p (note that sign(sij) = — sign(s2j) in the 
SL colhnear region, since S12 < 0). 



7 Squared amplitudes and cross sections 

The perturbative QCD computation of cross sections (and related physical observables) 
requires the evaluation of the square of the matrix element A4{pi,p2, ■ ■ ■ ,Pn) and its in- 
tegration over the phase space of the final-state partons. In this section we consider the 
collinear limit of squared amplitudes. In particular, we are interested in the implications of 
violation of strict collinear factorization at the level of squared amplitudes and, possibly, 
of cross sections. 



7.1 The collinear behaviour of squared amplitudes 

We consider the squared matrix element, |A^P, summed over the colours and spins of the 
external QCD partons (see Eq. ([7])): 



|-M'=^''^'-'^^'^^'-(Pi,P2,...) ^ 

{ci} {si} 

Using the notation in colour+spin space (see Eq. ([8])), |A^p can be written as 

\Mipi,P2,...,Pn)\'^ = {M{puP2,...,Pn)\M{puP2,...,Pn)) ■ (165) 



\Mipi,p2, ...)\'^Y.I1 [M''"'''''-iPi,P2, ■ ■ ■ )]^ M'^''^'■■■■''^''^'■■ip^,p2, . . . ) . 

(164) 



The all-order singular behaviour of |A^P, in a generic kinematical configuration of 
m collinear partons with momenta {pi, . . . ,Pm}, is obtained by squaring the generalized 
factorization formula in Eq. fl78|) . We have 



\M\^ c^{M \ Pipi,...,pm;P;Pm+i,...,Pn) \M) , (166) 
where the matrix P is the square of the all-order splitting matrix Sp : 

P = [Sp]^ Sp . (167) 
The loop expansion of the squared splitting matrix P is 

p = p{0,K) + p(l,K) + p(2,K) _^ _ _ _ ^ (^Qg) 

where pC^'^) (with A; = 0, 1, 2, ... ) are the renormalized perturbative contributions. Insert- 
ing Eq. ( I83l) in Eq. fll67p . we obtain the expression of P^'^'-^) in terms of the perturbative 
contributions to Sp : 

pi^'R) = (^Sp^my 5p(o,ii) ^ (169) 
p(i.«) = 5p(^'^) + h.c. , (170) 
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pirn = (^Sp^^y + [(^5p(°'«))^ h.c. 

where the abbreviation 'h.c' means hermitian conjugate. 



(171) 



The perturbative (loop) expansion of the all-order factorization formula (11661) is ob- 
tained by using Eqs. ( I82l) and (11681) . Considering the expansion up to the two-loop level, 
we have 

I A^(0'«) |2 ^ {M^'^""^ I P(°'«) \M^'''') , (172) 



{0,R) 



(173) 



+r(_/^(o,i?)|_^{2,K)^^ ^^^^ 
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M 



+ 



{M 



P(i'«) |A7^^'"^) + c.c. 



(174) 



where the abbreviation 'c.c' means complex conjugate. The tree-level factorization for- 
mula (11721) depends on p(°'^). The one-loop factorization formula (I173p also depends on 
p(i:^)^ whereas P^^'^) enters the two-loop factorization formula (I174p . The expressions 
in Eqs. (I172p . (I173p and (I174p directly correspond to the terms that contribute to the 
order-by-order perturbative calculation of cross sections. 



The kernel P on the right-hand side of the generalized factorization formula (I166P is a 
matrix in colour -|-spin space. The matrix acts on the vector space of the reduced matrix 
element M. {P,Pm+i, ■ ■ ■ ,Pn)- The dependence of P on the spin and colour indices can be 
denoted in the following explicit form: 



(r' r' c' I 

\"-' ) "-'m+l' ■ ■ ■ i^n\ 



P(p 



l5 



) PrTi) -P) Pm+l ) • • • ) Pn) \ C; Cm-\-l ; • • • ; Cn) 



(175) 



where Cm+i, ■ ■ ■ ,Cn are the colour indices of the non-collinear partons (with momenta 
Pm+i, ■ ■ ■ ,Pn) in whereas c and s are the colour and spin indices of the parent collinear 



parton (with momentum P) in (the indices 4 



m+l5 



,c'c' and s' refer to the vector 



space of the complex conjugate matrix element \ ). The matrix structure of Eq. (I175p 
follows from Eq. (I167p and from the matrix structure of the splitting matrix Sp. We briefly 
comment on the dependence of Eq. fll75p on the spin and colour indices, in turn. 

According to the generalized factorization formula (!78|) . the splitting matrix Sp is 
independent of the spin of the non-collinear partons: indeed, Sp only depends on the spin 
indices of the collinear partons and of the parent collinear parton. Since the right-hand side 
of Eq. (11671) (or, equivalently, the left-hand side of Eq. f ll66p ) involves the (implicit) sum 
over the spins of the collinear partons, the squared splitting matrix P can only depend on 
the spin indices, s and s', of the parent collinear parton (as explicitly denoted in Eq. (I175P ). 
The parent collinear parton can be either a quark (antiquark) or a gluon, and we recall 
(see, e.g., Ref. [11]) that the spin structure of [P l^/g is different in these two cases. Fermion 
helicity is conserved by QCD radiation from massless quarks (antiquarks). Therefore, if 
the parent collinear parton is a fermion (quark or antiquark), we can consider the helicity 
basis in spin space and the squared splitting matrix P turns out to be diagonal in this basis: 
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actually, due to parity invariance, we simply have [P]^/^ oc Ss's- An analogous reasoning 
cannot be applied if the parent collinear parton is a gluon. In the gluon case, [P]^/^ has a 
non-trivial dependence on the gluon spin indices s and s' . This dependence has to carefully 
be taken into account (see, e.g., Refs. [35| and [511 [52]) to achieve the cancellation of IR, 
soft and collinear divergences in calculations of cross sections. Moreover (see Refs. [531 El]), 
the dependence of [P]^/^ '-"^ gluon spin indices can lead to physically-observable and 
logarithmically- enhanced ejects in specific kinematical configurations. 

The colours of the collinear partons are (implicitly) summed on the right-hand side of 
Eq. fll67p (or, equivalently, on the left-hand side of Eq. fll66p ) and thus, in general, the 
squared splitting matrix P depends on the colour indices of the parent collinear parton 
and of the non-coUinear partons (as explicitly denoted in Eq. fll75p ). In the TL collinear 
region, strict factorization is valid: since Sp is independent of the colours of the non- 
collinear partons, the squared matrix P is also independent of them and it can only depend 
on the colours, c and c', of the parent collinear parton. Owing to colour conservation, this 
residual colour dependence is, however, trivial: it is diagonal and simply proportional to 
5c'c- In summary, considering the TL collinear region, P is proportional to the unit matrix 
in colour space, and the matrix structure in Eq. (11751) simplifies as follows: 

V{pi,...,Pm]P) , Sc'cSc'c^.^ ... 6c'c„ , (TL coll. lim.) , (176) 



J s's 



where the matrix Vs's only depends on the spin indices of the parent collinear parton. The 
strictly factorized version of Eq. fll66p in the TL collinear region is 

\M\' \v{pi,...,p^;P)] {M\s'){s\M) , (TL coll. lim.) . (177) 

' L J s's 

s,s' 

In the SL collinear region, the simplified structure of Eqs. f ll76p and (I177P is valid only 
at the tree level (i.e. it applies only to p(0''^) in Eqs. fll69p and f ll72p ). Beyond the tree 
level, the SL collinear limit violates strict factorization, and the perturbative contributions 
p(i,-R) g^j^jj p(2,-R) j^qg_ p73p and f ll74p depend on the non-coUinear partons (as generically 
denoted in Eq. fll75p ). In the following subsections, we present explicit expressions for the 
one- loop and two- loop collinear matrices, p(^'^) and p(^'^), and we discuss the structure of 
the terms that produce violation of strict collinear factorization at the squared amplitude 
level in the SL collinear region. 



7.2 Two-parton collinear limit of squared amplitudes 

We consider the two-parton collinear splitting A — > A1A2 in a generic (TL or SL) kine- 
matical region. The corresponding tree-level collinear matrix p(0''^) = p('''^)(pi,p2; -P) is 
strictly factorized. It can be computed (see Eq. ( 1169P ) by squaring the splitting matrices in 
Eqs. ( |T2l) -( lT5|) . The four- dimensional expression for p(°'^) is well known p[5], and the ex- 
plicit d- dimensional expressions in various variants of dimensional regularization are given 
in Ref. [15]. 

The one-loop collinear matrix p(^''^) is computed by inserting the splitting matrix 
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5p(i'-f^) of Eq. ([n2D in Eq. ffTTUD . We straightforwardly obtain the result 



where the IR divergent factor Ip\e) is 



Sp 



Ic (e) 



{0,R) 



Sp 



{1,R) 
H 



+ h.c. 



(178) 



h.c. . 



(179) 



Inserting the explicit expression of the colour operator I(j (see Eq. (!93|) ) in Eq. f ll79p . the 
factor Jp turns out to be a c-number (more precisely, Jp is simply proportional to the 
unit matrix in colour space). We explicitly obtain the following expression: 



r(i), 



e) 



as(/i^) 1 ~ 
27r 2 



Si2 - iO 



+ c.c. 



I ^ ((^12 - Ci - Cs) + ^ (^7i2 - 7i - 72) 

C12 + Ci - C2) /i?(e; ^i) + fCi2 + C2-CA fnie; z^) 



1 

e L 



. (180) 



An important conclusion to be drawn from the result in Eq. f ll78p is that p(^'^) is strictly 
factorized, despite the fact that the corresponding one-loop splitting matrix Sp^^'^^ violates 
strict coUinear factorization in the SL collinear region. The contribution Sp^^'^^ in Eq. ( l92l) 
is strictly factorized, and the factorization breaking terms of Sp^^''^^ are entirely embodied 

in . However, these terms are antihermitian and, therefore, they cancel (see Eq. (11791) ) 
in the computation of the squared splitting matrix P^^'^). The absence of factorization 
breaking terms in p(^'^) implies that, in the specific case of the two-parton collinear limit, 
the factorization structure of Eq. f ll77p is valid up to the one-loop level (this factorization 
structure can be implemented in the one-loop formula fll73p ). 

The result in Eqs. (11780 and (11800 also suggests a practical recipe to compute p(^'^) (in 
both the TL and SL regions) by-passing the violation of strict collinear factorization at the 
amplitude level. The recipe is: consider the expression of Sp^^'^^ in the TL collinear region 
(this expression is strictly factorized, but it depends on the function /(e; x) that is ill-defined 
in the SL collinear region), replace f{e;x) with its (well-defined) real part fii{e;x) (see 
Eq. (IMP ), and use the corresponding expression of Sp^^'^^ to compute p(^'^). This practical 
recipe, which gives the correct result in Eq. (I178p . coincides with the 'effective prescription' 
proposed in Refs. [27j (see Sect. 7.4 therein) to perform the analytic continuation of p(^'^) 
from the TL into the SL collinear regions. An extension of this practical recipe from one- 
loop to two- loop level is not feasible. Indeed, as shown below, the two- loop matrix p(2.^) 
is not strictly factorized in the SL collinear region. 

The one- loop SL prescription (recipe) that we have just illustrated was used in practice 
in the actual NNLO computation of Ref. [55j. This prescription is also consistent with the 
one-loop calculation of the initial-initial three-parton antenna functions (see Eq. (3.19) in 
Ref. [56]). Incidentally, we note that the initial-initial antenna functions derived in Ref. [56] 
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are based on explicit one-loop computations of specific squared amplitudes with n = 3 
QCD partons (e.g., the squared amplitude of the DY subprocess qq — )• 'j*g). Therefore, the 
universality (process independence) of those initial-initial three-parton antenna functions 
is eventually a consequence of the strict factorization of the one-loop coUinear matrix p(^'^) 
for the SL coUinear limit of m = 2 partons. 

The right-hand side of the two-loop coUinear relation f ll74p has four contributions. The 
first three contributions depend on P^'^'^) and P^^'-^) (which are strictly factorized), and the 
last contribution depends on P^^'-^). The two-loop coUinear matrix P^^'-^) can be computed 
by inserting Eqs. ( 192|) . fll44p and f ll45p in Eq. fll7ip . Performing some algebraic operations, 
we obtain 



+ Sp 
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(2) fin. 



where /p^(e) is given in Eq. fllSOp . and the two-loop IR divergent factor /p''(e) is 
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Since /p (e) is a c-number, the first three terms on the right-hand side of Eq. ( 118ip are 
strictly factorized. The last two terms, instead, lead to violation of strict collinear fac- 



(2) fin. 



^-j ..... - (<2\ 

torization. The last term, which depends on Sp , is IR finite. The operator Ap (e) 
is IR divergent and it originates from the colour operator AJ^ on the right-hand side of 
Eq. (11451) : we have 

;i83) 



Ag^(e) 

(2) 



Ag^(e) + h.c. . 



Inserting the explicit expression of A^ (e) (see Eq. (11460 ) in Eq. (11830 . we obtain 
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Note that the argument of the logarithm in Eq. (11840 can be replaced by using the approx- 
imation 

SjpSf,pZiZ2 Sji Sk2 Sj2 Ski 

— ~ — ~ — , (185) 

SjkSl2 SjkSi2 SjkSi2 

which is valid in the collinear limit {pi ~ ZjP, with i = 1, 2). 

The expression (I144p for Sp^"^'^^ includes several terms that violate strict collinear 
factorization, and most of them cancel in the computation of P^^'^). Many factorization 

breaking terms of Sp*-^'^^ have a one-loop origin and are included in the operator /|^\ These 
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one-loop terms do not appear in pC^'^); their cancellation is due to the iterative dependence 
of Sp^'^'^^ on Sp^^'^^ and J^"*, and to the fact that the factorization breaking part of is 
antihermitian. The two-loop factorization breaking operator has instead a hermitian 
component that leads to the operator Ap (see Eqs. f ll83p and f ll84p ) and, thus, to ensuing 
factorization breaking terms in P^^.^). Additional (though IR finite) factorization breaking 

contributions to P'-^'-"-' are produced by Sp 

The IR divergent colour operator Ap (e) is non-abelian and it is non- vanishing only 
in the SL collinear region. The main features of Ap are similar to those of the operator 

- — ' ^2') ' — ' (2\ ' — ' ('2") 

A^ (see the related discussion in Sect. 16.21) . In particular, A^ and, hence, Ap vanish 
in the case of the SL collinear limit of the amplitudes that are involved in lepton-hadron 
DIS. Moreover, it is important to note (as discussed below) that Ap gives a vanishing 
contribution to p(^'^) in the case of scattering amplitudes Ai{pi, . . . ,Pn) with n < 4 QCD 
partons (and an arbitrary number of colourless external legs). Therefore, to detect the 

^ (2) 

effect of Ap at the squared amplitude level, we have to consider the SL collinear limit of 
amplitudes with n > 5 QCD partons. 

^ (2) 

The explicit expression of A^"^ for a generic matrix element with n = 4 QCD partons (if 
n = 3, A^ vanishes trivially, as noticed in Sect. 16. 2p is given in Eq. f ll48p . Independently 

~ (2) 

of the specific kinematical configuration, colour conservation implies that A^ and, hence, 

~ (2) 

Ap are proportional to a single colour charge operator; we have 

(6) oc Uc T| , (n = 4) . (186) 

We note that, in the case of the two-parton collinear limit, the tree- level splitting matrix 
5'p(o>^) fulfils the following colour charge relation: 

(5p(°'«)(pi,P2; P)^ fabcT^n Sp^''^^{p,,p2; P) = (187) 
Therefore, Eqs. f ll86p and f ll87p imply: 

Ag) 5p(°'^) = , (n = 4) , (188) 

and, hence, Ap^ does not contribute to p(^'-^) (see Eq. fll8ip ) if n = 4. 

The relation (11870 follows from colour conservation. The proof is very simple. Using 
Eq. (HT]) and the fact that the Casimir factor C12 is a real c-number {CI2 = C12)) we obtain 

[T^(Ti + T2)2] 5p(°'^) = (T{=Ci2-Ci2Tn = . 

(189) 

Moreover, using elementary colour algebra, the commutator in Eq. (11890 gives 

% [T^(Ti + T2)2] =2fabcT^T^ . (190) 

Inserting Eq. ffT90|) in Eq. ffT89|) . we get the result in Eq. (fTSTD . 

^ (2) 

The operator Ap contributes to the SL collinear limit of the squared amplitudes (with 
n > 5 partons) that are involved in hadron-hadron hard-scattering processes. We explicitly 
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consider a typical example with n = 5 partons: the process 'parton + parton — )■ 3 partons' 
(see Fig. H]) or, more generally, 'parton + parton — )■ X+ 3 partons' (X denotes non-QCD 
particles, e.g. a vector boson), where one of the final-state partons is collinear to one of 
the initial-state partons. We specify the parton momenta as follows: the two initial-state 
partons have momenta —pi and —p^ (i.e., the 'energies' p^ and are negative), and the 
three final-state partons have momenta P2,P'i and p^ (i.e., P2,P^ and p^ are positive). Since 
Z2 < < Zi and S45 > 0, there are only two different colour operators that contribute to 
the expression fll84p of Ap'*: 
colour conservation (i.e., T5 = 



these operators are fabc T2T^T^ 



and fabc T^T^T§. Using 



-(Ti + T2 + T3 + T4)) and the identity f lTSOD . we obtain 



r rpa rpb rjic 

Jabc -'2 -'S -'S 
(2) 



/rpa rpb rpc 
abc J 2 -'a -'4 



/rpa rpb rpc 
afec J 2 3 1 



{n = 5) 



:i9i] 



and we can write Ap in terms of the colour operators fabc Tj^T^TI and fabc Tf T2 T|. 
From Eq. f ll84p . we thus obtain 

-2e 



Ag^(e) 



TT 



-Sl2 



fabcT^T'Tl In 



■S34 ggp 
■535 ■S4P 



1 



2e' 



e V ^1 



+ fabcT^T^T^ X (■■■) . 



(192) 



On the right-hand side, we have not written the explicit expression for the coefficient 
of the colour operator fabcTiT2T^: indeed, as a consequence of the relation f ll87p . this 
colour operator does not contribute to the collinear matrix P^^.^) (gee Eq. f llSip ). The term 
proportional to fabc T| in Eq. (fT92|) (see Fig. Hj) produces a non- vanishing factorization 
breaking contribution, (5p(°'^))t A.f Sp^°'^\ to P^^.^). 




Figure 4: Squared amplitude with n = 5 QCD partons in parton-parton hard scattering 
(the dashed line cuts the final-state partons). Representative colour structure of non-abelian 
factorization breaking correlations that accompany the two-parton SL collinear limit (pi ~ 
2;jP, i = 1,2) at two-loop order. 

We note that, in the case of a generic scattering amplitude with n QCD partons, the 
colour structure of the contribution of Ap^ to P^^'-^) (see Eq. fllSip ) can be simplified. We 
have 

Ag)(e) = A^^e) p(°'«) , (193) 
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(2) 

where the colour operator Ap acts on the colour space of the n — 1 external QCD partons 
of the reduced matrix element Ai (we recall that p(0'^) is strictly factorized and, hence, 
it is proportional to the unit matrix in colour space). The simplified colour structure in 
Eq. (]193p follows from the fact that we can explicitly evaluate the action of the colour 
charge Tj {i = 1,2) of the coUinear partons on the tree- level two-parton splitting matrix 
^p(o,-R)_ Using basic colour algebra relations, we straightforwardly find (z = 1,2) 

(5p(°'^)(pi,P2;P))^ Tt Sp^''^\p,,p,;P) = TJ^^ P(0'^) , (194) 

where the colour operator T^^^ p is proportional to the colour charge of the parent coUinear 

parton A with momentum P. The explicit expression of T^^^-^p depends on the flavour of 
the coUinear parton i and of the parent collinear parton. Considering the flavour structure 
of the various splitting processes A — )■ Ai A2, we have 

^'^-f. T,,,f = ^T,, (195) 



f f 9 (/ = g or g) : 


Tu)f = 




9^99- 


'^(9) a = 




g^qq : 


^(9)9 = 





2Cf ^ ' 2Cp 



(196) 



T, + d,) , T(,-), = i (T,-d,) ,(197) 

where Tf and Tg are the customary colour charges of a fermion and a gluon, while dg is 
the gluon colour matrix in a symmetric 'octet' state and its matrix elements are defined as 

id\^= d,ac , datc = 2TT{t''{t\f}) , (198) 

and dabc is the fully-symmetrized trace of f^t^f^. The explicit expression of the colour 

(2) 

operator Ap in Eq. (I193P is 



A?\e) = {'^^\ i^Y'^^fa.cY. E ^(:)pT,''T,^e(-..)sign(.,,)e(-.,,j 



(S^pSj^pZiZ2 

X In I — 

Sjk Sl2 



2e2 e 



i=l,2 j,k=3 

' In ^ 



(199) 



This expression is directly obtained by inserting Eqs. fll84p and fll94p in the left-hand side 
of Eq. ffT93|) . 



7.3 Multiparton collinear limit of squared amplitudes 



In this subsection we consider the multiparton splitting A — )■ Ai . . . Am of m {m > 3) 
collinear partons. The corresponding tree- level collinear matrix p('''^) fulfils strict fac- 
torization. The explicit expressions of p(°'^) for all the flavour configurations of m = 3 
collinear partons were computed in Refs. |T8l |19l [11] . 

The one-loop collinear matrix P^^'-^) is obtained by inserting Eqs. (196 p and fl97p in 
Eq. ([170]) . We have 



Sp 



(0,R) 



+ 



Sp 



Sp 



(l)fin. 



+ h.c. 



(200) 
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where /^p is given in terms of the colour operator 1^^^^ in Eqs. (p7I) and (p9l) : 

I^U^) = itU^) + h.c. . (201) 
Using Eq. (!99|) . the exphcit expression of /^p is 



^7^ 



(202) 



The factorization breaking contribution A^'^^ (see Eq. f llOOp ) to I^^c antihermitian 
and, thus, it cancels in Eq. fl20ip . The operator /^p depends on the colour charges of the 
coUinear partons, but it is independent of the non-coUinear partons. Therefore, the first 
term on the right-hand side of Eq. f l200p does not violate strict coUinear factorization. This 
term produces a strictly factorized and IR divergent contribution to p(^'^). 

The second term on the right-hand side of Eq. f l200p is IR finite, since it depends on the 
IR finite contribution Sp^^^^'^' to the one-loop splitting matrix. As mentioned in Sect. 15. 3^ 
^p(i)fin. pQj^i^g^jj^g factorization breaking terms in the SL coUinear region and, therefore, 
their interference with Sp^^'^^ can produce violation of strict collinear factorization at 
the squared amplitude level. In this respect, the main difference between the two-parton 
and multiparton collinear limits is that in the latter case several different colour structures 
contribute to the splitting matrix Sp already at the tree level (in the two-parton case, Sp^^^ 
involves a single colour structure: see Eq. ( 1TB]) ). The result in Eqs. f l200p and fl202p for the 
multiparton collinear limit shows that the IR divergent part of p(^'^) is strictly factorized, 
while the IR finite part can contain terms that produce violation of strict factorization in 
the squared amplitudes at the one-loop level (see Eq. f ll73p ). 

The two-loop collinear matrix P^^-^) (^^ee Eq. fll7ip ) of the multiparton collinear limit 
can be computed by using Eqs. ( 1^ . (|57p . fll56p . fll57p and the explicit form of the colour op- 
erators I'^^Q, and I^^Q, (see Eqs. f M?]) and f ll58p ). Performing some straightforward (though, 
slightly cumbersome) algebraic operations, we can write the result in the following form: 

p(2.R)^p(2,iJ)^p(2,iJ) ^ ^203) 

where the contribution P^^'^-* is strictly factorized, while the term Pfff^"* includes all the 
contributions that violate strict collinear factorization (this term also includes additional 
contributions that are strictly factorized). 

The strictly-factorized term P\ ' does not depend on the non-collinear partons. Its 
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explicit expression is 



,(2,ii) 



+ 



where -T^p is the one-loop operator in Eq. (12021) . and the colour operator /, 



27r \ e 

2^\ 2 



K J«p(2e 
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(2) 




(0,R) 



(1) 

mC,f. 



from /^^(^ (see Eq. (!99|) ) by removing its factorization breaking part A, 



(1) . 

mC- 



(204) 
is obtained 

(205) 



The contribution P^^f^"* on the right-hand side of Eq. (12031) has the following expression: 



>(2,i?) 
n.f. 



Sp 



{0,R) 



+ 



Sp 



{0,R) 



J2;2) 



r(l) 



e)Sp 



{0,R) 



:^p{e)Sp^''>^^^- + Sp'^'- + Sp 



-(2) div. 



(2) fin. 



+ {Sp 



(l)fin. 



Sp 



(l)fin. 



where the colour operator A^p'' is 



(2; 2) 



'e) + h.c. 



(1) 



+ h.c. 



(206) 



(207) 



The last term on the right-hand side of Eq. fl206p is IR finite, since it is given by the 
square of the IR finite contribution Sp''^^^^' to the one- loop splitting matrix (see Eq. ( |96|) ). 
As already recalled, Sp^^^^^' contains factorization breaking terms in the SL collinear re- 



gion. 



The term in the square bracket on the right-hand side of Eq. 



is IR divergent. This 

term depends on Sp''^^ on the colour operator /^p in Eq. (I202p . and on the contributions 
5p(2)fin- (^ggg ]T;q_ (I156P ) and Sp^ (see Eq. (11570 ) to the two-loop splitting matrix 
5'p(2,R)^ As stated in Sect. 16.31 Sp^'^^^^^' contains factorization breaking contributions that 
are IR divergent at the level of single poles 1/e. The operator /^p contains also a double- 
pole term with a c- number coefficient (see Eq. f l202p ) and, thus, it produces the following 
IR divergent contribution to the square-bracket term in Eq. (I206p : 



+ 0{l/e) 



(208) 



We note that the double-pole term in Eq. (I208p is simply proportional to the factorization 
breaking contribution to the one- loop collinear matrix p(i'^) (see Eq. (I200p ). 
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The colour operator A^'p in Eq. (120 7p depends on the two-loop factorization break- 



ing term A^',^^ (see Eqs. (I158P and (11591) ). on the one- loop operator J^p(e), and on the 
one- loop factorization breaking term A^''^ (see Eqs. (j99l) and (llOOp ). The commutator on 
the right-hand side of Eq. (120 7p originates from a non-abelian interference of two types of 
one-loop contributions: the one-loop absorptive (antihermitian) contribution to the scat- 
tering amplitude (the absorptive interaction involves a coUinear and a non-collinear parton) 
and the one- loop radiative (hermitian) contribution to the complex conjugate amplitude 
(the radiative interaction involves two coUinear partons). Using Eqs. (llOOp and (I202p . we 

(2- 2) 

evaluate the commutator and we obtain the following explicit expression of A^'p : 



(1) 



^(2; 2), 



X 



27r 



abc 



ieC j,keNC 



T^T^n <S>i~z^) sign(s,,) Qi-s,k) In 



- E E 



T^^T^T^'Qi-Zi)Qi-s^e)signis,,) In 



i/eC j€NC 



^jP ^kP 
Sjk A'-o 



2 



. (209) 



Note that the argument of the logarithms in the square bracket depends on the scale /^o, 
which is arbitrary (we can possibly set fio = fi). However, the expression (I209p is actually 
independent of /io; the independence of /io directly follows from the relation (I162p (to apply 
Eq. (11620 . we note that the explicit constraint su < can be removed from the second 
term in the square bracket of Eq. (12090 : indeed, the sum of the terms with su > gives a 
vanishing contribution to Eq. (12090 ). 

The operator A^p^ is IR divergent and proportional to the double pole 1/e^. The first 
term in the square bracket of Eq. (I209p involves non-abelian contributions that correlate 
the colour charges of two non-collinear paitons (and a coUinear parton). In Eq. (12060 . these 
contributions produce factorization breaking terms that cannot be cancelled by the IR di- 
vergent term of Eq. (I208p (indeed, as mentioned in Sect. 15. 3[ Sp^^^^^' embodies factorization 
breaking correlations with a single non-collinear parton). Therefore, the explicit expres- 
sion (I209p of A^p^ shows that the two-loop multiparton coUinear matrix P^^'^) necessarily 
includes non-abelian contributions that lead to violation of strict coUinear factorization at 
the squared amplitude level. 

A special exception to this conclusion about violation of strict coUinear factorization 
regards the case of the SL coUinear limit in lepton-hadron DIS. In the DIS case, the operator 

[2' 2") 

A^'p effectively takes a form that is independent of the non-collinear partons. To show 
this, we consider a DIS matrix element A^(pi, . . . ,Pm,Pm+i, ■ ■ ■ ,Pn)- the initial-state parton 
has 'outgoing' momentum —pi {p^ < 0), the momenta of the final-state coUinear partons 
are P2, ■ ■ ■ ,Pm, and the momenta of the final-state non-coUinear partons are Pm+i, ■ ■ ■ ,Pn- 
If j, k G NC, we have Sjk > 0: thus, the correlation terms of Eq. (I209p that depend on two 
non-collinear partons vanish. As for the remaining correlation terms of Eq. (I209p . the sign 
of Sij is independent of j if j G NC (actually, Sij > if i G C and Zi < 0); therefore, we 
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can perform the sum over j by using colour conservation 
final expression of ^^'p lepton-hadron DIS is 



'jeNc'^j ~ 



ErecTr). The 




X J2 T^nY.T^Q{-z,)Q{-su) In 



(-^) ' ^^^S)- (210) 



This expression is independent of the non-coUinear partons and, thus, it effectively has a 
strictly-factorized form. This form is a consequence of a colour coherence mechanism (due 
to DIS kinematics and colour conservation). 

In generic kinematical configurations (typically, those that occur in hadron-hadron 
hard-scattering processes), the two-loop multiparton collinear matrix P^^'-^) is not strictly 
factorized. We note that the two-loop collinear formula f ll74p for the squared amplitudes 
includes factorization breaking contributions that are due to both P^^'-^) and P^^'-^). The 
factorization breaking contributions that are due to P^^'-^) involve correlations with one 
non-coUinear parton, while the contributions due to p(^'^) also involve correlations with 
two non-coUinear partons. 

7.4 Cross sections and violation of strict collinear factorization: 
some remarks 

The applicability of perturbative QCD to the calculation of cross sections for hard-scattering 
processes is based on the universal (process-independent) factorization theorem of mass sin- 
gularities [29]. According to this factorization picture, the sole uncancelled IR divergences 
that are eventually encountered in the computation of inclusive partonic cross sections are 
due to partonic states whose momenta are collinear to the momenta of the colliding par- 
tons or of triggered final-state partons; these uncancelled divergences are factorizable in 
a process-independent form and, therefore, they can be removed by a formal redefinition 
('renormalization') of the 'bare' parton densities and parton fragmentation functions. 

The violation of strict collinear factorization at the level of squared amplitudes certainly 
challenges the validity of universal mass- singularity factorization at the cross section level. 
We present some comments and remarks on this issue. 

We first present a general comment. The singular behaviour of the two-parton collinear 
limit at one-loop order is one of the key ingredients that are used to handle IR divergences 
and mass singularities and to cancel the IR divergences in perturbative QCD computations 
of hard-scattering processes at the NNLO (see, e.g., Refs. [SZlESlESlEDlEIlESlEilE^and 
references therein). In this context, it is reassuring that the one-loop two-parton collinear 
matrix p(^'^) is strictly factorized even in the SL collinear region (see Eqs. (I178l) - fll80p ). 
This result guarantees that the extension of NNLO methods from lepton-lepton collisions 
to lepton-hadron and hadron-hadron collisions does not involve additional conceptual dif- 
ficulties related to the violation of strict collinear factorization. 
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In the following, rather than considering the issue of mass-singularity factorization in 
completely general terms, we limit our discussion to the simplest case in which our study 
has definitely uncovered the presence of strict-factorization breaking effects at the squared 
amplitude level. We thus consider the inclusive production of a high-p^ hadron or jet (with 
at least one recoiling jet) in the collision between two high-energy hadrons. 

At the leading order (LO) in QCD perturbation theory, this production process is con- 
trolled by the square of the tree-level amplitude (with n = 4 partons) of the corresponding 
partonic subprocess 

parton + parton — > 2 partons . (211) 

Part of the higher-order QCD corrections to this '2 — 2' partonic subprocess are obtained 
by considering the squared amplitude of the process 

parton + parton —> 3 partons , (one low— final — state parton) , (212) 

in the kinematical region where one of the three final-state partons (the 'low-p^ parton') 
is coUinear to one of the initial-state partons. The IR divergences produced by the phase 
space integration over the SL coUinear region of the low-pT final-state parton have to be 
factorized with respect to the corresponding LO '2 — i- 2' partonic subprocess, and the IR 
factor has to be strictly factorized, namely, it has to be independent of the kinematical and 
colour flow structures of the LO partonic subprocess in Eq. f l21ip . As explicitly shown in 
Sect. \7.'2\ the squared amplitude of the process in Eq. f l212p violates strict factorization in 
the SL coUinear region. This violation of strict factorization starts at the two-loop level (see 
Eqs. fllSip and f ll92p . and Fig.H]) and, therefore, it may invalidate the factorization theorem 
of mass singularities starting at the N^LO in QCD perturbation theory. The validity of 
the factorization theorem can be recovered only if the N^LO factorization breaking effects 
produced by the two-parton SL coUinear limit of the process in Eq. fl212p are cancelled by 
corresponding factorization breaking effects produced by other partonic subprocesses. An 
explicit quantitative proof of such a cancellation would represent a highly non-trivial check 
of the validity of the factorization theorem beyond the NNLO. 

We continue our discussion of the cancellation mechanism at the qualitative level and, 
eventually, we shall identify a sole additional partonic configuration that can lead to fac- 
torization breaking correlations with the same structure as in Eq. fll92p (see Fig. H]). The 
contribution of this partonic configuration (after combination with all the other IR diver- 
gent terms) can cancel the violation of strict coUinear factorization that is produced by the 
subprocess in Eq. f l212p . 

To reach this conclusion, we consider the various high-order subprocesses that lead to 
the N^LO IR divergent contributions to the cross section controlled by the LO process in 
Eq. fl21ip . and we group them in three classes: i) virtual and soft-parton subprocesses, ii) 
collinear-parton subprocesses, and iii) mixed soft-coUinear parton subprocesses. 

The class i) contains the virtual three-loop corrections to the '2 — > 2' process in 
Eq. (12 lip and the (3 — A;)-loop corrections to the processes 'parton + parton — )■ 2 par- 
tons -|- k soft partons' with k = 1,2, 3. Owing to kinematics, the IR divergences produced 
by this class cannot cancel the factorization breaking effects due to the two-loop corrections 
of the process in Eq. fl212p . In fact, the virtual and soft-parton radiation that accompanies 
the basic '2 — )■ 2' partonic process cannot match the kinematics of the process 'parton 
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+ parton — )■ 3 partons' in the region where the final-state colhnear parton has a large 
longitudinal-momentum fraction (e.g., p2 — Z2P with \z2\ ~ |-2i| ~ in Eqs. fllSip and 
(11921) . and in Fig. H]). The subprocesses of the class i) can (partly) cancel the IR diver- 
gences of the subprocess in Eq. (12121) only in the kinematical region where the final-state 
collinear parton is also soft (e.g., at the kinematical endpoint 2:1 = 1 — 1 in Eq- (I18ip ). 

The class ii) contains the subprocess in Eq. fl212p at the two-loop level and two other 
contributions: the tree-level subprocess 'parton + parton — )■ 5 partons' with three low-p^ 
final-state partons, and the one-loop subprocess 'parton + parton — )■ 4 partons' with two 
low-pT final-state partons. The subprocess 'parton + parton — )■ 5 partons' produces IR 
divergent terms from the tree-level SL collinear limit of m = 4 partons: this tree-level 
collinear limit is strictly factorized. The subprocess 

parton + parton — )■ 4 partons (two low— final— state partons) (213) 

produces IR divergent terms from the one-loop SL collinear limit of m = 3 partons: this 
collinear limit violates strict factorization but, as discussed in Sects. 15. Ij 15.31 and [7.31 (see 
Eq. fl200p ). the corresponding factorization breaking terms involve correlations with a single 
non-collinear parton. Since the factorization breaking effects due to the two-loop subprocess 
in Eq. (12120 also involve correlations with two non-collinear partons (see Eq. (I192p and 
Fig.Hj), these effects cannot be fully cancelled by the contributions of the other subprocesses 
in the class ii). 

The class Hi) contains two tree-level subprocesses: the subprocess 'parton + parton 
— )■ 4 partons + 1 soft parton', which produces IR divergent terms from the SL collinear 
limit of m = 3 partons (two of the non-soft final-state partons have low Pt), and the 
subprocess 'parton + parton — i- 3 partons + 2 soft partons', which contributes through the 
SL collinear limit of m = 2 partons (one of the non-soft final-state partons has low p^)- 
Owing to colour coherence (see Sects. 3.4 and 3.5 in Ref. [11]), the singular collinear factors 
are strictly factorized in these mixed soft-coUinear limits of tree-level QCD amplitudes. The 
remaining contribution to the class Hi) is due to the subprocess 

parton -|- parton — )■ 3 partons -|- 1 soft parton , (one low— Pt final — state parton) , 

(214) 

in the kinematical region where one of the three non-soft final-state partons is collinear to 
one of the initial-state partons (Fig. [5]). The subprocess in Eq. (12140 contributes to one-loop 
order. The structure of the absorptive part of the one-loop interaction limits the coherent 
action of the two collinear partons in the SL region (the mechanism is analogous to that 
discussed in Sect. 14.50 . As a consequence, the one-loop mixed soft-collinear limit of the 
process in Eq. fl214p leads to a singular SL collinear factor that violates strict factorization. 
Moreover, this singular factor can also produce correlations with two non-collinear partons 
(see Fig. [5]), whose structure is analogous to that in Eqs. (11810 and (11920 (see Fig. H]). 

In summary, the two-loop radiative corrections to the subprocess in Eq. (12120 and the 
one-loop radiative corrections to the subprocesses in Eqs. fl213p and (12140 produce violation 
of strict factorization in the SL collinear region. The factorization breaking effects are due 
to the singular partonic configurations with m = 2 collinear partons, m = 3 collinear 
partons, and m = 2 collinear partons plus one radiated soft parton. The factorization 
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Figure 5: Squared amplitude for parton-parton hard- scattering with three QCD partons and 
one soft gluon in the final state (the dashed line cuts the final-state partons). Representative 
colour structure of non-abelian factorization breaking correlations that accompany the mixed 
soft-collinear limit (with two collinear partons, pi ~ ZiP, i = 1,2, in the SL region) at one- 
loop order. 

breaking terms involve correlations with one and two (in the case of the subprocesses in 
Eqs. (12121) and (I214p ) hard non-collinear partons. 

In the inclusive hadropro duct ion of a single high-pT jet (or hadron), this violation of 
strict collinear factorization leads to N'^LO contributions that are separately IR divergent; 
the total contribution of these IR divergent terms should vanish to guarantee the validity 
of the factorization theorem of mass singularities. 

Even if the IR cancellation occurs, the violation of strict collinear factorization leads 
to residual IR finite terms. These terms can produce observable contributions that are 
logarithmically enhanced (and, therefore, large) in particular kinematical configurations 
where virtual (see Eq. (I212p ) and real (see Eqs. fl213p and (I214p ) radiative corrections are 
highly unbalanced. For instance, in these kinematical configurations the cancellation of the 
IR pole term of (9(l/e^) in Eq. f ll92p can lead to a residual double-logarithmic contribution. 
The logarithmically-enhanced terms due to the violation of strict collinear factorization 
have a distinctive signature: the factorization breaking correlations with the non-collinear 
partons (see, e.g., Eq. f ll92p and Figs. H] and [5]) produce 'entangled logarithms', namely, 
logarithmic terms whose coefficients get tangled up with the colour flow and kinematical 
structure of the lowest-order hard-scattering subprocess (e.g., Eq. fl21ip ). In the following, 
we mention some specific examples of processes that can exhibit entangled logarithms at 
various perturbative orders. 

The partonic subprocess in Eq. f l21ip also controls the LO inclusive hadropro duction of 
a pair of nearly back-to-back high-p^ hadrons or jets. In the kinematical region where the 
total transverse momentum of the pair is small (e.g., much smaller than the invariant 
mass of the pair), the perturbative QCD calculation leads to large contributions that are 
proportional to powers of IuQt- Owing to the violation of strict collinear factorization in 
the subprocesses of Eqs. (I212p - fl214p . part of these logarithmic contributions can be due 
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to entangled logarithms (the entangled logarithms are absent if the small-Qy triggered 
final-state system includes no QCD partons [6^^ 'SS]) that first appear at the N^LO. The 
presence of these small-Qr entangled logarithms is also related to issues that arise in the 
context of factorization of transverse- momentum dependent distributions [301 EH ES] . 

The perturbative QCD calculation of the inclusive hadropro duct ion of three high-pj- jets 
is also affected by large logarithms, Inp™™, in the kinematical region where the transverse 
momentum p™™ of the lowest-p^ jet is small (e.g., much smaller than the transverse mo- 
menta of the other two jets). In this three-jet production process, the partonic subprocess 
in Eq. fl212p enters at the LO, and the partonic subprocesses in Eqs. fl213p and f l214p first 
contribute at the NLO. The two-loop (in the case of the subprocess in Eq. f l212p ) and one- 
loop (in the case of the subprocesses in Eqs. f l213p and f l214p ) violation of strict coUinear 
factorization can produce small-p™" entangled logarithms starting from the NNLO pertur- 
bative calculation. 

Another example of entangled logarithms is represented by the super- leading ('non- 
global') logarithmic terms discovered [33] in the calculation of the cross section for the 
hadroproduction of a pair of jets with a rapidity gap between them. These super-leading 
logarithms occur at the N^'LO in QCD perturbation theory. 



8 Summary 

We have studied the singular behaviour of QCD scattering amplitudes in the kinematical 
configurations where the momenta of two or more external partons become collinear. We 
have shown that, beyond the tree-level, strict (process-independent) collinear factorization 
is violated in the SL collinear region. We have introduced a generalized form of collinear 
factorization (see Eq. f l33|) in Sect. 14. ![ and Eqs. ( ITHj) and f l79|) in Sect. 15. 2p of the all-order 
scattering amplitudes, which is valid in both the TL and SL collinear regions. In the SL 
region, the singular collinear factor retains some process dependence, since it depends on 
the momenta and colour charges of the non-collinear partons. Owing to this colour depen- 
dence, the formulation of collinear factorization directly in colour space, through the use of 
collinear splitting matrices, is particularly suitable. An equivalent formulation in terms of 
colour subamplitudes and collinear splitting amplitudes is feasible (see Appendix |A]) . In the 
TL collinear region, strict factorization is recovered because of colour coherence. In the SL 
collinear region, colour coherence is limited by the causality structure of long-range gauge 
interactions (roughly, the distinction between initial-state and final-state interactions), and 
this produces absorptive contributions that eventually originate in the strict-factorization 
breaking phenomenon (see Sect. 14.5]) . 

In the case of the two-parton SL collinear limit, we have computed the ((i-dimensional) 
one- loop splitting matrix (Eqs. f l34p and f l45p ) to all orders in the dimensional regularization 
parameter e. The SL result explicitly shows that it cannot be obtained from the previously 
known TL result by using crossing symmetry. The strict-factorization breaking terms are 
antihermitian ('imaginary') at one-loop level. The expressions in Eqs. f l52]) and f l67]) . which 
refer to the SL collinear limit in hadron-hadron collisions, are the simplest examples of 
terms that violate strict collinear factorization. At two-loop level (Sect. 16. 2p . we have 
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explicitly computed all the IR divergent contributions (i.e., the e poles) to the splitting 
matrix. The two-loop result shows the presence of both hermitian ('real') and antihermitian 
terms that violate strict collinear factorization. 

We have derived the structure of the IR divergences of the multiparton collinear limit 
to all orders in the loop expansion (Sect. 16. ip . We have explicitly computed one-loop 
(Sect. 15. 3p and two-loop (Sect. 16. 3p IR divergent contributions to the SL collinear limit of 
m > 3 partons, thus extending the results of the two-parton collinear limit. 

The SL collinear limit of the scattering amplitudes that are involved in lepton-hadron 
DIS (namely, parton radiation collinear to the sole initial-state parton) is a special case, 
since all the non-collinear partons are produced in the final state. The one-loop expression 
of the two-parton SL splitting matrix 'effectively' takes a strictly-factorized form, in which 
there is no explicit dependence on the non-collinear partons (see Sect. 14. 4p . This effective 
strict factorization of the SL collinear limit can be a general feature of the DIS kinematics 
(there are no initial-state interactions between collinear and non-collinear partons). The 
one-loop and two-loop terms of the multiparton splitting matrix that we have explicitly 
computed fulfill this effective strict factorization in the DIS kinematical configuration. 

In hadron-hadron collision configurations, the two-loop SL splitting matrix has fac- 
torization breaking terms that are definitely non-abelian (see Sects. 16.21 and 16. 3p . These 
terms involve correlations between the momenta and colour charges of three partons and, 
in particular, between a collinear parton and two non-collinear partons. These non-abelian 
factorization breaking effects appear in the SL collinear limit of two-loop scattering ampli- 
tudes with at least five external legs and, in particular, n > 4 external QCD partons. 

Owing to their absorptive origin, strict-factorization breaking effects partly cancel at 
the level of squared amplitudes. Nonetheless, strict factorization is certainly violated in 
the SL collinear limit of two-loop squared amplitudes for parton-parton hard scattering 
with at least three final-state partons (see Sect. 17.21 and, in particular, Eq. (I192p ). These 
factorization breaking effects have consequences in perturbative QCD calculations of hard- 
scattering cross sections in hadron-hadron collisions (Sect. I7.4p . The violation of strict 
factorization affects the non-abelian structure of logarithmically-enhanced terms at NNLO 
and higher orders; it has implications on various factorization issues and, in particular, 
it challenges the validity of mass-singularity factorization in jet and hadron production, 
starting from the N'^LO. 

A Appendix: Multigluon amplitudes and the struc- 
ture of the coUinear hmit 

In this Appendix we consider the specific case in which the matrix element Ai{pi,p2, ■ ■ ■ , Pn) 
is a multiparton scattering amplitude with n external gluons. We recall the decomposition 
of the matrix element in colour subamplitudes, and we illustrate the behaviour of the colour 
subamplitudes in the collinear limit. 

The colour indices of the n external gluons of A^(pi, . . . ,p„) are denoted by ai, . . . , a^. 
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At the tree level the pure multigluon amplitude Ai^^^ (see Eq. ([T])) can be expressed as 

follows [m 

. . = (72)" Tr . .T") ^(°)(l,...,n)+ non-cyclic perms. , (215) 

where t'^^'s are colour matrices in the fundamental representation (see Eq. ( [TTj) ). and the 
sum of terms on the right-hand side extends over the (n — 1)! non-cyclic permutations of 
the set {1,2, . . . ,n} of the external legs. This is the customary colour decomposition in 
colour subamplitudes. Each tree-level colour subamplitude, ^'■^''(l, . . . ,n), is independent 
of the colour indices, and it embodies the kinematical dependence on the momenta and 
spin polarizations of the external gluons. Note, however, that the colour subamplitudes 
are colour ordered, namely, the functional form of A^^\l, . . . ,i, . . . , j, . . . ,n) depends on 
the specific ordering (1, . . . , . . . , j, . . . , n) of the external legs in the argument of A^^^ 
(e.g., A^^\l, . . . ,i, . . . , j, . . . ,n) 7^ A^^\l, . . . , j, . . . ,i, . . . ,n)). The overall normalization 
of the colour subamplitudes in Eq. fl215p is adjusted to that used in Refs. [121 ES]; the 
only difference between A^^\l, . . . , n) and the corresponding subamplitude A^^'^^{1, . . . , n) 
in Refs. ^T2\ [25] is due to the overall factor g^''^, which we have included in the definition 
of ^(°)(l,...,n). 

At the one- loop level the colour structure of the multigluon amplitude Ai^^^ (see Eq. ([1])) 

is [Sa [H 125] 

M^^^ {pi,...,Pn) = 1(^2)" Tr(^'^^..^'^") ^(^)(1, ... ,n) + non-cyclic perms. | 

+ double trace terms . (216) 

The term in the curly bracket has the same colour structure of Eq. f l215p . The subamplitude 
^'^^''(1, . . . ,n), which is called leading-colour subamplitude (or primitive amplitude), is the 
one-loop analogue of the tree-level colour subamplitude A^'^\l, . . . ,n). The remaining 
terms on the right-hand side of Eq. ( I216p are proportional to colour factors that involve 
the product of two traces of matrices, namely, Ti {t^'^ . . .t^'') Ti {t"-'=+^ . . .t^"). Each 
double trace is multiplied by a corresponding kinematical factor called subleading-colour 
subamplitude (or subleading-colour 'partial amplitude'). The subleading-colour partial 
amplitudes are in fact not independent of the leading-colour subamplitudes ^^^-"(l, . . . ,n); 
rather, they can be expressed as sum over permutations of the arguments of the latter. 
Owing to this linear dependence, it suffices to examine the coUinear limit of the leading- 
colour subamplitudes. In particular, we do not explicitly consider subleading-colour partial 
amplitudes and their contribution to Eq. f l216p . 

As in the case of the tree- level colour decomposition in Eq. f l215p . our normalization 
of the one-loop subamplitude A^^\l, . . . ,n) is adjusted to that used in Refs. [121 125]. To 
be precise, at the one-loop level the authors of Refs. [121 125] decompose the QCD leading- 
colour subamplitude in two terms, y4n~'°°'''^' and ^^~^°°p[^/^]: the subscripts [1] and [1/2] 
respectively refer to the contribution of a gluon and a quark circulating in the loop. The 
relation between our ^^^^(1, . . . ,n) and the subamplitudes of Ref. [12] is ^'^^^(l, . . . ,n) = 
g^[N, {l,...,n) + Nf A^'-^^ ^'^'^ (1, . . . , n)] . 

To present the coUinear behaviour of the colour subamplitudes, we first explicitly relate 
the colour matrix Sp of the splitting process g — t- gig2 to its colour stripped component. 
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the splitting amplitude Split (see Eq. f lT^ ). At the tree level we define: 

W,a.;a)(^^^^^.p) ^ ^ ^ j^^^^^ SpU^'^p^ , p,; P) , (217) 

where the explicit expression of Split can be extracted by direct inspection of Eq. f|T5|) . 
Analogously, at the one-loop level we consider the splitting matrix Sp^^ (see Eqs. ( 128|) and 
and we define the splitting amplitude Split^^ as follows: 

Spi\%T"'''\p^^P^-^P) ^ ^ ' f^i^^a Split5;^(pi,p2;P) . (218) 

We note that, according to the notation used throughout this paper, the splitting ampli- 
tudes on the right-hand side of Eqs. (I217P and ( 1218^ are still spin matrices: they act onto 
the spin (helicity) space of the two coUinear gluons (with momenta pi and P2) and their 
parent gluon (with momentum P). 

We recall that the colour subamplitudes ^"^"^^(l, . . . , n) and ^'^^-'(l, . . . , n) in Eqs. (12151) 
and (I216P are colour ordered and, therefore, their kinematical structure depends on the 
specific ordering (1, . . . ,n) of the gluon momenta of the external legs. This features has 
consequences on the coUinear behaviour of the colour subamplitudes. Considering the 
coUinear limit of the momenta pi and p2, we can distinguish two types of configurations, 
according to whether the two coUinear gluons are adjacent or not adjacent in the argu- 
ment of the colour subamplitude. The colour subamplitudes A{.., 1, .., 2, ..), where the two 
coUinear momenta are not adjacent, are not singular in the coUinear limit. The colour 
subamplitudes ^(.., 1, 2, ..), where the two coUinear momenta are adjacent, are singular in 
the coUinear limit. 

At the tree level, the singular behaviour of the multigluon colour subamplitudes with 
adjacent coUinear legs is given by the following (helicity space) factorization formula: 

(. . . , fc, 1, 2, J, . . . ) ^ Split {p^,p2■, P) A^'^ (. . . , fc, P, . . . ) . (219) 

The colour subamplitude on the right-hand side has n — 1 external legs. It is obtained from 
the colour subamplitude on the left-hand side by replacing the two coUinear gluons with 
a single gluon (with momentum P). The relative ordering of the non-coUinear legs in the 
argument of the colour subamplitudes is left unchanged in going from the left-hand to the 
right-hand sides. We have explicitly introduced the labels k and j of two non-coUinear legs 
to remark the unchanged ordering of the non-coUinear legs. Note, however, that the tree- 
level splitting amplitude Split*^"^ is universal; it depends on the momenta (and helicities) 
of the two coUinear gluons and of the parent gluon, and it has no dependence on the non- 
coUinear legs. Moreover, the splitting amplitude Split '^"^ on the right-hand side of Eq. f l219p 
controls the singular behaviour of the coUinear splitting subprocess g — gig2 also in the 
case of colour subamplitudes with both gluons and quark-antiquark pairs in the external 
legs. 

At the one-loop level, the singular behaviour of the leading-colour subamplitudes with 
adjacent coUinear legs is given by the following g'enerafced factorization formula (in helicity 
space): 

^«(...,A;,l,2,j,...) ^ Sp\it^'\p^,p2;P) A^%..,k,P,j,..^) 

+ Split^'\p,,p,,p2,Pj;P) A^'\...,k,P,j,...) , (220) 
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and the one-loop splitting amplitude Split ^^-^ is 

Split«(Pfc,Pi,P2,P,;P) = Split^Hpi,p2;P) + /c7(Pfc,Pi,P2,P,;P) Split(°)(Pi,P2; ,(221) 

where the 'rational' part Split^'' is defined in Eq. (12181) and the transcendental function 
that multiplies Spht^°^ is 

Ic{Pk,Pi,P2,Pj;P) = glcr (^—^ji2~~^ 

fie; zi - zOsfci) + /(e; z^ - iOs.s) 1 1 • (222) 



1 1 

" ^ ~ 7 

The arguments (. . . , A;, 1, 2, j, . . . ) and {. . . , k, P, j, . . .) of the colour subamplitudes in the 
factorization formula (12201) are the same as those in Eq. (12191) . 

The tree-level factorization formula (12191) is a well-known result [16]. At the one-loop 
level, considering the TL collinear limit (i.e. su > 0), the momentum fractions zi and Z2 
are positive, and we can remove the iO prescriptions on the right-hand side of Eq. (I222p : 
therefore, Split^^-* is universal (i.e., independent of pk and pj), and Eqs. (I220l) - (l222p give 
the one- loop factorized result derived in Refs. fTIj [23] . The structure and the explicit form 
of Eqs. (I220p - (I222I) in the case of the SL collinear limit (i.e. Si2 < 0) is a new result, which 
derives from the application of the general results in Eqs. (133!) . (!34|) and (135|) (or fHSj) ) to 
pure multigluon scattering amplitudes at one-loop order. 

The main new feature of the SL collinear limit is that the one-loop splitting amplitude 
Split^^) in Eq. fl220D is not 'universal', since it also depends on the non-collinear partons. 
More precisely, Split*-^^ depends on the two non-collinear gluon legs k and j that are adjacent 
(colour connected) to the two collinear gluons in the one-loop leading-colour subamplitude 
^(i)(. l,2,j, ...). The dependence (see Eqs. (I22TD and (l222|) ) is due the signs of 
Ski = '2pk ■ Pi and Sj2 = 2pj ■ p2, which control the imaginary part of the analytic functions 
/(e; zi - iOski) and /(e; Z2 - iOs^a) in Eq. ([222D. 



B Appendix: The IR divergences of the two-loop split- 
ting matrix 

In this Appendix we illustrate the perturbative expansion of the IR factorization formula 
(I137p . We recall that Eq. (I137P refers to the collinear limit of m (m > 2) parton momenta 
in amplitudes with n [n > m) external QCD partons. 

The splitting matrix operators / and I in Eqs. (I139p and (I140p have the following 
renormalized perturbative expansions: 

/(e) = /(i)(e) + j(2) + ^ (223) 

7(e) = Y'\e)+t'\e) + ... . (224) 

The perturbative contributions I^'^^ and I^'^^ are obtained from the corresponding scattering 
amplitude operators, and I^, according to Eqs. (I135P and (I136p . To be precise, we 
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use the relations ( 11151) and (11161) at one-loop order {k = 1), and analogous relations at 
two- loop order {k = 2). Performing the collinear limit (as specified in Eq. (IllSp ) of the 
expression (11051) . we obtain 



j(i)( 



2n 2 

1 



i=i ^ ^ jeNC 



iO 



- - ^ T, ■ T, In 



-sg - iO 



jeNC 



(225) 
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where we have also used colour conservation (X^igc-^* ~ ~ XljgArc -^i )• Considering 
Eq. ( I105p with the replacement — t- A^, and applying the definition in Eq. ( I116p . we 
obtain 



j(l)/ X ^ «s(^^) 1 
27r 2 



1 



1 



- + - - - In 



-s-p - iO 



. (226) 



Since the two-loop operator if) (or I^) is simply related to if) (or J-jl^) by the expression 



(I119p . the operators J*-^-* and P^' are 
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(228) 



The computation of the perturbative expansion of Eq. ( 11371) is elementary and straight- 
forward. We simply illustrate few intermediate steps. To perform the expansion, we find 
it convenient to rewrite Eq. (I137P in the following equivalent form: 

Sp= [l-V(e) V-i(e)] 5p + V(e) V"'(e) . (229) 
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The two-loop perturbative expansion of the operator 1 — VV ^ is obtained by using 
Eqs. (II39D, (din]), (I223D and (1221; we have 



l-V(e) V-He) 



= J«(e)-T^'\e 



+ 0{al) . (230) 



Using Eqs. (11381) . (11401) and (12241) . the two- loop expansion of the second term on the right- 
hand side of Eq. f l229|) is 



V(e)5p'^"-V '(e 



(l)fin. 



+ Sp 



(2) fin. 



+ 



Sp 



(l)fin. 



{l + 0(ai)} . 

(231) 

Note that in Eq. f l23ip we have used the fact that Sp^^'^^ is strictly factorized and, hence, 

it commutes with I (e) (as already pointed out in Eqs. (11131) and (11141) ). Inserting the 
perturbative expansions (183!) . (I230p and (I23ip in Eq. (I229p . we reobtain the one- loop re- 
lations (!96|) and (!97|) . and we directly obtain the two-loop relations (11560 and (I157p . The 
three contributions on the right-hand side of Eq. (I157p depends on the perturbative terms 

I^^\e) and / (e) {k = 1,2). The one-loop operator I^^^{e) is given in Eq. (11171) . while 



I^^U^) and Sp 



(2) div 



are given by the following expressions: 
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We discuss the structure of Eqs. (I232p and (I233p . in turn. 

Using elementary algebra, the expression (12320 can be rewritten as follows: 

1 
2 



/(^)(e) 
/«(el 



f{i) 



+ A 



(2; 2), 

rnC ^ 



where we have defined 



The form on the right-hand side of Eq. (I234p is useful to express the operator / 



terms of the one- loop operator I^^^ = /'■^■' — I^^^ (see Eq. (11171) ). Owing to the relations 

(12271) and fl228|) . the linear combinations 2 /^^^ + [J^^^j^ and 2 7^^^ + [T^^^]^ are proportional 

to I^^^ and /'■^^ respectively. Therefore, by simple inspection of Eqs. (I227p and (I228p . we 
see that Eq. (I234p corresponds to the result reported in Eq. (I158p . 

The two-loop factorization breaking operator A,^'^ in Eq. (I158p originates from the 
colour-matrix commutator in Eq. (12350 (the commutator trivially vanishes in QED and 
any abelian theories). The commutator can be rewritten as follows: 



(2) 
mC 
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Since the operator J*'^'' does not depend on the colour matrices of the coUinear partons, 
to obtain Eq. (12361) we have exploited the fact that only A^^^(e) (the factorization break- 
ing part of I^j^c(^) -^'l- dsn])) contributes to the commutator. Inserting the explicit 
expressions (llOOp and (I226p in Eq. (12361) and computing the colour-charge commutator, we 
straightforwardly obtain the result in Eq. (I159p . 

In the TL collinear region, the one-loop splitting matrix Sp^^''^^ and its IR finite part, 
5'p(i)fin.^ in Eq. ( l96i) are strictly factorized, and they do not depend on the colour matrices 
of the non-coUinear partons. Therefore, the commutator in Eq. (I233p vanishes in the TL 
collinear limit. In the SL collinear region, Sp^^'^^^' depends on the colour matrices of the 
non-coUinear partons (as mentioned in Sects. 15.11 and 15.31 Sp^^'^ has a linear dependence 

on these colour matrices) and, therefore, the commutator term Sp ^ in Eq. (I233p does 

not vanish. Since Sp^^''^^' is IR finite, the IR divergent part of Sp^'^^'^"' contains only single 

poles 1/e (the coefficient of the double pole in the expression (I226P of I^^^ (e) is a c- number 
and, hence, it gives a vanishing contribution to the commutator in Eq. (I233P ). 

In the specific case of the SL collinear limit of m = 2 partons, we explicitly know 
Sp^^'^^ to all orders in e (see Sect. 14. ip . This information can be exploited to extract 
5'p(i)fin- g^nd^ then, to explicitly compute Sp ^ in Eq. (I233p . The multiparton collinear 
operator I^^^q {k = 1,2) is denoted by I^2c ^^^^ of m = 2 collinear partons. From 

the expressions in Eqs. ( p9|) and (llOOp . we have 

It: I yt t 
+ - (Clin Cain |z2|)--Ti-T2 lnf ~"7~'° 

27r " 1 
+ z — ^ J] ■ T, e(-zO sign(s,,) \ , (237) 

^ i=3 i=l,2 J 

whereas I^2c obtained from Eq. f ll58p by the replacements I^^q — )■ l''2c ^^'^ ^^mc ~^ 
^fc^. Using Eq. ^ and, then, Eqs. (El and OHZD, we have 

5p(i)«"- = - 5^(1)^^^- = fc\^) Sp^''""^ + Sp^i'""^ - 4^i(e) 5p(°'^) . (238) 



Sp 



Inserting Eq. (I238p in Eq. (I233p . we obtain 

Y'\e) , 7?(e) - 7« (e) ] ^ A^^^'\e) , (239) 

where we have used the fact that I^^^ (e) commutes with both Sp^^'^^ and Sp^^'^^ , and we 
have defined the colour operator A^2C^ ■ Using the explicit expressions of Iq and I^^, 
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in Eqs. 



and fl237p . we evaluate the commutator in Eq. 



and we find 



I (e) , Ic (e) 
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Therefore, in the specific case of the two-parton collinear limit, the computation of Sp^'^^'^"' 
(see Eqs. f ll56p and f ll57p ) is explicitly completed in the form 
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where /111, J*-^^ 
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2(^, J.2C cind given in Eqs. f l237p . (11580 and (12400 . In particular, the 

('2,' 2^ fS' ______ 

factorization breaking operators A2 and A2 in Eqs. (I159p and (I240p can be combined 
by defining 
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and we obtain 
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We note that the expression of Sp*-^'^-* presented in Sect. 16.21 (see Eq. f ll44p ) has a form 
that differs from the expression in Eq. (12411) . However, these two expressions are completely 
equivalent since they lead to the same IR divergent contribution to Sp^"^'^^ (the differences 

~{2)fin 

in the IR finite part can be absorbed in the definition of Sp ). To be precise, comparing 
Eqs. dm]) and we have 

5p(2,«)_5p(2)div. ^ (1^^ {e) - I^Ue)) Sp^^'^^ 

+ (?c (e) - If ) (e) - Af^^)(e)) + Sp'^'"''- , (244) 

and it can be explicitly checked that the expression on the right-hand side is IR finite if 
e — )■ 0. This explicit check, which is left to the reader, can be performed by using the 
expressions in Eqs. dUSD, (^M>, dlSSD, (11521), (ESZD and flMID . 

We briefly illustrate the derivation of the relations in Eqs. (I162p and (I163p . 

We first note that the explicit constraint sjk < can be removed from the expression on 
the left-hand side of Eq. (11620 . Indeed, the sum of the terms with sjk > gives a vanishing 
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contribution to that expression, as explicitly shown by the following relation: 
fabc ^ T^T^ sign(sij) Q{sjk) = ^ fabc ^ T^T^ [ sign(sij) - sign(sa.) ] Q{sjk) = . 

j, k j, k 

j^k jy^k 

(245) 

Here, we have simply replaced the factor sign(sjj) Q{sjk) with its antisymmetric part with 
respect to the exchange j ^ k (the symmetric part gives a vanishing contribution, since 
fabcTjTj: is antisymmetric under the exchange j O k) and, then, we have used the fact 
that Sij and Sik have the same sign if sjk > 0. Considering the left-hand side of Eq. f ll62p 
and removing the constraint sjk < 0, we thus obtain 



fabc E T^^TJTk sign(s,,) e{-sj,) K = fabc E sign(s,,) h, 

ieC j,keNC iGC j,k&NC 

j^k j^k 



i&C j&NC l&C 



where we have performed the sum over k by using the colour conservation relation 

Y Tl = -YTt-T^, . (247) 



k^NC teC 
k+3 



Owing to the algebraic identity fll50p , the final expression in Eq. fl246p exactly corresponds 
to the expression on the right-hand side of Eq. 



The left-hand side of Eq. f ll63p can be written as 

- fabc J2 E T^^T^T^^ignis^p) Qi-z,) , (248) 

i,ieC jeNC 

where we have simply used the collinear approximation Sij ~ Zi s-p. The only difference 
between the right-hand side of Eq. fll63p and the expression fl248p is due to the presence 
of the explicit constraint su < 0; however, this apparent difference is harmless, since the 
sum of the terms with su > gives a vanishing contribution to the expression f l248p . The 
vanishing of this sum follows from the relation 

fabc Yl T^'T^Qi~z,)Qis,e) = ^fabc Y [ e(-^,) - e(-^o ] e(s,,) = o . (249) 

Here (analogously to the procedure used in Eq. fl245p ). we have simply replaced the factor 
Q{—Zi) Q{sii) with its antisymmetric part with respect to the exchange i -{^ I and, then, 
we have used the fact that 0(— Zj) = Q{—zi) if > (i.e., ZiZi > 0). 



C Appendix: TL collinear limit, strict factorization 
and requirement of IR consistency 

In this Appendix we consider the TL collinear region, and we discuss how a strictly- 
factorized splitting matrix can be recovered from the all-order IR structure presented in 
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Eq. fll37p . Then, we present explicit expressions for the IR structure of the TL sphtting 
matrix. 

In the multiparton TL coUinear hmit, the all-order splitting matrix Sp is strictly factor- 
ized (see Eq. ( 18T|) ) and, thus, it does not depend on the colour matrices of the non-coUinear 
partons. On the contrary, owing to Eq. (11361) . the matrix structure of V only depends on 
the colour matrices of the non-collinear partons. Therefore, Sp commutes with V and, 
inverting Eq. (I137p . we obtain 

= V~\e) Sp V(e) = V-\e) V(e) Sp . (250) 

This relation can be inverted to give 

Sp = V'\e) V(e) Sp^""- . (251) 

Note that, in the derivation of Eqs. (I137P and (I25ip . we have not used the property that 
Sp^'^' is strictly factorized (in general, the assumption that Sp^'^' is strictly factorized 
is not valid, even in the TL case). Nonetheless Sp^'^' is IR finite and, therefore, the IR 
divergent contributions to Sp are produced by the operator V V on the right-hand side 
of Eq. (I25ip . Since Sp and, hence, its IR divergent terms are strictly factorized, Eq. (125 ip 
enforces a constraint on V V: the IR divergent part of the operator V V has to be 
strictly factorized. More precisely, this operator must have the following form: 

V"'(e) V(e) = VTL(e) Vfin.(e) , (TL coll. lim.) , (252) 

where the IR divergent operator VxL(e) is strictly factorized, whereas the operator Vgn (e) 
is IR finite (but it is not necessarily strictly factorized) and its form is such that the matrix 
Vfin.(e) Sp^""^' is strictly factorized. This constrained structure of Eq. f l252p guarantees that 
the right-hand side of Eq. (125 ip and, hence, the splitting matrix Sp are strictly factorized. 
Note that the IR divergent operators V and V separately depend on the momenta and 
colour charges of the non-collinear partons. This separate dependence is certainly con- 
strained by Eq. (I252p (in particular, the dependence largely cancels in Eq. (I252p ). since the 
IR divergent operator Vtl is strictly factorized and, thus, it is completely independent of 
the non-collinear partons. 

The constrained structure of Eq. (I252p has been derived by using two properties of mul- 
tiparton QCD scattering amplitudes: their strict factorization in the multiple TL collinear 
limit of m partons (Eqs. f lTHj) and flHTj) ) and their IR structure (according to Eq. fll24p or, 
equivalently, Eq. (I127P ). This constrained structure can be regarded as a requirement of 
consistency between the TL collinear limit and the IR properties of the QCD amplitudes. 
If at some perturbative order Eq. (I252p is not valid, either strict TL collinear factorization 
is violated or the IR structure in Eqs. (I124p and (I127P is not valid at the corresponding 
perturbative order. 

In the case of m = 2 collinear partons, the presence of a valuable IR consistency 
constraint from the behaviour of the QCD amplitudes in the TL collinear limit was pointed 
out in Ref. [7] and exploited also in Ref. [8]. The discussion in this subsection on the TL 
collinear limit generalizes the discussions in Sect. 5 of Ref. [7] and Sect. 4 of Ref. [8]; our 
generalization deals with the extension to m (m > 3) collinear partons and to a generic 
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all-order form of the IR operator Yni^) (or lAf(e)) of the scattering amplitudes. Note that 
both the operators V and V in Eq. (12521) originate from VM(e) in Eq. (11271) (the form of 
\jY follows from Vjvf by simply reducing the total number of the external parton legs in 
A4) through the coUinear-limit procedure in Eqs. (I135P and (I136p . Therefore, the strict 
factorization requirement in Eq. (I252p eventually constrains the colour and kinematical 
structure of the IR operator VM(e) at arbitrarily-high perturbative orders [71 

To sharpen our all-order discussion of the multiparton TL collinear limit, we consider 
the case in which the IR operator VA/(e) has a 'minimal form', which includes only the 
terms proportional to the IR poles 1/e'' with k > 1, whereas additional terms of order 
e°, e, and so forth are absent (we do not specify the colour and kinematical dependence 
of Vm)- The corresponding exponentiated operator in Eq. (11291) is denoted by /Af,cor(min) 
and, analogously, the exponentiated splitting matrix operators in Eqs. (I139p and (I140p are 
denoted by /cor(min) and /cor(mm), respectively. We also define the following operator: 

-fTL,cor(e) = -'^cor(min) (e) — -fcor(min) (c) ; (253) 

which, by definition, also has a minimal form. Having specified these definitions, we sharply 
state our main conclusion: the strict factorization of the splitting matrix Sp is equivalent to 
the requirement of strict factorization of both the collinear matrix Sp^^' and the operator 
-'^TL.cor- The proof of this statement is given below. 

If Sp is strictly factorized, the strict factorization of Sp^'^' is a simple consequence of 
the minimal form of Ym- Indeed, owing to the relations (I135P and (I136p . if VA/(e) has a 
minimal form, V and V also have a minimal form: this implies that the IR finite operator 
Vfin. in Eq. (I252p is trivial (i.e., Vfin.(e) = 1) and, therefore, that the IR finite matrix Sp^^' 
in Eq. (12510 is strictly factorized. The strict factorization of Sp^"^' then implies that Sp^^' 

commutes with V and, using Eq. (I137p . we get 

5p = V(e) V"'(e) Sp^""- . (254) 

Comparing Eqs. (125 ip and (I254p . we conclude that the operators V and V ^ commute. 
This commutation property is valid also for the corresponding exponentiated operators 
-fcor(min) (e) and -rcor{min) (e) in Eq. (I253p and, therefore, the operator Vtl in Eq. (I252p has 
the following minimal form: 



V 



TL(e) = V ^(e) V(e) = expj /cor(mm)(e) - i"cor(min)(e)} 

= exp|/TL,cor(e)} , (TL coll. lim.) . (255) 

Here, we have used the property e^e^ = e^~^^, which is valid if A and B are commuting 
matrices. Since Vtl is strictly factorized, Eq. (I255P finally implies that the exponentiated 
operator /TL,cor(e) of Eq. (I253p is strictly factorized. 

If Sp^^' and the operator /TL,cor(e) in Eq. (I253P are strictly factorized, we can easily 
proof that Vtl and, hence, Sp are strictly factorized. Indeed, if /TL,cor in Eq. (12530 is 
strictly factorized, this operator commutes with /cor(min)- Therefore, we have 

[-^TL,cor(^) ? -^cor(min) (^) ] [ -^cor(min) (^) -^cor(min) (^) ) -^cor(min) (^) ] 

~ [ -^cor(min)(f) > -^cor(min) (^ ) ] ; (256) 
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namely, /cor(min) and Jcor(min) are commuting operators. This commutation property leads 
to the explicit expression f l255p of the operator Vtl, and this expression is evidently strictly 
factorized. 

In summary, the validity of strict factorization in the TL collinear limit requires that 
the IR operator /TL.cor in Eq- (12531) is independent of the non-collinear partons. Since 
/cor(min) (e) and -rcor{min) (c) dcrivc from the collinear limit (see Eqs. (I135p . (I136p . (I139p and 
(I140p ) of the scattering amplitude operator Ym{^) (see Eq. f ll27p ). this requirement directly 
constrains the colour and kinematical structure of the exponentiated operator /Af,cor(min) 
(see Eq. (I129P ) in its minimal form [TJ [8]. 

Moreover, using Eqs. (I25ip . (I252p and (I255p . the IR structure of the splitting matrix Sp 
for the multiparton TL collinear limit can be presented in the all-order form of Eqs. (11411 - 

In the derivation of Eqs. (II 41 p - ( 11431) . we have used a minimal form of the scattering 
amplitude operator Va/ (e) in Eqs. (ll27p -f lT29|) . This minimal form is explicitly presented in 
Eqs. (fT05|) and (1TT9|) at one-loop and two-loop orders, respectively (we recall that I^j^j ^^j. — 

I^^j and 2lf) = 21^'^) + (/^j )^ ). Therefore, using Eq. (]253p . we can explicitly compute 
the perturbative expansion of /xL.cor, Vtl and Itl up to the two-loop order. To be precise, 
we define the renormalized loop expansion of Itl: 

lTL(e) = lS(e) + lS(6) + ... , (257) 

and we straightforwardly find 

+ ^lS(2e) 

(258) 

where the explicit expression of the one-loop operator Itl(c) exactly corresponds to the 
expression (llOip of the IR operator I^^(j{e) in the TL collinear region. The one- loop 
perturbative expansion of the TL factorization formula (I142p gives Eqs. (l96il and ( 1971) with 
the obvious identification of li^L with The two-loop perturbative expansion gives 

= iW(g) +l5?i(e) + 5p(2)fi- ^ coll. lim.) , (259) 

where iTL(e) is explicitly presented in Eq. ( I258p . In the case of m = 2 collinear partons, 
the results in Eqs. (I14ip - (ll43p and Eqs. (I257p -( l259l) exactly correponds to those presented 
in Refs. 

In this Appendix we have explicitly derived Eqs. (I14ip - (ll43p by considering IR diver- 
gent operators that have a minimal form. However, Eq. (I14ip has the same IR factorization 
invariance as Eq. ( ]127p . Specifically, the right-hand side of Eq. ( I14ip is invariant under the 
joint transformations (redefinitions) VTL(e) — >■ VTL(e) Vfin.(e) and Sp^^' -> Vfr^^(e) Sp^^' , 
where Vgn. (e) is an invertible IR finite operator that is strictly factorized. Since this invari- 
ance can be used to redefine Vtl and Sp^^' , the structure of Eqs. (I14ip - (ll43p is actually 
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valid independently of the minimal form of the IR divergent operator Vtl (or, Itl and 

-fxL.cor)- 
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Note added 



After the completion of the present paper, J. Forshaw, M. Seymour and A. Siodmok 
made an important observation [66]. We would like to thank Jeff, Mike and Andrzej for 
discussions and the communication of their result before its publication. The observation 
regards the expectation value of the two-loop operator Ap (e) (see Eq. fll84p ). which gives 
the (IR dominant) factorization breaking contribution to the squared splitting matrix p(2.^) 
(see Eq. 01811) ) for the two-parton coUinear limit. The expectation value onto the reduced 
matrix element A4 is 

{M\ Ag)(6) 5p(°'«) \M) = (AT^^''^^ I Ag)(e) \M^'''') 



+ 



(^Sp^o,R)y ^(2)(^) Spio,R) 17^(0'^)) + e.c. 



higher orders , (260) 



where the right-hand side corresponds to the perturbative (loop) expansion of Ai {Ai = 

-i-T(0,R) -ir-ri^'R) \ 

M +M +...). 

The key observation [BB] is that, at the lowest-order level (i.e., considering the expec- 
tation value onto A4 ' , as given in the first term on the right-hand side of Eq. (12601) ) 
the expectation value vanishes in pure QCD (i.e., if the lowest-order reduced matrix el- 
ement Ai is produced by tree- level QCD interactions). We refer to Ref. [66] for the 
explanation and discussion of this effect. 

Note that Ap is not vanishing. The lowest-order vanishing of the expectation value 
in Eq. (I260p can indeed be avoided by changing (and properly choosing) the lowest-order 

matrix element A4^^'^\ For instance, we can consider tree- level quark-quark scattering 
produced by electroweak interactions (with CP- violating electroweak couplings and/or fi- 
nite width of the Z and bosons), or we can supplement tree- level QCD scattering with 
one-loop (pure) QED radiative corrections. Therefore, as a matter of principle, it remains 
true that the operator Ap explicitly uncovers two-loop QCD effects that lead to violation 
of strict collinear factorization at the squared amplitude level. In particular, the conceptual 
discussion presented in Sect. 17.41 (and briefly recalled at the end of Sects. [1] and [8]) continues 
to be valid, although the lowest-order partonic subprocess ' parton -|- parton — )■ 2 partons ' 
(see Eq. f l21ip ) that is used as starting point of the discussion has to be interpreted in the 
generalized sense mentioned at the beginning of this paragraph. Of course, the fact that 
the lowest-order partonic subprocess in Sect. 17.41 is not due to tree-level QCD interactions 
reduces the possible phenomenological consequences of the effects discussed in that section. 

In the context of pure QCD, the expectation value in Eq. fl260p is not vanishing at 
the next-to-lowest order (i.e., the second term on the right-hand side of Eq. f l260p ). which 

is obtained by the one- loop QCD correction A4 ' to the reduced matrix element A4. 
Therefore, Ap'' certainly contributes, through the two-loop splitting matrix p(2'^)^ to the 
SL collinear limit of three-loop QCD squared amplitudes. At the three-loop level, the 
violation of strict collinear factorization produced by p(^'-^) (e.g., Ap^) joins additional 
factorization breaking effects that are produced by the three-loop collinear matrix p(3.^). 

This combined occurence of p(^'^) and p(^'^) in the collinear limit of QCD squared 
amplitudes at the three-loop level has prompted our preliminary investigation of three-loop 
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effects. The analysis and results of the present paper can be readily and straightforwaldy 
extended to higher-loop orders by exploiting the exponentiated structure of the leading 
(i.e., 0{a'^/e"'~^^)) and next-to-leading (i.e., C(ag/e"')) IR divergences of multiparton QCD 
scattering amplitudes (see Sect. 16.11 and Refs. [50l [H [2])- In the case of the SL collinear 
limit of two partons, we know the complete e dependence of the one-loop splitting matrix 
Sp^^'^^ and, therefore, we have explicit control of three-loop factorization breaking effects 
starting from their dominant IR divergent terms of (9(ag/e^). Computing P|^^^^^ (the part 
of p(3'-^) that violates strict collinear factorization) up to the accuracy of (9(ag/e^), we 
obtain 



p(3..) ^ ^^^(0,.) t 1 1^^, (e) - ^ i b,^ Ag) (el 



+ 



1 

6 



r^'(e)-/e(e) - h.c. , AV'(6) 



5p(o,iJ) + o(^^^ , (261) 



where 1^ , Ip^ and I^^^ are the customary one-loop operators used throughout the paper 
(see Eqs. (ES]), fll79p and f l226p ) and bo is the first-order coefficient of the QCD /5 function 
(see Eq. fl52]) ). The three- loop expression in Eq. (126 ip is controlled by the iterated action 
of Ap without (first term on the right-hand side of Eq. (I26ip ) and with (second term on 
the right-hand side of Eq. (126 ip ) an additional colour commutator with one-loop terms. 

We note that the three- loop factorization breaking contribution from p(2'^) (the second 
term on the right-hand side of Eq. (I260p ) cannot be cancelled by the factorization breaking 
contribution from P^^ (e.g., the lowest-order expectation value of the second term on the 
right-hand side of Eq. (126 ip ). In particular, we remark that in pure QCD each of these 
two terms produces an IR divergent three-loop contribution of order ctg/e^. In the case of 
scattering amplitudes with n = 5 QCD partons, the colour correlation structure of these 
two non-cancelled factorization breaking contributions is analogous to the commutator 
structures that were found [33j in the N^LO computation of super-leading logarithms. 
This observation is consistent with the common physical mechanism that originates the 
violation of strict collinear factorization and the emergence of super-leading logarithms in 
'gaps-between-jets' cross sections (as mentioned at the end of Sect. 17.41 and remarked in 
Ref. 
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